Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,007,7870CT72.4% 155.1 / 70.4 105R116R (CGC→CGTfliSflagellar protein potentiates polymerization
Reads supporting (aligned to +/- strand):  ref base C (28/1);  new base T (0/76);  total (28/77)
Fisher's exact test for biased strand distribution p-value = 1.19e-24
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 105 total aligned reads displayed.

AAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTCCCC  >  W3110S.gb/2007732‑2007825
                                                       |                                      
aaGCCTATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1220739/61‑1 (MQ=255)
 aggcAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1465719/57‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCATTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:396506/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAATCATTAATGCGTAatat                                   <  1:589390/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGTATTAATGCGTAatat                                   <  1:2372370/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCGTTAATGCGTAatat                                   <  1:1977321/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2083046/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2694184/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2646300/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2600178/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:764651/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2406159/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2305763/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2108612/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:876548/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2799769/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2055599/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2033862/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2020790/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1993794/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1938331/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:193093/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:326171/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:743185/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:605086/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:561554/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:531072/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:49757/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:491183/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:385096/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:372021/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2744796/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:296910/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:287121/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2843527/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2838777/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2808161/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:974375/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:760277/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1878065/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1469983/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1451870/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1502705/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1511230/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:967626/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1411277/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:157018/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1309606/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1584927/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1602861/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1621196/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1257631/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1237587/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1641064/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1214832/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1212657/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1190364/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1127286/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1118608/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:111394/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1711972/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:962880/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1726312/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:947062/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1785147/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1861874/60‑1 (MQ=255)
 aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGCAatat                                   <  1:1875042/60‑1 (MQ=255)
  gCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2360074/59‑1 (MQ=255)
  gCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:753272/59‑1 (MQ=255)
     aaTTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2298140/56‑1 (MQ=255)
                cGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:757859/45‑1 (MQ=255)
                cGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2756414/45‑1 (MQ=255)
                cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCTCAATATTGCCGATGCCTGGAAag                  >  1:1627468/1‑62 (MQ=255)
                cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGaa                    >  1:833124/1‑60 (MQ=255)
                cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:1746627/1‑62 (MQ=255)
                cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:2825140/1‑62 (MQ=255)
                cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:1002871/1‑62 (MQ=255)
                cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:511621/1‑62 (MQ=255)
                cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:1662989/1‑62 (MQ=255)
                cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:1649668/1‑62 (MQ=255)
                cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:2696340/1‑62 (MQ=255)
                cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:1340671/1‑62 (MQ=255)
                    gTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2249473/41‑1 (MQ=38)
                      ctcCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCg             >  1:2247651/1‑61 (MQ=255)
                      ctcCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGt            >  1:1693150/1‑62 (MQ=255)
                      ctcCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTCGAAAGAGTCGt            >  1:2412116/1‑62 (MQ=255)
                       tcCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1105115/38‑1 (MQ=38)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTActct      >  1:2159450/1‑58 (MQ=255)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACtctc     >  1:872763/1‑59 (MQ=255)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTcccc  >  1:2444306/1‑62 (MQ=255)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTcccc  >  1:1542222/1‑62 (MQ=255)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTccc   >  1:799636/1‑61 (MQ=255)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTccc   >  1:212659/1‑61 (MQ=255)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTccc   >  1:2107793/1‑61 (MQ=255)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTccc   >  1:886946/1‑61 (MQ=255)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTccc   >  1:917565/1‑61 (MQ=255)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTccc   >  1:1890145/1‑61 (MQ=255)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTccc   >  1:1714548/1‑61 (MQ=255)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTccc   >  1:1004/1‑61 (MQ=255)
                                gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGCAAAGAGTCGTTAc         >  1:2775415/1‑55 (MQ=255)
                                                       |                                      
AAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTCCCC  >  W3110S.gb/2007732‑2007825

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