Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,396,8680TG73.6% 149.3 / 50.0 87T261T (ACA→ACCrngribonuclease G
Reads supporting (aligned to +/- strand):  ref base T (0/23);  new base G (0/64);  total (0/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTC  >  W3110S.gb/3396814‑3396897
                                                      |                             
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTTGCGTCCGGTGTa                          <  1:978655/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:414175/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:595912/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:594694/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:828471/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:547920/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:530845/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:830506/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:474809/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:41674/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:826263/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:402809/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:399632/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:379319/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:353931/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:343401/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:84023/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:603365/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1029758/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:620573/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:645743/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:650048/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:657417/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:657528/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:664028/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:664575/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:688055/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:697528/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:708134/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:825908/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:806181/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:749419/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:798041/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1376373/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1132263/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1193386/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1195082/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1198430/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1208440/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1227942/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1232473/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1089660/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1256097/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1068013/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:953379/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1278927/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:129999/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:900967/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:303882/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:874695/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:183170/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:16586/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:164681/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:889741/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1457983/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1387472/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1432159/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1407554/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:949256/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCAGGTGTa                          <  1:619982/60‑1 (MQ=255)
cgctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTCGCGTCCGGTGTa                          <  1:836363/60‑1 (MQ=255)
 gctcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:791667/59‑1 (MQ=255)
   tcgctGAATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:1365456/57‑1 (MQ=255)
                   tCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCGGTGTa                          <  1:309504/41‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:973007/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:989985/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:957518/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:914848/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:998340/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:1432289/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:1036632/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:1081894/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:1172900/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:1251359/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:1263484/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:1265209/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:1395544/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:1409415/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:745245/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:146171/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:233310/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:328861/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:402382/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:513432/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:56817/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:611186/61‑1 (MQ=255)
                       cATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCAtctc  <  1:729265/61‑1 (MQ=255)
                                                      |                             
CGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTC  >  W3110S.gb/3396814‑3396897

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