Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,996,2240GA73.0% 274.1 / 117.6 226Q330Q (CAG→CAAprlColigopeptidase A
Reads supporting (aligned to +/- strand):  ref base G (11/50);  new base A (165/0);  total (176/50)
Fisher's exact test for biased strand distribution p-value = 8.84e-40
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-03
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 226 total aligned reads displayed.

CGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGC  >  W3110S.gb/3996183‑3996275
                                         |                                                   
cGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTtac                                   >  1:278279/1‑60 (MQ=255)
cGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTtac                                   >  1:841556/1‑60 (MQ=255)
cGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTtac                                   >  1:83137/1‑60 (MQ=255)
cGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTtac                                   >  1:314025/1‑60 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCg                            >  1:1043268/3‑60 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:324664/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:519729/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:505061/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:493213/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:469876/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:456595/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:408001/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:406470/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:394325/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:381965/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:35265/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:333495/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:576911/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:30033/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:271866/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1011149/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:249485/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:221193/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:210829/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:730559/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:982544/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:966467/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:957411/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:9503/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:918552/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:8947/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:860503/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:782650/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:758604/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:738100/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:545551/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:667057/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:656192/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:630820/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:624961/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:610437/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:599281/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:591291/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:26701/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:553210/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:120986/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1424127/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1383120/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1370522/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:206566/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1342186/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1326336/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1306170/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1290178/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1272092/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1257264/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:123473/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1353851/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1201803/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1168473/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1159435/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1139717/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1127797/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1108410/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1078373/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1073635/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1035974/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1025667/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:201466/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1629800/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1651326/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1612597/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1584999/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:171319/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1558436/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1546110/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1535020/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1521102/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1507836/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1478455/3‑61 (MQ=255)
       tcgCCAAAGCCGAATTTGGCGTCGATGAGTTACAACCGTGGGATATCGCTTACTACAGCGa                           >  1:1464293/3‑61 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:940267/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:227563/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:93008/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:1063862/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:905239/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:19837/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:1640847/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:1604579/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:273663/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:485172/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:1232464/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:378540/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:1335346/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:1377152/61‑1 (MQ=255)
                                tGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:1436641/61‑1 (MQ=255)
                                 gtgTTGCAGCCGTGGGATATCGCTTACTAGAGCGAAATACATAAACAGCACCTCTACAGc  <  1:439854/58‑1 (MQ=255)
                                 gAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:805794/60‑1 (MQ=255)
                                 gAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:795929/60‑1 (MQ=255)
                                 gAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGc  <  1:714709/60‑1 (MQ=255)
                                         |                                                   
CGTGCCTTTGCCAAAGCCGAATTTGGCGTCGATGAGTTGCAGCCGTGGGATATCGCTTACTACAGCGAAAAACAAAAACAGCACCTCTACAGC  >  W3110S.gb/3996183‑3996275

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: