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,5330CT91.7% 24.4 / ‑6.2 12G52G (GGC→GGTflgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (5/6);  total (6/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGT  >  NC_000913/1138286‑1138739
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                               
aGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGACAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTgg                                                                                                                                                                                                              >  1:203229/1‑251 (MQ=255)
aGTAGGAACATTGATAACCTGTTCTGAATAACTCAAGGCGGGCGGGGGGCTGCCAATAATACTCTGTAATTGAAGGCTTATAAGAAACCCCCATGCCCAGCTTGATTAATAACGCCAGGAGCGGACGGAGCGCGGCCCAGGCGGCTTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATGCCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTgg                                                                                                                                                                                                              <  1:110374/251‑1 (MQ=255)
                                         gcggGTCGCTGCCGATAATCCTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTa                                                                                                                                                                     <  1:254893/251‑1 (MQ=255)
                                                  tGCCGCTAATACTCTGTAATTGACGGCTTATAAGGCACCTCCATGTCCAGCATGATTAATCACGCAATGAGCGTACTGAACGCGGCTCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATAACCGCCAAACCACTATTATGTCGCAGGCCAATCGCACGTTGGGCTCTGGCGGTTGGGTTGGCACAGGTGTCTACCTTTCTTGTGAGCTGCGGGTGTATGATGCggt                                                                                                                                                            >  2:261437/1‑249 (MQ=255)
                                                                                                                      gAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATg                                                                                         <  2:205825/250‑1 (MQ=255)
                                                                                                                            acggAACGCGGCCCAGGCGGCGTTAAATACGGAAAGTAATAATATCACCAGCTAAAGCGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCGGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTAGGATGCGTTTAGTTCCATCCAGTTACGTGCGGCGCAGACGCAAAGGAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGaaaa                                                                                  <  2:76467/248‑1 (MQ=255)
                                                                                                                                aaCGCGGCCCACGCGGCGTTGAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGGGTCTATGTTTCTGGTGGGCAGCGTGTGGATGATGCGTTTCTTACCAACCCGTTACGGGCGGGGCCGACGCAAAGTAGAGGTCTGGCTGCCCGCTCTGAGCCGATGTCGCAAATCgg                                                                              >  2:156997/1‑250 (MQ=255)
                                                                                                                                     ggCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAata                                                                        <  2:219029/251‑1 (MQ=255)
                                                                                                                                                                            cAGCGATAACTTTGCCGGATATACCCGCCAAACCACGATGATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAAGCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGc                                 <  1:142997/251‑1 (MQ=255)
                                                                                                                                                                                            ggATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGTTTCTGGTGGGCAGCGTGAGGATGATGCGTTTATTACCAACCAGTTAAGGGCGGCGCAGACGCAAAGTAGGGGTCTGAATGCCCGCTATGAGCAGATGTCGCAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGGTTTCTTCAcc                   >  1:105306/1‑249 (MQ=255)
                                                                                                                                                                                                          aaCCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGTGTCTACGGTTCTGGTGTGCAGCGTGCGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCCATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCCCCAGCCTGCAAACGCTgg   >  1:27437/1‑251 (MQ=255)
                                                                                                                                                                                                           aCCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGTTGGGTTGGCAATGGGGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGAAGCAAAGTAGCGGTCTGACTGCCCGCTTTGAGCAGATGTCGAAAATCGAACAAATGCTCTCCACCAGTACCTCCTCGCCGGCCACACCGATGCAGGATTTTTTCACCCACCTGCAAACGCTggg  >  2:75753/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                               
AGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAATAACGCCATGAGCGGACTGAACGCGGCCCAGGCGGCGTTAAATACGGCAAGTAATAATATCTCCAGCTATAACGTTGCCGGATATACCCGCCAAACCACTATTATGGCGCAGGCCAATAGCACGTTGGGCGCTGGCGGCTGGGTTGGCAATGGTGTCTACGTTTCTGGTGTGCAGCGTGAGTATGATGCGTTTATTACCAACCAGTTACGTGCGGCGCAGACGCAAAGTAGCGGTCTGACTGCCCGCTATGAGCAGATGTCGAAAATCGACAATATGCTCTCCACCAGTACCTCTTCGCTGGCAACACAGATGCAGGATTTCTTCACCAGCCTGCAAACGCTGGT  >  NC_000913/1138286‑1138739

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