Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 17,340 T→C intergenic (‑382/‑148) mokC ← / → nhaA regulatory protein for HokC, overlaps CDS of hokC/pH‑dependent sodium/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL60617,3400TC98.3% 241.8 / ‑3.8 60intergenic (‑382/‑148)mokC/nhaAregulatory protein for HokC, overlaps CDS of hokC/pH‑dependent sodium/proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (33/26);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

TTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACA  >  REL606/17267‑17409
                                                                         |                                                                     
ttaatCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCa                                                                     >  1:916880/1‑76 (MQ=255)
 taatCTCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCaa                                                                    >  1:192101/1‑76 (MQ=255)
  aatCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCaaa                                                                   >  1:1046302/1‑76 (MQ=255)
  aatCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCaaa                                                                   >  1:1281944/1‑76 (MQ=255)
   atCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCaaaa                                                                  <  2:1009557/76‑1 (MQ=255)
      ccTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAAc                                                               >  1:749593/1‑76 (MQ=255)
       cTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACt                                                              <  1:1233436/76‑1 (MQ=255)
       cTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACt                                                              <  2:1303281/76‑1 (MQ=255)
       cTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACt                                                              >  1:459631/1‑76 (MQ=255)
       cTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACt                                                              <  2:1318582/76‑1 (MQ=255)
        tCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTc                                                             <  1:60290/76‑1 (MQ=255)
         cAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCa                                                            >  2:1153219/1‑76 (MQ=255)
             aaTCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCAtctc                                                        >  1:1055846/1‑76 (MQ=255)
                ccTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCAcc                                                     >  1:831793/1‑76 (MQ=255)
                ccTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCAcc                                                     >  1:13692/1‑76 (MQ=255)
                ccTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCAcc                                                     >  1:1321858/1‑76 (MQ=255)
                 cTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCg                                                    <  2:1405924/76‑1 (MQ=255)
                 cTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCg                                                    <  1:586547/76‑1 (MQ=255)
                  tCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGc                                                   >  1:934429/1‑76 (MQ=255)
                    gccTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTg                                                 <  1:1622367/75‑1 (MQ=255)
                    aCCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTg                                                 >  1:1517741/1‑76 (MQ=255)
                     ccTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGa                                                >  1:734166/1‑76 (MQ=255)
                        tAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATgg                                             >  2:294160/1‑76 (MQ=255)
                        tAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATgg                                             >  1:1721991/1‑76 (MQ=255)
                         aaGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGgc                                            >  2:1655136/1‑76 (MQ=255)
                           gCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGgcgc                                          <  1:1022575/76‑1 (MQ=255)
                            cTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCa                                         <  1:61822/76‑1 (MQ=255)
                            cTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCa                                         <  1:1092416/76‑1 (MQ=255)
                            cTATGATTATCTGGGCCTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCa                                         >  2:1665882/1‑76 (MQ=255)
                                 aTTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAAttc                                    >  2:1304999/1‑76 (MQ=255)
                                    aTCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAAttcttc                                 <  2:149642/76‑1 (MQ=255)
                                      cTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCaa                               <  2:1442741/76‑1 (MQ=255)
                                        gggCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATg                             >  1:533352/1‑76 (MQ=255)
                                        gggCTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATg                             <  2:383042/76‑1 (MQ=255)
                                           cTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGAt                          <  2:912362/76‑1 (MQ=255)
                                           cTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGAt                          <  2:876653/76‑1 (MQ=255)
                                           cTTAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGAt                          <  2:61049/76‑1 (MQ=255)
                                             tAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCg                        >  2:112936/1‑76 (MQ=255)
                                             tAGGGTCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCg                        >  2:804608/1‑76 (MQ=255)
                                                  tCACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTaaaa                   >  1:1681507/1‑76 (MQ=255)
                                                   cACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTaaaaa                  >  2:1766896/1‑76 (MQ=255)
                                                   cACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTaaaaa                  <  1:122807/76‑1 (MQ=255)
                                                   cACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTaaaaa                  <  1:1559842/76‑1 (MQ=255)
                                                   cACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTaaaaa                  >  1:1597037/1‑76 (MQ=255)
                                                   cACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTaaaaa                  <  1:30783/76‑1 (MQ=255)
                                                   cACTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTaaaaa                  <  1:469769/76‑1 (MQ=255)
                                                     cTCGTGAGCGCTTACAGCCGCCAGAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAAc                >  2:1162947/1‑76 (MQ=255)
                                                     cTCGTGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAAc                >  2:26676/1‑76 (MQ=255)
                                                         tGAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGt            >  2:776716/1‑76 (MQ=255)
                                                          gAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAgat           >  1:621320/1‑74 (MQ=255)
                                                          gAGCGCTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGtt           >  1:1246151/1‑76 (MQ=255)
                                                             cgcTTACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTAtt        <  1:758176/76‑1 (MQ=255)
                                                                ttACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCt     >  1:210393/1‑76 (MQ=255)
                                                                ttACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCt     <  2:189894/76‑1 (MQ=255)
                                                                ttACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCt     >  1:1522394/1‑76 (MQ=255)
                                                                 tACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTa    <  2:1083775/76‑1 (MQ=255)
                                                                 tACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTa    <  2:797273/76‑1 (MQ=255)
                                                                 tACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTa    <  2:424551/76‑1 (MQ=255)
                                                                 tACAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTa    <  2:40438/76‑1 (MQ=255)
                                                                  aCAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTac   >  1:954141/1‑76 (MQ=255)
                                                                   cAGCCGCCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTaca  >  1:294630/1‑76 (MQ=255)
                                                                         |                                                                     
TTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACA  >  REL606/17267‑17409

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: