Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb872,0630TC73.3% 333.8 / 120.6 240R225R (CGT→CGCyliCpredicted peptide transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (0/64);  new base C (0/176);  total (0/240)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 240 total aligned reads displayed.

TGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGG  >  W3110S.gb/872026‑872107
                                     |                                            
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:532390/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:718736/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:705667/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:1341696/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:653401/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:629777/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:1399884/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:1415846/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:625183/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:1389529/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:1522997/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:515197/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:1610517/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:179841/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:43971/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:233076/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:315433/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:280889/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:93765/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:1103086/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:894050/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:1079945/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:1061420/62‑1 (MQ=255)
tGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:1242826/62‑1 (MQ=255)
tGAGAGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:1469417/62‑1 (MQ=255)
 gAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:549016/61‑1 (MQ=255)
 gAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:450478/61‑1 (MQ=255)
 gAGCGAAACCTGGGGTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGAcc                      <  1:377399/61‑1 (MQ=255)
                    cTCAAACACGGTCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1474903/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:48571/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:966446/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:928850/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:499731/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:372543/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:955201/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:477868/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:976729/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:425343/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:406371/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:383834/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:833561/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:754747/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:731536/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:768385/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:777536/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:703237/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:791650/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:673314/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:798895/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:543084/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:20989/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:858239/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:871917/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:884651/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:62250/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:593674/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:586302/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1182955/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:272061/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1351557/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1327772/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1298748/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1296564/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1287458/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1232206/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1226243/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1214347/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:119153/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1356929/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:114607/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1139074/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1129504/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1096010/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1084606/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1073940/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1046027/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1034475/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1024651/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1511338/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1004260/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:204361/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:186141/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:183529/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:170498/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1619542/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1369717/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1574485/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:149501/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:235795/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:254873/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:148332/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1457423/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:14508/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1438424/62‑1 (MQ=255)
                    cTCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1432897/62‑1 (MQ=255)
                     tCAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:158527/61‑1 (MQ=255)
                      cAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:763769/60‑1 (MQ=255)
                      cAAACACGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:352580/60‑1 (MQ=255)
                           acGGGCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACACTTTgg  <  1:227114/55‑1 (MQ=255)
                              ggCTACGCAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTgg  <  1:1254244/52‑1 (MQ=255)
                                     |                                            
TGAGCGAAACCTGGGTTGTCCTCAAACACGGGCTACGTAACGCGATGATCCCGGTAGTGACCATGATGGGCTTACAGTTTGG  >  W3110S.gb/872026‑872107

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