Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,133,116 (A)6→7 coding (543/1209 nt) flgE → flagellar hook protein FlgE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,133,1101.A94.1% 54.8 / ‑2.9 17coding (537/1209 nt)flgEflagellar hook protein FlgE
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (8/8);  total (9/8)
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

GATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTG  >  NC_000913/1133032‑1133176
                                                                               |                                                                   
gATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACa                                                                    >  2:205749/1‑80 (MQ=255)
                aaTTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAAGGTTCGGTg                                                     <  2:205714/79‑1 (MQ=255)
                 aTTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAAGGTTCGGTGACt                                                  >  2:205733/1‑81 (MQ=255)
                  ttCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAAGGTTCGGTGACTg                                                 <  1:205715/81‑1 (MQ=255)
                     cAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAAGGTTCGGTGACTGttt                                              <  1:205716/81‑1 (MQ=255)
                          aTCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAAGGTTCGGTGACTGTTTTCGAc                                         <  2:205717/81‑1 (MQ=255)
                             cGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAAGGTTCGGTGACTGTTTTCGACAGt                                      <  2:205719/81‑1 (MQ=255)
                             cGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAAGGTTCGGTGACTGTTTTCGACAGt                                      <  2:205718/81‑1 (MQ=255)
                              gCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAAGGTTCGGTGACTGTTTTCGACAGt                                      >  2:205735/1‑80 (MQ=255)
                              gCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAAGGTTCGGTGACTGTTTTCGACAGTc                                     <  1:205720/81‑1 (MQ=255)
                                ttCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAg                                   >  1:205745/1‑81 (MQ=255)
                                       tGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAAGGCTCGGTGACTGTTTTCGACAGTCAGGGTAATg                            >  2:205751/1‑81 (MQ=255)
                                              ccATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAGAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTcatga                     >  1:205742/1‑81 (MQ=255)
                                                     gcgcCAGCAATGCGGATAGCTATAACAAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGc              >  1:205741/1‑81 (MQ=255)
                                                      cgcCAGCAATGCGGATAGCTATAACAAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCg             <  1:205721/81‑1 (MQ=255)
                                                                gCGGATAGCTATAACAAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTtgt   >  1:205744/1‑81 (MQ=255)
                                                                   gATAGCTATAACAAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTtgtg  >  1:205752/1‑79 (MQ=255)
                                                                   gATAGCTATAACAAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTtgtg  >  1:205736/1‑79 (MQ=255)
                                                                               |                                                                   
GATGCAGATCAACCTGAATTCCAGTGATCCGCTTCCTACTGTTACGCCATTCAGCGCCAGCAATGCGGATAGCTATAACAAAAAAGGTTCGGTGACTGTTTTCGACAGTCAGGGTAATGCTCATGACATGAGCGTCTACTTTGTG  >  NC_000913/1133032‑1133176

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