Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,413,3700AT73.0% 61.8 / 25.2 37N679Y (AAT→TAT) rnrexoribonuclease R, RNase R
Reads supporting (aligned to +/- strand):  ref base A (10/0);  new base T (0/27);  total (10/27)
Fisher's exact test for biased strand distribution p-value = 2.87e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCA  >  W3110S.gb/4413308‑4413428
                                                              |                                                          
gcttCTTCGTCCGTCTGGACGTCTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1929809/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1745916/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:273406/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1413344/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1452678/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1497554/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:149860/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:806186/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1555119/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1594178/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:2441442/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1491291/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1750044/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:2415082/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1777877/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1826760/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:698909/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1850654/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:2093660/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:2079962/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1181819/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1054302/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1096647/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:5821/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1112090/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1383525/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:497890/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:464571/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:411005/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1340869/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:134688/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:407181/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:610953/64‑1 (MQ=255)
gcttCTTCGTCCGTCTGGACGACTTGTTCATTGATCGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1559777/64‑1 (MQ=255)
gcttCTTCGTCCGGCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:105745/64‑1 (MQ=255)
gcttCTTCGTCCGGCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1074402/64‑1 (MQ=255)
 cttcttCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:2489959/63‑1 (MQ=255)
 cttcttCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1836340/63‑1 (MQ=255)
                ggACGACTTGTTCTTTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1056276/48‑1 (MQ=255)
                   cgacTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1877145/45‑1 (MQ=255)
                         gTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACaa                                                           <  1:1171478/39‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAATGCCTGATGGGGGa              <  1:42337/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:287412/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:2112621/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:672385/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:2142828/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:2235061/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:591523/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:2467459/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1025523/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1228645/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1234156/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1188834/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1352070/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1366145/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1111165/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1402463/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:204131/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1638602/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:788328/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1758204/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1805216/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1817054/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:702337/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1831980/65‑1 (MQ=255)
                                            cATGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:191693/65‑1 (MQ=255)
                                             aTGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:636034/64‑1 (MQ=255)
                                              tGTCTCTTCGCTGGACTATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGa              <  1:1319353/63‑1 (MQ=255)
                                                        cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATcc          >  1:1382007/1‑57 (MQ=255)
                                                        cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCa  >  1:983848/1‑65 (MQ=255)
                                                        cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCa  >  1:307032/1‑65 (MQ=255)
                                                        cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCa  >  1:2538284/1‑65 (MQ=255)
                                                        cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCa  >  1:2118392/1‑65 (MQ=255)
                                                        cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCa  >  1:183456/1‑65 (MQ=255)
                                                        cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCa  >  1:1822863/1‑65 (MQ=255)
                                                        cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCa  >  1:1639480/1‑65 (MQ=255)
                                                        cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCa  >  1:1503819/1‑65 (MQ=255)
                                                        cTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCa  >  1:1481397/1‑65 (MQ=255)
                                                              |                                                          
GCTTCTTCGTCCGTCTGGACGACTTGTTCATTGATGGTCTGGTCCATGTCTCTTCGCTGGACAATGACTACTATCGCTTTGACCAGGTAGGGCAACGCCTGATGGGGGAATCCAGCGGCCA  >  W3110S.gb/4413308‑4413428

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: