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

GCATCGCGGCAACCACATAATTTCCTGTGGCCTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGCCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATAT  >  NC_000913/1130557‑1131019
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                         
gCATCGCGGCAACCACATAATTTCCTGTGGCCTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTa                                                                                                                                                                                                                       <  1:70334/251‑1 (MQ=255)
                  aaTTTCCTGTGGCCTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTTCTCATCGAAAACAATGGCTGCTCGCATAGCGGGAGGAGATTGGGGGCTGTACGAATAACAACAAGAACCTCATCCCCTACCCCCGCGAGCTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGGTTGCTCGTCGGTCTTATCGGACTGAACAgg                                                                                                                                                                                                    >  1:48516/1‑250 (MQ=255)
                              ccTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTTAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATc                                                                                                                                                                                        <  1:46375/251‑1 (MQ=255)
                                 gTACATTAACCTGTAAATATCGTTTGGCGTTACCGCAGCGTGCCAACACATTCACATTGCACCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAAACTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAAGTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGcc                                                                                                                                                                                     <  1:192409/251‑1 (MQ=255)
                                        aaCCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGACTTATTGCTCATCGAAAGCAATAGCTGCTCGCATGGGGGGAGTAGATTTGGCGCTGTAAGAATAAAAACACGGAACTCATCAATTACCCCCGCGAGTTGGGTGCTAAAAAAATTTGTGAAATGGGGTGTGGGATTGTTCGTCGTACTTTATGGGCTAGAACGTATCGCGGTGATCGACACGCtcc                                                                                                                                                                              >  2:28212/1‑249 (MQ=255)
                                               aaaTATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGATGCTCGCATGGCGGTAGAAGATTGGGGGCTGTACGAAAAGAAAAACGAACCTCATCACTTAACCCCGCGGGTTTGGCGCTTAAAAAATTGTGCAATTGCGCTGTGAGAATGCTCGTCGTACATAAAGGACTGAACAGGATCGCGATGATCGCCACGCCACGGTttat                                                                                                                                                                       >  1:63080/1‑251 (MQ=255)
                                                   atcttttGTCGTTACCGCAGCGTGCCAACACATTCACATGGCCCCACAGGCGGGAATTATTGCGCAGCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATGGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTAGCGCGATGATCGCCACGCTACGTTTtattat                                                                                                                                                                    <  2:118742/246‑1 (MQ=255)
                                                          gTCGTTACCGCAGCGTGCCAACACATTCACATTGCCTCACAGGCGGGAATTAGTGCTCATCGAAAGCAATGGCTGCTCGCATGGCTGTAGTAGATTGGGCGCTGTTCGAATAGAAACACGAACCTCATCACTTTCCCCCGCGAGTTGGGCTCGAAAAAAGTGGTGCAATTGCGAGGTGAGATTGCTCGGCGTGCGTAACGGACTTAACAGTAGCGCTATGATCGCCACGGTACGTTTTATTATCAGCAtt                                                                                                                                                             <  2:260483/250‑1 (MQ=255)
                                                          gTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACCAACCTCATCACTTAACCCCCCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGGGAGATTGCTCGTCGTACTTAAAGGACTGAACAGTATCGCGATGATCGCCCCGCTCCGTTTTTTTATCAGCAttt                                                                                                                                                            >  1:32743/1‑251 (MQ=255)
                                                                                                                       gAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGAACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCCTTTTCGCCCCCAGCCATTTCTACAACGTGACTTGTACCTGTCCGCAATGACCATCAACGGCAt                                                                                                >  2:194508/1‑250 (MQ=255)
                                                                                                                         aaGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATaa                                                                                              >  1:79982/1‑250 (MQ=255)
                                                                                                                                          aTGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCg                                                                            <  2:32743/251‑1 (MQ=255)
                                                                                                                                                  agtagATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTcc                                                                    >  1:72623/1‑251 (MQ=255)
                                                                                                                                                      gattgagagcTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGAGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAAAAGTATCGCGATGATCGCCACGTTACGTTTTATTATCAGAATTTTCTCCCCCAGCCATTTCTACAACGTGAATTGCACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTccgc                                                                  >  2:76004/9‑249 (MQ=255)
                                                                                                                                                              gcTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAACTGCGCTGCGAGATCGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGACCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCACCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAAcg                                                          >  1:255299/1‑249 (MQ=255)
                                                                                                                                                                   aCGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGCCTGAACAGTATCGCGATGATCGCCACACTAAGTTTTATTATCAGCATTTTCCCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAACTAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGt                                                   >  1:194832/1‑251 (MQ=255)
                                                                                                                                                                    cGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGt                                                   >  1:37639/1‑250 (MQ=255)
                                                                                                                                                                         aGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAgg                                              <  1:265127/250‑1 (MQ=255)
                                                                                                                                                                                            cTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGCGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGAACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGa                           >  1:204492/1‑250 (MQ=255)
                                                                                                                                                                                                 ccccGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAAGTTTg                     >  2:209999/1‑251 (MQ=255)
                                                                                                                                                                                                               cgcTAAAAAATGTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTACCGGACTGAACAGTGTCGCGATGATCGCCACGCTCCTTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGGGAATTGTTCCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCggagga       <  2:63080/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaaaaGTTGTGCAATTGCGATGTGAGATTGCTCGTCGTACTTACCGGACTGAACAGTCTCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGCTAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGatat  >  2:132345/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                         
GCATCGCGGCAACCACATAATTTCCTGTGGCCTGTACATTAACCTGTAAATATCGTTTGTCGTTACCGCAGCGTGCCAACACATTCACATTGCCCCACAGGCGGGAATTATTGCTCATCGAAAGCAATGGCTGCTCGCATGGCGGTAGTAGATTGGGCGCTGTACGAATAGAAACACGAACCTCATCACTTACCCCCGCGAGTTGGGCGCTAAAAAAGTTGTGCAATTGCGATGTGAGATTGCTCGCCGTACTTAACGGACTGAACAGTATCGCGATGATCGCCACGCTACGTTTTATTATCAGCATTTTCGCCCCCAGCCATTTCTACAACGTGAATTGTACCTGTCCGCAATGACCATCAACGGCATAAATAGCGACCCATTTTGCGTTTATTCCGCCGATAACGCGCGCGTAAAGGCATTTAAGCTGATGGCAGAATTTTGATACCTGCGGAGGAGATAT  >  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: