Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,138,533 C→T G52G (GGC→GGT flgK → flagellar hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,138,5330CT97.0% 91.9 / ‑4.3 34G52G (GGC→GGTflgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/13);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

TGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGC  >  NC_000913/1138290‑1138767
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
tGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCggttgggttg                                                                                                                                                                                                                                            >  1:303568/1‑251 (MQ=255)
                     gAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTTTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGTATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCt                                                                                                                                                                                                                       <  1:176691/250‑1 (MQ=255)
                        tAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCAATGGTGTCTACGTTTTTGgtg                                                                                                                                                                                                                   >  1:326269/1‑251 (MQ=255)
                                               cccGATAATACTTTGTAATTGAAGGCTTATAAGGAACCTCCATGGCCATCTTTATTAATAACGACATGAGCGGACTGAGCGCGGCCCAGGCGGCGTAAAATACGGAAAGAAATATTATCTTCAGCTAAAACTTTGCCGGAGATCCCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGttt                                                                                                                                                                                          <  2:271803/250‑1 (MQ=255)
                                                   ataataCTCTGTAATTGAAGTTTTAAGAGGACCCTCCATGTCAAGCTTGATTAATAACGCCATGAGCGGACTGACCGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATGATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTttatta                                                                                                                                                                                      <  2:161799/251‑1 (MQ=255)
                                                    taataCTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGCTTGGCAATGGGGTCTACTTTTCTGGTGGGCAGCGGGAGTATGAGTGGTttatta                                                                                                                                                                                      >  2:111277/1‑250 (MQ=255)
                                                       taCTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGATTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGTGGTGGGCATGGTGTATACGTTTCTGGTGTGCGGCGAGGGGATTATGCGTttattaaca                                                                                                                                                                                   >  1:199378/1‑247 (MQ=255)
                                                             atAATGGAAGGATTATAAGGACCCCCCATGTCCAGCGTGATTAATAACGCCATGAGCGGAATGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTGGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGGGCAGCGTGAGTATGATGCGTTTATTACCAACCAGtt                                                                                                                                                                            <  2:326269/250‑1 (MQ=255)
                                                                              aGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCAATTGGGTCCACGTTTCTGTTGTGTACAGTGGGTGTTATGCGTTTATTACCAACCCGTCTCTTGCGGCGCAGACGc                                                                                                                                                          >  2:215545/1‑251 (MQ=255)
                                                                                   ccTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGGGTCCACGTTTCTGGTGTGCCGCGGGAGTATGGTGGGTTTATTACCAACCAGTTACGTgcggcggcggcgcca                                                                                                                                                         >  2:499534/1‑239 (MQ=255)
                                                                                        aggTCCAGCGTGATTAATAACTCCTTGATCGGATTGAACGCGCCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCAATAACGTTGCAGTATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCgg                                                                                                                                                <  1:272966/249‑1 (MQ=255)
                                                                                                       aataaCGCCATGAGCGGACTGAACGCGGCCCAGCCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTGGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCt                                                                                                                                  <  1:355095/250‑1 (MQ=255)
                                                                                                              ccATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATACTATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGATGGCAATTGGGTCTACGGTTCTGGTGTGCATCGGACGTATGATGAGTTTTTTACCTACCAGTTACCTGCGGCGCCAACACAACGTAGCGGTATGACCGCCCGCttagcgc                                                                                                                            >  1:61415/1‑243 (MQ=255)
                                                                                                               cATGAGAGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTGGGTTGGCAATGGAGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCCGACGCAAAGTTGCGGTCTTACTGCCCGcatagcgcca                                                                                                                          >  2:333849/1‑241 (MQ=255)
                                                                                                                                    ccAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAAACACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGCGTCATAGTTTCTGGTAGGCAGAGTGAGTATGGTGAGTTTATTACCAACCAGTTACGTGCGGCGCAGACGACAAATAACGGTCTGACTGCACGCGCTGAGAAAGTGTCGAAAATCGACAAAATGc                                                                                                    >  1:334221/1‑251 (MQ=255)
                                                                                                                                      aGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCa                                                                                                                                                                                                               <  1:437286/145‑1 (MQ=255)
                                                                                                                                            cgTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGGTATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGACGCTGACGGGTGGGGTGGCTAGTGTGTCTCCGTTTGTGGGGTGCGGCGTGGGTAAGATATATTTATTGCAAAACCGTTAAGTGCGGCGCCGGCGCCAAGTAGCGGGCTGTCTGCCCGCTTAGAGAAGATGGCGAACATCgtagaaaaggcagtcacc                                                                                            >  2:104467/1‑234 (MQ=255)
                                                                                                                                            cgTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCAGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGGGTCTACGGTTCTGGTGTGCAGCGTGAGTATGGTGCGTTTATTACCAACCCGTTACGTGGGGGGCACACGCAAAATAGCGATCTGGCTGCCCGCCATAAGCAGATGTCGAAAAtaggaaataagctccccccc                                                                                            >  2:305259/1‑232 (MQ=255)
                                                                                                                                                     cGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCAAAGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCCGTTACGTGCGGCGCCGACGCAAAGTAGCGGGCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACctct                                                                                   >  1:486025/1‑251 (MQ=255)
                                                                                                                                                             aataatATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGGGTATGATGCGTTTATTACCAACCCGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTTACTGCCCGCTCTGAGCAGATGTCGAAAATTGAACATATGCTCTCCACCCATACCTCTTCGCTGGc                                                                           >  1:395120/1‑251 (MQ=255)
                                                                                                                                                                       ccAGCTATAACGTTGCCGGATATCCCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATg                                                                 <  1:147666/251‑1 (MQ=255)
                                                                                                                                                                         aGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGGGAGgatgat                                                                                                                                                                                                >  1:457886/1‑123 (MQ=255)
                                                                                                                                                                                aCGTTGCCGAATATACCCGCCAAACCACTATGATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATttc                                                        <  1:104547/251‑1 (MQ=255)
                                                                                                                                                                                      ccGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGAAAAAGTGTCTGCGTTTCTGGGGTGCAGCGGGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGAATGCCCGCTATGAGCAGATGTCGTAAATCGACAATATGCTCTCCACCAGAACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTTAcc                                                 >  1:325665/1‑251 (MQ=255)
                                                                                                                                                                                          atatACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGGGCAGCGGGAGGATGATGCGTTTATTACCACACAGTTACGTGCGGCGCAGACGCACAGTAGAGGTCTGACTGCCCGCCATGAGCAGATGTCGAAAATCAACACTACTATCTCCCCCAGTAACTCTTCCCTGGCAACACAGATGCAGGCTTTCTTTACCaccc                                             >  2:109591/1‑248 (MQ=255)
                                                                                                                                                                                               cccGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGt                                                                                                                                                                             <  1:158226/120‑1 (MQ=255)
                                                                                                                                                                                               cccGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCAATGGTGTCTACCGTTCTGGGGTGCCGCGTGAGTATGGTGGGGTTATTACCAACCAGt                                                                                                                                                                             >  2:158226/1‑120 (MQ=255)
                                                                                                                                                                                                 cGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAAc                                     <  1:408216/251‑1 (MQ=255)
                                                                                                                                                                                                  gtcAAACCTCTTTTATGGCGCAGGCCAATAGCACGTGGGGCGCTGGCGGTTGGGGTGGCAATGGTGGCGACGGTTCTGGTGGGCAGCGGGAGTAGGATGCGTTTATTAACATCCAGTTACGTGCGGCGCAGACGCATAGTAGCGGTCTGACTGCACGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTACCCCATTTCCTCTTCGCTGGCAACACGGATGCAGGATTTCTCCACCAGCCTGCaaccg                                    <  2:489243/249‑4 (MQ=255)
                                                                                                                                                                                                       aCCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTggg                               >  1:361648/1‑250 (MQ=255)
                                                                                                                                                                                                                atGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTATCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAAcg                       <  1:7982/250‑1 (MQ=255)
                                                                                                                                                                                                                          ccAATAGCACGTTGGGCGCTGGCGGTTGGGGTGGCACTGGCTTCTACGCTTATGGTGTGAAGCGGGAGTATAATGCGTTTATTACCACCCAGGTACGTTCGGCGCAGGCGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAAATGACAATATGCTCTCCACCAGTCACTCTTCGCTGGCACCACAGATTCAGGATTTTTTCACCAGCCTGCAAAAGATGGAGCGTAACGCGGGAAACCCg           >  1:417519/1‑251 (MQ=255)
                                                                                                                                                                                                                           cAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACCTTTCTGGTGGGAAGCGGGGGTATGATGGGGTTATTAACAACCAGTTACGTGCGGCGCAGAAGCAAAGGAGCGGTCTGACTGCCCCCTACGAGCAGGTGTCGGAAAACGACAATATGCTCTACACCCGTACCTCTTCGCTGGCAACACAAATGCAGGATTTCTTCACCAGCCTGCAAACGCTGATGAGATACGCGGTAGACTCgg          >  1:266904/1‑251 (MQ=255)
                                                                                                                                                                                                                            aaTAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAAAGGTGTCTAAGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTTTTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGGCAATATGCTCTCCACCAGTACCTCTTCGCTGGCACCACAGATGCAGGATTTCTTCACCCGCCTGCAAACGCTGGGGGGGAACACGGAAGACCCGGc         >  2:270109/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cGTTGGGCGCTGGCGGTTGGGCTGGCAATGGTGTCTTCGTTTCTGGTGTGCCGCGTGAGGATGATGCGTTTATTACCAACCAGTTACGCGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCCGATGTCGAAAATCGACAATATGACCTCCACCAGTACCTCTTCGCTGGCCACACAGTTGCAGGTTTTCTTCAACCACCTGCCAAACCTGGTGGGTAACCAGGAAGCACCCGCAgcgccc  >  2:167382/1‑249 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  
TGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGTGAGTAACGCGGAAGACCCGGCAGCGCGC  >  NC_000913/1138290‑1138767

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