Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,130,803 C→T A21T (GCG→ACG)  flgA ← flagellar basal body P‑ring formation protein FlgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,130,8030CT95.5% 57.5 / NA 22A21T (GCG→ACG) flgAflagellar basal body P‑ring formation protein FlgA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/7);  minor base G (0/1);  total (14/8)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

GCATCGCGGCAACCACATAATTTCCTGTGGCCTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGCCGTACTTAACGGACT‑GAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATAT  >  NC_000913/1130557‑1131019
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                         
gCATCGCGGCAACCACATAATTTCCTGTGGCCTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGTCGTa                                                                                                                                                                                                                       <  1:70334/251‑1 (MQ=255)
                  aaTTTCCTGTGGCCTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTTCTCATCGAAAACAATGGCTGCTCGCATAGCGGGAGGAGATTGGGGGCTGTACGAATAACAACAAGAACCTCATCCCCTACCCCCGCGAGCTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGGTTGCTCGTCGGTCTTATCGGACT‑GAACAgg                                                                                                                                                                                                    >  1:48516/1‑250 (MQ=255)
                              ccTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTTAGATTGCTCGTCGTACTTAACGGACT‑GAACAGTATCGCGATGATc                                                                                                                                                                                        <  1:46375/251‑1 (MQ=255)
                                 gTACATTAACCTGTAAATATCGTTTGGCGTTACCGCAGCGTGCCAACACATTCACATTGCACCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAAACTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAAGTGCGATGTGAGATTGCTCGTCGTACTTAACGGACT‑GAACAGTATCGCGATGATCGcc                                                                                                                                                                                     <  1:192409/251‑1 (MQ=255)
                                        aaCCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGACTTATTGCTCATCGAAAGCAATAGCTGCTCGCATGGGGGGAGTAGATTTGGCGCTGTAAGAATAAAAACACGGAACTCATCAATTACCCCCGCGAGTTGGGTGCTAAAAAAATTTGTG‑AAATGGGGTGTGGGATTGTTCGTCGTACTTTATGGGCTAGAAC‑GTATCGCGGTGATCGACACGCtcc                                                                                                                                                                              >  2:28212/1‑249 (MQ=255)
                                               aaaTATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGATGCTCGCATGGCGGTAGAAGATTGGGGGCTGTACGAAAAGAAAAACGAACCTCATCACTTAACCCCGCGGGTTTGGCGCTTAAAAA‑ATTGTGCAATTGCGCTGTGAGAATGCTCGTCGTACATAAAGGACT‑GAACAGGATCGCGATGATCGCCACGCCACGGTttat                                                                                                                                                                       >  1:63080/1‑251 (MQ=255)
                                                   atcttttGTCGTTACCGCAGCGTGCCAACACATTCACATGGCCCCACAGGCGGGAATTATTGCGCAGCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATGGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACT‑GAACAGTAGCGCGATGATCGCCACGCTACGTTTtattat                                                                                                                                                                    <  2:118742/246‑1 (MQ=255)
                                                          gTCGTTACCGCAGCGTGCCAACACATTCACATTGCCTCACAGGCGGGAATTAGTGCTCATCGAAAGCAATGGCTGCTCGCATGGCTGTAGTAGATTGGGCGCTGTTCGAATAGAAACACGAACCTCATCACTTTCCCCCGCGAGTTGGGCTCGAAAAAA‑GTGGTGCAATTGCGAGGTGAGATTGCTCGGCGTGCGTAACGGACT‑TAACAGTAGCGCTATGATCGCCACGGTACGTTTTATTATCAGCAtt                                                                                                                                                             <  2:260483/250‑1 (MQ=255)
                                                          gTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACCAACCTCATCACTTAACCCCCCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGGGAGATTGCTCGTCGTACTTAAAGGACT‑GAACAGTATCGCGATGATCGCCCCGCTCCGTTTTTTTATCAGCAttt                                                                                                                                                            >  1:32743/1‑251 (MQ=255)
                                                                                                                       gAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGAACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACT‑GAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCCTTTTCGCCCCCAGCCATTTCTACAACGTGACTTGTACCTGTCCGCAATGACCATCAACGGCAt                                                                                                >  2:194508/1‑250 (MQ=255)
                                                                                                                         aaGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACT‑GAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATaa                                                                                              >  1:79982/1‑250 (MQ=255)
                                                                                                                                          aTGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACT‑GAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCg                                                                            <  2:32743/251‑1 (MQ=255)
                                                                                                                                                  agtagATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACT‑GAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTcc                                                                    >  1:72623/1‑251 (MQ=255)
                                                                                                                                                      gattgagagcTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGAGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACT‑GAAAAGTATCGCGATGATCGCCACGTTACGTTTTATTATCAGAATTTTCTCCCCCAGCCATTTCTACAACGTGAATTGCACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTccgc                                                                  >  2:76004/9‑249 (MQ=255)
                                                                                                                                                              gcTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAACTGCGCTGCGAGATCGCTCGTCGTACTTAACGGACT‑GAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGACCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCACCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAAcg                                                          >  1:255299/1‑249 (MQ=255)
                                                                                                                                                                   aCGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGCCT‑GAACAGTATCGCGATGATCGCCACACTAAGTTTTATTATCAGCATTTTCCCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAACTAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGt                                                   >  1:194832/1‑251 (MQ=255)
                                                                                                                                                                    cGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACT‑GAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGt                                                   >  1:37639/1‑250 (MQ=255)
                                                                                                                                                                         aGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACT‑GAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAgg                                              <  1:265127/250‑1 (MQ=255)
                                                                                                                                                                                            cTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGCGATTGCTCGTCGTACTTAACGGACT‑GAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGAACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGa                           >  1:204492/1‑250 (MQ=255)
                                                                                                                                                                                                 ccccGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACT‑GAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAAGTTTg                     >  2:209999/1‑251 (MQ=255)
                                                                                                                                                                                                               cgcTAAAAAA‑TGTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTACCGGACT‑GAACAGTGTCGCGATGATCGCCACGCTCCTTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGGGAATTGTTCCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCggagga       <  2:63080/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaaaa‑GTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTACCGGACT‑GAACAGTCTCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGCTAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGatat  >  2:132345/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                         
GCATCGCGGCAACCACATAATTTCCTGTGGCCTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAA‑GTTGTGCAATTGCGATGTGAGATTGCTCGCCGTACTTAACGGACT‑GAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATAT  >  NC_000913/1130557‑1131019

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