Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,964,425 G→C 17.3% A209G (GCC→GGC)  flhA ← putative flagellar export pore protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,964,4250GC17.3% 141.9 / 26.6 58A209G (GCC→GGC) flhAputative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base C (6/4);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

TAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCA  >  NC_000913/1964326‑1964523
                                                                                                   |                                                                                                  
taACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATgg                                                                                                    <  1:102301/100‑1 (MQ=255)
taACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATgg                                                                                                    <  1:1288737/100‑1 (MQ=255)
 aaCTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGc                                                                                                   >  2:1824820/1‑100 (MQ=255)
 aaCTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGc                                                                                                   <  1:274102/100‑1 (MQ=255)
  aCTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCa                                                                                                  <  1:576598/100‑1 (MQ=255)
   cTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCAt                                                                                                 >  2:1419544/1‑100 (MQ=255)
       ccgccgCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATcgcc                                                                                             <  1:756343/100‑1 (MQ=255)
        cgccgcGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATcgccg                                                                                            >  2:1908755/1‑100 (MQ=255)
         gccgcGTGTCCCATTCTCATGCCATGGTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATcgccgg                                                                                           <  2:1716320/100‑2 (MQ=255)
         gccgcGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATcgccgg                                                                                           <  1:2156223/100‑2 (MQ=255)
          ccgcGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCcgcg                                                                                          <  1:2089776/100‑1 (MQ=255)
          ccgcGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGCCATcgccggg                                                                                          <  1:494231/100‑3 (MQ=255)
          ccgcGTGTCCCATGCGCAGGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCcgcg                                                                                          <  1:1292432/100‑1 (MQ=255)
             cgTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAAc                                                                                       <  2:2176033/100‑1 (MQ=255)
                 tCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAac                                                                                   <  1:1957595/100‑1 (MQ=255)
                  cccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAact                                                                                  <  1:2018487/100‑1 (MQ=255)
                  cccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAact                                                                                  >  1:2303174/1‑100 (MQ=255)
                  cccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAac                                                                                   >  1:1661538/1‑99 (MQ=255)
                   ccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAactt                                                                                 >  2:2219212/1‑100 (MQ=255)
                   ccATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAactt                                                                                 <  2:734001/100‑1 (MQ=255)
                    cATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAactta                                                                                >  1:1493483/1‑100 (MQ=255)
                       gCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAacttactt                                                                             >  2:64856/1‑100 (MQ=255)
                          cATGCCATGTTGCAGCACGCCGACCAGCAACCCGTCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGcc                                                                          <  2:12304/100‑1 (MQ=255)
                             gCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCg                                                                       >  2:417248/1‑100 (MQ=255)
                              ccATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGt                                                                      >  1:1284074/1‑100 (MQ=255)
                                    tGCAGCACGCCGACCAGCAACCCGCCAACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTACCCCGTCCAt                                                                  >  1:1454431/1‑98 (MQ=255)
                                       agcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGc                                                             <  1:2383297/100‑1 (MQ=255)
                                       agcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGc                                                             <  2:1445085/100‑1 (MQ=255)
                                        gcaCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGcc                                                            <  1:2219238/100‑1 (MQ=255)
                                         caCTCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGATTCGTGGCGATGGCATCGACGCGAACAAACTTACTTGCCCCGTCCATTGAGCCg                                                           <  2:843317/100‑1 (MQ=255)
                                         caCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCg                                                           <  2:746729/100‑1 (MQ=255)
                                          aCGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGt                                                          >  1:1022042/1‑100 (MQ=255)
                                            gCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTaa                                                        >  1:688936/1‑100 (MQ=255)
                                                 ccAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGTGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATc                                                   >  1:2384320/1‑100 (MQ=255)
                                                  cAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCg                                                  <  1:1714509/100‑1 (MQ=255)
                                                  cAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCg                                                  <  2:1853655/100‑1 (MQ=255)
                                                   aGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCgg                                                 <  2:569275/100‑1 (MQ=255)
                                                    gCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGc                                                >  2:500576/1‑100 (MQ=255)
                                                        cccGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACCAAATTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTcc                                            >  2:1688938/1‑100 (MQ=255)
                                                         ccgccgACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCt                                           <  2:520460/100‑1 (MQ=255)
                                                          cgccgACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTg                                          <  2:2361753/100‑1 (MQ=255)
                                                            ccgACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAg                                        <  1:1174086/100‑1 (MQ=255)
                                                                    gTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACATCg                                >  1:322416/1‑100 (MQ=255)
                                                                      tAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCgg                               >  2:2315095/1‑99 (MQ=255)
                                                                            ccATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAgcggcg                        >  1:2050285/1‑100 (MQ=255)
                                                                                gatgaGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGtttt                    >  1:169375/1‑100 (MQ=255)
                                                                                 atgaGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGttttt                   >  1:1681045/1‑100 (MQ=255)
                                                                                     ggATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGc               >  1:1624026/1‑100 (MQ=255)
                                                                                      gATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGcc              <  1:2390245/100‑1 (MQ=255)
                                                                                       ttcgcGGCGATGCCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCt             >  1:231346/5‑100 (MQ=255)
                                                                                       aTCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCt             <  1:924097/100‑1 (MQ=255)
                                                                                         cgcGGCGATGCCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCa           >  2:1788826/3‑100 (MQ=255)
                                                                                         cgcGGCGATGCCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCa           >  1:1131268/3‑100 (MQ=255)
                                                                                         cgcGGCGATGCCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCa           >  2:1788662/3‑100 (MQ=255)
                                                                                         cgcGGCGATGCCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCa           >  1:1923195/3‑100 (MQ=255)
                                                                                          gcGGCGATGCCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCAt          >  1:1265759/2‑100 (MQ=255)
                                                                                          ccGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCan          <  1:772670/100‑2 (MQ=255)
                                                                                          ccGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCAt          >  1:605619/1‑100 (MQ=255)
                                                                                                aTGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCa     <  2:1934526/99‑1 (MQ=255)
                                                                                                 tGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCAcc   <  1:777210/100‑1 (MQ=255)
                                                                                                  ggCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACca  >  2:1397238/1‑100 (MQ=255)
                                                                                                  ggCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACca  <  2:987954/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TAACTTTCCGCCGCGTGTCCCATGCTCATGCCATGTTGCAGCACGCCGACCAGCAACCCGCCGACAATGTTAATGACCATGATGAGGATCCCGGCGATGGCATCGCCGCGAACAAACTTACTTGCCCCGTCCATTGAGCCGTAAAAATCGGCTTCCTGAGTCACTTCGGAGCGGCGTTTTTTCGCCTCATCTTCACCA  >  NC_000913/1964326‑1964523

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: