Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,669,8380TC73.9% 78.9 / 28.6 46V92A (GTT→GCT) hcaC3‑phenylpropionate dioxygenase, predicted ferredoxin subunit
Reads supporting (aligned to +/- strand):  ref base T (0/12);  new base C (22/12);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 8.61e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCTGAAAACGGGGAAAGCGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTC  >  W3110S.gb/2669772‑2669903
                                                                  |                                                                 
gCCTGAAAACGGGGAAAGCGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGCTGa                                                                <  1:1908152/70‑1 (MQ=255)
 ccTGAAAACGGGGAAAGCGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGCTGaa                                                               <  1:3796274/70‑1 (MQ=255)
 ccTGAAAACGGGGAAAGCGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGCTGaa                                                               <  1:2259952/70‑1 (MQ=255)
 ccTGAAAACGGGGAAAGCGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGCTGaa                                                               <  1:2799659/70‑1 (MQ=255)
                       tgcctgccCGCCACCGATCCGCTCACTACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACt                                         <  1:3817623/70‑1 (MQ=255)
                       tgcctgccCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACt                                         <  1:534984/70‑1 (MQ=255)
                       tgcctgccCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACt                                         <  1:2442821/70‑1 (MQ=255)
                       tgcctgccCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACt                                         <  1:1755922/70‑1 (MQ=255)
                       tgcctgccCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACt                                         <  1:3401237/70‑1 (MQ=255)
                       tgcctgccCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACt                                         <  1:3360360/70‑1 (MQ=255)
                       tgcctgccCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACt                                         <  1:2607573/70‑1 (MQ=255)
                       tgcctgccCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACt                                         <  1:3088771/70‑1 (MQ=255)
                       tgcctgccCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACt                                         <  1:2203897/70‑1 (MQ=255)
                       tgcctgccCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACt                                         <  1:2193427/70‑1 (MQ=255)
                             cccGCCACCGATCCGCTCACCACTTATCCAGTACACGt                                                                   >  1:2270505/1‑38 (MQ=255)
                             cccGCCACCGATCCGCTCACCACTTATCCAGTACACGt                                                                   >  1:3009812/1‑38 (MQ=255)
                             cccGCCACCGATCCGCTCACCACTTATCCAGTACACGt                                                                   >  1:2930790/1‑38 (MQ=255)
                             cccGCCACCGATCCGCTCACCACTTATCCAGTACACGt                                                                   >  1:3327903/1‑38 (MQ=255)
                             cccGCCACCGATCCGCTCACCACTTATCCAGTACACGt                                                                   >  1:3455288/1‑38 (MQ=255)
                             cccGCCACCGATCCGCTCACCACTTATCCAGTACACGt                                                                   >  1:4004401/1‑38 (MQ=255)
                             cccGCCACCGATCCGCTCACCACTTATCCAGTACACGt                                                                   >  1:876287/1‑38 (MQ=255)
                             cccGCCACCGATCCGCTCACCACTTATCCAGTACACGt                                                                   >  1:2265758/1‑38 (MQ=255)
                             cccGCCACCGATCCGCTCACCACTTATCCAGTACACGt                                                                   >  1:1363878/1‑38 (MQ=255)
                             cccGCCACCGATCCGCTCACCACTTATCCAGTACACGt                                                                   >  1:20818/1‑38 (MQ=255)
                             cccGCCACCGATCCGCTCACCACTTATCCAGTACACGt                                                                   >  1:1823328/1‑38 (MQ=255)
                             cccGCCACCGATCCGCTCACCACTTATCCAGTACACGt                                                                   >  1:1983743/1‑38 (MQ=255)
                                                  tcTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACt                                         <  1:2366959/42‑1 (MQ=255)
                                                        ccAGTACACGTTGAAGGTGGTGACATTTTCATCGACt                                         <  1:3640868/37‑1 (MQ=255)
                                                         cAGTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAAc       >  1:3854276/1‑70 (MQ=255)
                                                         cAGTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAAc       >  1:3650331/1‑70 (MQ=255)
                                                         cAGTACACGCTGAAGGTGGTGACATTTTCATCGACTTACAGGAGGCGCAGCCATGAGCGATCTGCATaa        >  1:214806/1‑69 (MQ=255)
                                                           gTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATc               <  1:1231258/60‑1 (MQ=255)
                                                           gTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATc               <  1:863905/60‑1 (MQ=255)
                                                           gTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATc               <  1:820761/60‑1 (MQ=255)
                                                           gTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATc               <  1:3712314/60‑1 (MQ=255)
                                                           gTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATc               <  1:1835710/60‑1 (MQ=255)
                                                           gTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATc               <  1:2057878/60‑1 (MQ=255)
                                                           gTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATc               <  1:2612117/60‑1 (MQ=255)
                                                           gTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATc               <  1:1054863/60‑1 (MQ=255)
                                                           gTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAg    >  1:4091665/1‑71 (MQ=255)
                                                           gTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAg    >  1:3072119/1‑71 (MQ=255)
                                                            tACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGt   >  1:2262904/1‑71 (MQ=255)
                                                            tACACGCTGAAGGTGGGGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGt   >  1:1359913/1‑71 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGt   >  1:3142816/1‑69 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc  >  1:2067188/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc  >  1:3806997/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc  >  1:1789428/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc  >  1:3639325/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc  >  1:3874976/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc  >  1:2749548/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc  >  1:3631622/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc  >  1:487866/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc  >  1:60542/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc  >  1:323010/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc  >  1:2223342/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc  >  1:1984016/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAAAGAGTc  >  1:3703675/1‑70 (MQ=255)
                                                              cacGCTGAAGGTGGTGACATTTTCATCGACTTACAGGAGGCGCAGCCATGAGCGATCTGCATAACGAg    >  1:3037496/1‑68 (MQ=255)
                                                                  |                                                                 
GCCTGAAAACGGGGAAAGCGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTC  >  W3110S.gb/2669772‑2669903

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: