Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,134,618 (G)5→6 intergenic (+61/‑111) flgF → / → flgG flagellar basal‑body rod protein FlgF/flagellar basal‑body rod protein FlgG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,134,6131.G93.8% 48.0 / ‑2.8 16intergenic (+56/‑116)flgF/flgGflagellar basal‑body rod protein FlgF/flagellar basal‑body rod protein FlgG
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (9/6);  total (10/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 = 8.75e-01

CAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAG  >  NC_000913/1134533‑1134686
                                                                                 |                                                                         
caaccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTAt                                                                            <  1:226950/81‑1 (MQ=255)
  accaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATgg                                                                         <  1:226951/81‑1 (MQ=255)
             tCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGGCAATc                                                                 <  1:226952/79‑1 (MQ=255)
                  gAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTa                                                         >  1:226978/1‑81 (MQ=255)
                       aaTTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGGCAATCGCAATTTATCGa                                                     >  1:226972/1‑81 (MQ=255)
                            aaaGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGGCAATCGCAATTTATCGAtttt                                                 <  2:226953/80‑1 (MQ=255)
                            aaaGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGGCAATCGCAATTTATCGATTTTg                                                >  2:226986/1‑81 (MQ=255)
                            aaaGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGGCAATCGCAATTTATCGATTTTg                                                >  2:226983/1‑81 (MQ=255)
                                     aTGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGGCAATCGCAATTTATCGATTTTGCGAGCACtt                                       >  2:226981/1‑81 (MQ=255)
                                        aCAAGTATAAGTTGCTCGATGCGCAAGTTTATCGGGTCTATGGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTg                                      >  2:226984/1‑79 (MQ=255)
                                         cAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAg                                   >  1:226979/1‑81 (MQ=255)
                                               taAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCgg                              <  1:226954/80‑1 (MQ=255)
                                                  gTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGAt                            >  1:226994/1‑79 (MQ=255)
                                                          aTGCGCAAGTTTATCGGGTCTATGGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTa                  <  2:226955/81‑1 (MQ=255)
                                                              gcAAGTTTATCGGGTCTATGGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTAcgcc              <  2:226956/81‑1 (MQ=255)
                                                               cAAGTTTATCGGGTCTATGGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTAcgccg             <  1:226957/81‑1 (MQ=255)
                                                                      aTCGGGTCTATGGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGCAGGCCGGATAAGGCGTTTACGCCGCATCCgg      >  2:226991/1‑81 (MQ=25)
                                                                           gTCTATGGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaag  >  1:226985/1‑80 (MQ=255)
                                                                                 |                                                                         
CAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAG  >  NC_000913/1134533‑1134686

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