Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,007,7690AG74.0% 171.5 / 63.1 104E110E (GAA→GAGfliSflagellar protein potentiates polymerization
Reads supporting (aligned to +/- strand):  ref base A (27/0);  new base G (0/77);  total (27/77)
Fisher's exact test for biased strand distribution p-value = 1.53e-25
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 117 total aligned reads displayed.

GCGCTTGCTGCAAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTCCCC  >  W3110S.gb/2007721‑2007825
                                                |                                                        
gcgcTTGCTGCAAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTaa                                             >  1:2635012/1‑62 (MQ=255)
gcgcTTGCTGCAAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTaa                                             >  1:2375957/1‑62 (MQ=255)
gcgcTTGCTGCAAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTaa                                             >  1:1703872/1‑62 (MQ=255)
gcgcTTGCTGCAAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTaa                                             >  1:512777/1‑62 (MQ=255)
gcgcTTGCTGCAAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTaa                                             >  1:1262983/1‑62 (MQ=255)
gcgcTTGCTGCAAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTaa                                             >  1:1229024/1‑62 (MQ=255)
gcgcTTGCTGCAAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTaa                                             >  1:1931300/1‑62 (MQ=255)
gcgcTTGCTGCAAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTaa                                             >  1:698717/1‑62 (MQ=255)
gcgcTTGCTGCAAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTaa                                             >  1:1623061/1‑62 (MQ=255)
gcgcTTGCTGCAAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTaa                                             >  1:2034422/1‑62 (MQ=255)
           aaGCCTATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1220739/61‑1 (MQ=255)
            aggcAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1465719/57‑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:2600178/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2744796/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2646300/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2694184/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:287121/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2492864/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:947062/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2055599/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2033862/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:49757/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:876548/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:764651/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:760277/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:501759/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2799769/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:326171/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:296910/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:974375/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2843527/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2808161/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1309606/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1641064/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1621196/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1584927/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:157018/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1511230/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1502705/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1451870/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1411277/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1993794/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1237587/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:1878065/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1785147/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:962880/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:193093/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:1726312/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2020790/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:967626/60‑1 (MQ=255)
            aGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGCAatat                                   <  1:1875042/60‑1 (MQ=255)
             gCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:753272/59‑1 (MQ=255)
             gCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2360074/59‑1 (MQ=255)
                aaTTTACGCAACGATGTCTCCGCAGTCGAAGAGGTGGAAGCATTAATGCGTAatat                                   <  1:2298140/56‑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:1002871/1‑62 (MQ=255)
                           cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:1662989/1‑62 (MQ=255)
                           cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:511621/1‑62 (MQ=255)
                           cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:1340671/1‑62 (MQ=255)
                           cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:1514652/1‑62 (MQ=255)
                           cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:1746627/1‑62 (MQ=255)
                           cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:2825140/1‑62 (MQ=255)
                           cGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAag                  >  1:1649668/1‑62 (MQ=255)
                                 ctcCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCg             >  1:2247651/1‑61 (MQ=255)
                                 ctcCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGt            >  1:1935422/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:1542222/1‑62 (MQ=255)
                                           gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTcccc  >  1:2444306/1‑62 (MQ=255)
                                           gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTccc   >  1:917565/1‑61 (MQ=255)
                                           gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTccc   >  1:212659/1‑61 (MQ=255)
                                           gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTccc   >  1:799636/1‑61 (MQ=255)
                                           gaagaaGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTccc   >  1:2107793/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)
                                                |                                                        
GCGCTTGCTGCAAGCCAATTTACGCAACGATGTCTCCGCAGTCGAAGAAGTGGAAGCATTAATGCGCAATATTGCCGATGCCTGGAAAGAGTCGTTACTCTCCCC  >  W3110S.gb/2007721‑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: