Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,495,9661.T62.0% 82.6 / 71.6 50coding (820/924 nt)ydcIpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (3/16);  new base T (17/14);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 7.97e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAG  >  W3110S.gb/1495896‑1496025
                                                                       |                                                           
atcaACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATg                                                              >  1:2199324/1‑71 (MQ=255)
  caACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGtt                                                           >  1:2672484/1‑70 (MQ=255)
  caACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGtt                                                           >  1:4100722/1‑70 (MQ=255)
  caACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGtt                                                           >  1:2266780/1‑70 (MQ=255)
          aaCCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATGGGTTCGCCTGCACCATGTTCCCGGaa                                                     >  1:3985329/1‑70 (MQ=255)
          aaCCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGaa                                                     >  1:2266169/1‑70 (MQ=255)
                             aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg                                 >  1:3323699/1‑71 (MQ=255)
                             aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg                                 >  1:202032/1‑71 (MQ=255)
                             aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg                                 >  1:2832614/1‑71 (MQ=255)
                             aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg                                 >  1:1515927/1‑71 (MQ=255)
                             aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg                                 >  1:3623910/1‑71 (MQ=255)
                             aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg                                 >  1:2286819/1‑71 (MQ=255)
                             aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg                                 >  1:2659274/1‑71 (MQ=255)
                             aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTgg                                >  1:2814217/1‑71 (MQ=255)
                             aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTgg                                >  1:2990393/1‑71 (MQ=255)
                             aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGACACCAGCGTgg                                >  1:3342588/1‑71 (MQ=255)
                                tACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGc                               >  1:604469/1‑70 (MQ=255)
                                tACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCa                              >  1:3215022/1‑71 (MQ=255)
                                tACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCa                              >  1:3073073/1‑71 (MQ=255)
                                tACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCa                              >  1:1363309/1‑71 (MQ=255)
                                tACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCa                              >  1:636618/1‑71 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTgga                               <  1:2442608/68‑2 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGc                               <  1:1536046/68‑1 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGc                               <  1:649081/68‑1 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAAGAGGCAGGGCCACCAGCGTGGCAt                             <  1:851075/69‑1 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCAt                             <  1:2682729/69‑1 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCAt                             <  1:820210/69‑1 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCAt                             <  1:1246941/69‑1 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCAt                             <  1:1286454/69‑1 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCAt                             <  1:3988610/69‑1 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCAt                             <  1:1388366/69‑1 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCAt                             <  1:3656852/69‑1 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCAt                             <  1:1430965/69‑1 (MQ=255)
                                  cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCAt                             <  1:1431657/69‑1 (MQ=255)
                                      ggTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACg                         >  1:3914715/1‑70 (MQ=255)
                                           gTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtc                   <  1:2823350/71‑1 (MQ=255)
                                           gTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtc                   <  1:1697375/71‑1 (MQ=255)
                                            tATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtcg                  <  1:1219474/71‑1 (MQ=255)
                                            tATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtcg                  <  1:2243241/71‑1 (MQ=255)
                                            tATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc                        <  1:1581956/64‑1 (MQ=255)
                                            tATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc                        <  1:1008290/64‑1 (MQ=255)
                                            tATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc                        <  1:537729/64‑1 (MQ=255)
                                            tATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc                        <  1:3639503/64‑1 (MQ=255)
                                            tATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc                        <  1:3005788/64‑1 (MQ=255)
                                                 cAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt             <  1:3903870/71‑1 (MQ=255)
                                                 cAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt             <  1:286187/71‑1 (MQ=255)
                                                 cAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt             <  1:1841289/71‑1 (MQ=255)
                                                 cAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt             <  1:3000188/71‑1 (MQ=255)
                                                 cAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt             <  1:3381154/71‑1 (MQ=255)
                                                 cAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt             <  1:1239624/71‑1 (MQ=255)
                                                           gCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc                        <  1:1610177/49‑1 (MQ=255)
                                                             cTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCt               >  1:1483486/1‑57 (MQ=255)
                                                             cTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAg  >  1:2453602/1‑70 (MQ=255)
                                                                 aCCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGc                               <  1:679084/37‑1 (MQ=37)
                                                                       |                                                           
ATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAG  >  W3110S.gb/1495896‑1496025

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: