Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,741,8790TC30.3% 151.5 / 94.7 109S30P (TCG→CCG) rplV50S ribosomal subunit protein L22
Reads supporting (aligned to +/- strand):  ref base T (22/54);  new base C (33/0);  total (55/54)
Fisher's exact test for biased strand distribution p-value = 2.07e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 109 total aligned reads displayed.

GTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACCAACAA  >  minE/2741823‑2741913
                                                        |                                  
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGCCCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:673421/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:2311128/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:3638589/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:1673868/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:1674753/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:3622431/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:2362480/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:183258/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:3602679/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:2392042/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:3310899/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:3035833/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:196243/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:2257607/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:2009980/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:2210214/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:2135987/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:2774884/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:958024/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:265024/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:110290/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:898338/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:1158870/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:2473843/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:2760597/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:125269/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:1292339/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:1360308/1‑59 (MQ=255)
gTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGCCggt                                >  1:1509279/1‑59 (MQ=255)
               aGGTTCGCCTTGTTGCTTACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2654458/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGGAAGTGTCGCAGGCTCTGGa                       <  1:1071727/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2958959/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:3043896/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:3047525/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2687513/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2872484/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:3490212/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2927283/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:390390/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:98614/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:911911/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:91140/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:798761/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:790856/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:666780/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:576280/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:54709/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:3122869/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:3614186/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:106995/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:3352413/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:328744/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:3245368/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:3202887/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:312965/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1941400/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1056350/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1078233/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:112392/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1238390/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1331193/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1558855/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1590460/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1662342/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1763821/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1774678/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1839197/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2637234/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1951934/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2100410/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2578427/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2404575/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2273488/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2103136/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2131758/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2218988/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2504587/55‑1 (MQ=255)
               aGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:2208430/55‑1 (MQ=255)
               aGGTGCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:604184/55‑1 (MQ=255)
                       cTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGa                       <  1:1072395/47‑1 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTg                >  1:407410/1‑47 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:1354735/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:1242531/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:2788626/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:1158806/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:981983/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:2996577/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:2329887/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:3153669/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:1670958/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:3512396/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:187269/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:1887944/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:3264463/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:1986207/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:3234384/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:3131728/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACcaacaa  >  1:3191010/1‑61 (MQ=255)
                              cTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGAAATTTTGACCTACACcaacaa  >  1:2523636/1‑61 (MQ=255)
                              cTGACCTGAATCGCGGTAAGAAAGTGTCGCAGGCTCTCGATATTTTGACCTACACcaacaa  >  1:2785709/1‑61 (MQ=255)
                                                        |                                  
GTTCTTCTGCTCAGAAGGTTCGCCTTGTTGCTGACCTGATTCGCGGTAAGAAAGTGTCGCAGGCTCTGGATATTTTGACCTACACCAACAA  >  minE/2741823‑2741913

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: