Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,866,4010GA67.5% 44.2 / 35.2 40Q337* (CAA→TAA) nlpDpredicted outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/13);  new base A (14/13);  total (14/26)
Fisher's exact test for biased strand distribution p-value = 1.09e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGAT  >  W3110S.gb/2866354‑2866467
                                               |                                                                  
ggTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt                                             <  1:776121/71‑1 (MQ=255)
                    ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATtg                                                   >  1:1097116/1‑45 (MQ=255)
                    ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATtg                                                   >  1:1112268/1‑45 (MQ=255)
                    ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATtg                                                   >  1:3552129/1‑45 (MQ=255)
                    ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATtg                                                   >  1:3734980/1‑45 (MQ=255)
                    ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATtg                                                   >  1:1583806/1‑45 (MQ=255)
                    ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt                                             <  1:850629/51‑1 (MQ=255)
                    ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt                                             <  1:3110721/51‑1 (MQ=255)
                    ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt                                             <  1:1992884/51‑1 (MQ=255)
                    ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt                                             <  1:1729207/51‑1 (MQ=255)
                    ttttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt                                             <  1:1725412/51‑1 (MQ=255)
                      ttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCa                                                      >  1:1388661/1‑40 (MQ=255)
                      ttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt                                             >  1:1478827/1‑49 (MQ=255)
                      ttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt                                             >  1:209454/1‑49 (MQ=255)
                      ttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt                                             >  1:2201844/1‑49 (MQ=255)
                      ttttGCCCCGCCTTAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTa                                           >  1:1390484/1‑51 (MQ=255)
                      ttttGCCCCGCCTAAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt                                             >  1:1632646/1‑49 (MQ=255)
                                   tAACTTCTTGTTATTCCCGGACCAGCATTGTGTCGt                                             <  1:3062960/36‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:4085314/35‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:3574827/35‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:2491828/35‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:3749229/35‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:1752391/35‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:3337507/35‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:419988/35‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:689165/35‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:1956167/35‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:1569855/35‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:141030/35‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:136057/35‑1 (MQ=255)
                                    aaCTTCTTGTTGTTCCCGGACCAGCATTGTGTCGt                                             <  1:106816/35‑1 (MQ=255)
                                    aaCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATgg                                        >  1:2660763/1‑40 (MQ=255)
                                    aaCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt         >  1:3293165/1‑71 (MQ=255)
                                    aaCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGtt         >  1:1282330/1‑71 (MQ=255)
                                     aCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGttt        <  1:1086840/71‑1 (MQ=255)
                                     aCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGttt        <  1:2978763/71‑1 (MQ=255)
                                     aCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGttt        <  1:2826075/71‑1 (MQ=255)
                                     aCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGttt        <  1:268498/71‑1 (MQ=255)
                                     aCTTCTTGTTATTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGttt        <  1:1355826/71‑1 (MQ=255)
                                       ttcttGTTATTCCCGGACCAGCATTGTGTCGTTATg                                         <  1:1526964/36‑1 (MQ=255)
                                            gttattCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat  <  1:3303581/70‑1 (MQ=255)
                                            gttattCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat  <  1:2728014/70‑1 (MQ=255)
                                            gttattCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat  <  1:4108541/70‑1 (MQ=255)
                                            gttattCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat  <  1:2394376/70‑1 (MQ=255)
                                            gttattCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTgatgat  <  1:2225912/70‑1 (MQ=255)
                                               |                                                                  
GGTGCTACCCATGGTCGCTATTTTTTGCCCCGCCTTAACTTCTTGTTGTTCCCGGACCAGCATTGTGTCGTTATGGGCGTAGGCACTCAGGTAATCATCATTATGTTTGATGAT  >  W3110S.gb/2866354‑2866467

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: