Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,885,9840CT53.6% 8.7 / 99.6 69intergenic (‑109/+250)ygcB/cysHconserved hypothetical protein, member of DEAD box family/3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base C (17/15);  new base T (23/14);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGA  >  W3110S.gb/2885928‑2886042
                                                        |                                                          
tCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTa                                               <  1:3122424/70‑1 (MQ=255)
tCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTa                                               <  1:3299920/70‑1 (MQ=255)
  gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAagaaga                                                    <  1:1627832/63‑1 (MQ=255)
      ctgctCATAAATATCGTATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:3961864/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAATAGACTAACATTc                                         >  1:904790/1‑70 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:2184187/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:1021587/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:3232935/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:3007955/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:3300946/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:2738422/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:3335559/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:3561323/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:2237636/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:2179759/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:2041625/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:966422/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:1393537/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:639841/1‑61 (MQ=255)
      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:624199/1‑61 (MQ=255)
                 tatCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTggg                             <  1:99170/71‑1 (MQ=255)
                 tatCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTggg                             <  1:3538124/71‑1 (MQ=255)
                 tatCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTggg                             <  1:1683514/71‑1 (MQ=255)
                 tatCGAATTACAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTggg                             <  1:2188497/71‑1 (MQ=255)
                          aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCg                         >  1:3300591/1‑66 (MQ=255)
                          aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGc                    >  1:247981/1‑71 (MQ=255)
                          aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGc                    >  1:4007275/1‑71 (MQ=255)
                          aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGc                    >  1:677941/1‑71 (MQ=255)
                          aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGc                    >  1:3243433/1‑71 (MQ=255)
                          aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGc                    >  1:982114/1‑71 (MQ=255)
                          aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGc                    >  1:1229351/1‑71 (MQ=255)
                          aCAGTGATATACACATGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGc                    >  1:488977/1‑71 (MQ=255)
                                atatACACATGTCTTGCCTCCCTGCAAAGAAgaga                                                  >  1:772538/1‑35 (MQ=255)
                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTc             <  1:63443/69‑1 (MQ=255)
                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTc             <  1:3213522/69‑1 (MQ=255)
                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTc             <  1:2968986/69‑1 (MQ=255)
                                     cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGt         >  1:2743928/1‑71 (MQ=255)
                                     cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGt         >  1:2569785/1‑71 (MQ=255)
                                     cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGt         >  1:3998488/1‑71 (MQ=255)
                                     cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGt         >  1:1636486/1‑71 (MQ=255)
                                     cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGt         >  1:3773277/1‑71 (MQ=255)
                                     cacaTGTCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGt         >  1:1389389/1‑71 (MQ=255)
                                           tCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGa      >  1:2847020/1‑68 (MQ=255)
                                           tCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGa      >  1:1257055/1‑68 (MQ=255)
                                           tCTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGa      >  1:3383953/1‑68 (MQ=255)
                                            cTTGCCTCCCTGTAAAGAAGAGACTAACATTCGGATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg   >  1:1880031/1‑70 (MQ=255)
                                            cTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTg       >  1:93743/1‑66 (MQ=255)
                                            cTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg   >  1:883920/1‑70 (MQ=255)
                                            cTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg   >  1:1615964/1‑70 (MQ=255)
                                            cTTGCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg   >  1:1839258/1‑70 (MQ=255)
                                               gCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACTCTGGCCTCAGGTTGATAGa  <  1:127457/68‑1 (MQ=255)
                                               gCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa  <  1:4081611/68‑1 (MQ=255)
                                               gCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa  <  1:1124730/68‑1 (MQ=255)
                                               gCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa  <  1:3136816/68‑1 (MQ=255)
                                               gCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa  <  1:2546636/68‑1 (MQ=255)
                                               gCCTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa  <  1:3092452/68‑1 (MQ=255)
                                                ccTCCCTGTAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa  <  1:2223481/67‑1 (MQ=255)
                                                      tGCAAAGAAGAGACTAACATTCGAATTTTTTGGGTTCGAACGCTGGCCTCa            <  1:803145/51‑1 (MQ=255)
                                                      tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCa            <  1:421973/51‑1 (MQ=255)
                                                      tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCa            <  1:3969835/51‑1 (MQ=255)
                                                      tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCa            <  1:244220/51‑1 (MQ=255)
                                                      tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCa            <  1:1492944/51‑1 (MQ=255)
                                                      tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCa            <  1:3621954/51‑1 (MQ=255)
                                                      tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCa            <  1:3602500/51‑1 (MQ=255)
                                                      tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCa            <  1:128954/51‑1 (MQ=255)
                                                      tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCa            <  1:3566153/51‑1 (MQ=255)
                                                      tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCa            <  1:967778/51‑1 (MQ=255)
                                                      tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCa            <  1:2677797/51‑1 (MQ=255)
                                                      tGCAAAGAAGAGACTAACATGCGAATTGTTTGGGTTCGAACGCTGGCCTCa            <  1:821246/51‑1 (MQ=255)
                                                        |                                                          
TCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGA  >  W3110S.gb/2885928‑2886042

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: