Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,971,4010AT58.8% 30.2 / 62.3 51V373E (GTG→GAG) ygeDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (13/8);  new base T (26/4);  total (39/12)
Fisher's exact test for biased strand distribution p-value = 5.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGATCCACAGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGTTT  >  W3110S.gb/2971341‑2971466
                                                            |                                                                 
cAGATCCACAGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGAt                                                         <  1:1091753/71‑1 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTa                                                    >  1:2873595/1‑70 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTa                                                    >  1:3140667/1‑70 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTAt                                                   >  1:2204382/1‑71 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTAt                                                   >  1:2340227/1‑71 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTAt                                                   >  1:2843753/1‑71 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTAt                                                   >  1:1720646/1‑71 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTAt                                                   >  1:3038463/1‑71 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTAt                                                   >  1:3371528/1‑71 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTAt                                                   >  1:138487/1‑71 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTAt                                                   >  1:794484/1‑71 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTAt                                                   >  1:900635/1‑71 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTAt                                                   >  1:1092200/1‑71 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTAt                                                   >  1:3861666/1‑71 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTAt                                                   >  1:3942333/1‑71 (MQ=255)
      cacaGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTAt                                                   >  1:2582395/1‑71 (MQ=255)
         aGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCt                                                     <  1:807734/66‑1 (MQ=255)
           cgcCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCt                                                     <  1:2765725/64‑1 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:2114114/1‑50 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:521076/1‑50 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:243104/1‑50 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:1405784/1‑50 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:2753841/1‑50 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:1886362/1‑50 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:1140471/1‑50 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:140597/1‑50 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:3024423/1‑50 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:3026780/1‑50 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:3731483/1‑50 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:3182353/1‑50 (MQ=255)
                  aTTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                            >  1:1561656/1‑50 (MQ=255)
                     gCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCATTACCGCCAGc                                      >  1:1893124/1‑69 (MQ=255)
                        aGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTa                                 <  1:3122483/71‑1 (MQ=255)
                        aGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCATTACCGACAGCGAGTa                                 <  1:2420065/71‑1 (MQ=255)
                           gcAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGAt                                                         <  1:1484212/44‑1 (MQ=255)
                           gcAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTATCAt                                                >  1:3738723/1‑53 (MQ=255)
                           gcAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTATCAt                                                >  1:24695/1‑53 (MQ=255)
                                 aGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTATCAt                                                >  1:488325/1‑47 (MQ=255)
                                                      atgGGCACGACCTGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGtt   <  1:3606316/71‑1 (MQ=255)
                                                      atgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGtt   <  1:3084868/71‑1 (MQ=255)
                                                      atgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGtt   <  1:1747444/71‑1 (MQ=255)
                                                      atgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGtt   <  1:3966599/71‑1 (MQ=255)
                                                      atgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGtt   <  1:712535/71‑1 (MQ=255)
                                                      atgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGtt   <  1:1286236/71‑1 (MQ=255)
                                                       tgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGttt  >  1:971615/1‑71 (MQ=255)
                                                       tgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGttt  >  1:1231710/1‑71 (MQ=255)
                                                       tgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGttt  >  1:89567/1‑71 (MQ=255)
                                                       tgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGttt  >  1:3407346/1‑71 (MQ=255)
                                                       tgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGttt  >  1:1602986/1‑71 (MQ=255)
                                                       tgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGttt  >  1:1816348/1‑71 (MQ=255)
                                                       tgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGttt  >  1:1896677/1‑71 (MQ=255)
                                                            |                                                                 
CAGATCCACAGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTGTTT  >  W3110S.gb/2971341‑2971466

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: