Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 433,668 C→T R55C (CGT→TGT)  hupB → HU, DNA‑binding transcriptional regulator, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606433,6680CT98.3% 205.7 / ‑3.7 59R55C (CGT→TGT) hupBHU, DNA‑binding transcriptional regulator, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/28);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

GCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAG  >  REL606/433570‑433764
                                                                                                  |                                                                                                
gCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAgtgt                                                                                                >  1:975684/1‑101 (MQ=255)
   gTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgct                                                                                             >  1:403694/1‑101 (MQ=255)
   gTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgct                                                                                             >  2:1376081/1‑101 (MQ=255)
   gTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgct                                                                                             <  2:1453656/101‑1 (MQ=255)
    tGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgctg                                                                                            >  2:491544/1‑101 (MQ=255)
      cGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGcc                                                                                          <  2:594396/101‑1 (MQ=255)
        ttAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCg                                                                                        >  1:1409063/1‑101 (MQ=255)
         tAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGCTAGGATTTGGTCCTTTTGCCGTTAAAGAGTGTGGTGcccct                                                                                       >  2:1224166/1‑99 (MQ=255)
           gATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTAc                                                                                     <  2:631959/101‑1 (MQ=255)
            aTGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACt                                                                                    >  2:63558/1‑101 (MQ=255)
             tGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTg                                                                                   >  1:1516729/1‑101 (MQ=255)
               cTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGc                                                                                 >  2:620182/1‑101 (MQ=255)
               cTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGc                                                                                 <  2:206367/101‑1 (MQ=255)
                tattatTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGATACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGcc                                                                                <  1:942116/101‑1 (MQ=255)
                     ttGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAc                                                                           <  1:1705013/101‑1 (MQ=255)
                        cTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCg                                                                        >  1:5984/1‑101 (MQ=255)
                           ccGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAg                                                                     <  2:143548/101‑1 (MQ=255)
                                aCTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCgg                                                                >  1:693664/1‑101 (MQ=255)
                                  tGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTa                                                              >  1:1200330/1‑101 (MQ=255)
                                   gAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTaa                                                             >  1:1581580/1‑101 (MQ=255)
                                   gAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTaa                                                             >  2:792700/1‑101 (MQ=255)
                                     aTCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAag                                                           <  1:1365121/101‑1 (MQ=255)
                                       ctctGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTagagag                                                         >  2:1080850/1‑101 (MQ=255)
                                       ctctGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAagag                                                         >  2:293333/1‑101 (MQ=255)
                                       ctctGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAagag                                                         >  2:145985/1‑101 (MQ=255)
                                       ctctGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAagag                                                         <  2:119825/101‑1 (MQ=255)
                                         ctGAAAGAAGGGGATGATGTAGCACTGTTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAt                                                       <  1:1646166/101‑1 (MQ=255)
                                          tGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATc                                                      <  2:809681/101‑1 (MQ=255)
                                           gAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCa                                                     >  2:174593/1‑101 (MQ=255)
                                           gAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCa                                                     <  1:901887/101‑1 (MQ=255)
                                           gAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCa                                                     >  1:1614503/1‑101 (MQ=255)
                                           gAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCa                                                     <  1:179525/101‑1 (MQ=255)
                                            aaaGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCAc                                                    <  1:903891/101‑1 (MQ=255)
                                              agaagGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCa                                                  >  1:487591/1‑101 (MQ=255)
                                                aagGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATc                                                >  1:851796/1‑101 (MQ=255)
                                                     gatgatGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCgctgc                                           <  1:1045636/101‑1 (MQ=255)
                                                          tgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTa                                      >  1:912393/1‑101 (MQ=255)
                                                          tgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTa                                      <  2:597169/101‑1 (MQ=255)
                                                              gCACTGGTAGGTTTTGGTACTTTTGTCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGt                                  <  2:1102776/101‑1 (MQ=255)
                                                              gCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGt                                  <  1:928218/101‑1 (MQ=255)
                                                              gCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTG‑CCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTa                                 >  2:1358132/1‑101 (MQ=255)
                                                                 cTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTAcc                               <  1:195598/101‑1 (MQ=255)
                                                                  tGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCg                              >  2:1175982/1‑101 (MQ=255)
                                                                   ggTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGa                             >  1:283477/1‑101 (MQ=255)
                                                                    gTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAg                            >  1:1010733/1‑101 (MQ=255)
                                                                         ttttGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTcc                       <  2:30040/101‑1 (MQ=255)
                                                                           ttGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGt                     <  2:975679/101‑1 (MQ=255)
                                                                            tGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTg                    >  1:1589300/1‑101 (MQ=255)
                                                                              gTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCa                  >  1:1646545/1‑101 (MQ=255)
                                                                               tACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAg                 <  1:620187/101‑1 (MQ=255)
                                                                                aCTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAgg                >  1:678148/1‑101 (MQ=255)
                                                                                aCTTTTGCCGTTAAAGAGGCTGCTTCCCGTACTGGCCGCAACCCGCAGACCGGCAAGGCGATCACCATCGCTGCTGCTAAAGCACGGAGCTTCCGTGCAgg                >  2:265837/1‑101 (MQ=37)
                                                                                   tttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTaa             <  2:151773/101‑1 (MQ=255)
                                                                                    ttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTaaa            <  1:1009959/101‑1 (MQ=255)
                                                                                    ttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTaaa            <  1:491549/101‑1 (MQ=255)
                                                                                     tGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAg           <  2:542238/101‑1 (MQ=255)
                                                                                     tGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAg           >  2:892793/1‑101 (MQ=255)
                                                                                     tGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGAAAAGAGATCACCATCGCTGCTGCTAAAGTCCCGACCTTCCGTGCAGGTaatg           <  2:403689/101‑3 (MQ=38)
                                                                                        cGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCAc        >  1:816299/1‑101 (MQ=255)
                                                                                         gTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGGAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGTGAAAGCACt       >  2:53083/1‑101 (MQ=38)
                                                                                           tAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGa     <  1:63558/101‑1 (MQ=255)
                                                                                            aaaGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGaa    >  1:168466/1‑101 (MQ=255)
                                                                                             aaGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGaaa   >  2:1488631/1‑101 (MQ=255)
                                                                                              agagTGTGCTGCCGGTACTGCCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAg  <  1:1349896/101‑1 (MQ=255)
                                                                                              agagTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAg  <  1:18519/101‑1 (MQ=255)
                                                                                                  |                                                                                                
GCCGTGCGTTAGATGCTATTATTGCTTCCGTAACTGAATCTCTGAAAGAAGGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCACCATCGCTGCTGCTAAAGTACCGAGCTTCCGTGCAGGTAAAGCACTGAAAG  >  REL606/433570‑433764

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑