Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,669,838 T→C 75.0% V92A (GTT→GCT)  hcaC → 3‑phenylpropionate dioxygenase, predicted ferredoxin subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,669,8380TC75.0% 32.8 / 8.6 20V92A (GTT→GCT) hcaC3‑phenylpropionate dioxygenase, predicted ferredoxin subunit
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (7/8);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGCCTGAAAACGGGGAAAGCGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTCCATTT  >  W3110S.gb/2669769‑2669908
                                                                     |                                                                      
tttGCCTGAAAACGGGGAAAGCGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGtt                                                                       >  1:421960/1‑71 (MQ=255)
   gCCTGAAAACGGGGAAAGCGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGa                                                                     <  1:9109/70‑1 (MQ=255)
        aaaaCGGGGAAAGCGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAggt                                                                 <  1:137010/69‑1 (MQ=255)
                   aGCGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGCTGAAGGTGGTGACATTTTCa                                                    >  1:309986/1‑71 (MQ=255)
                     cGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGCTGAAGGTGGTGACATTTTCATc                                                  <  1:1596999/71‑1 (MQ=255)
                         aTGCCTTCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGCTGAAgg                                                                  >  1:640888/1‑51 (MQ=39)
                         aTGCCTGCCCGCCACCGATCCGCTCACCGCTTATCCAGTACACGCTGAAgg                                                                  >  1:2117992/1‑51 (MQ=39)
                         aTGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGCTGAAgg                                                                  >  1:1613491/1‑51 (MQ=255)
                         aTGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGCTGAAgg                                                                  >  1:2367256/1‑51 (MQ=255)
                                                  accacTTATCCAGTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCg                       <  1:3133527/69‑1 (MQ=255)
                                                   ccacTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATc                    >  1:1144891/1‑71 (MQ=255)
                                                   ccacTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATc                    >  1:267108/1‑71 (MQ=255)
                                                   ccacTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATc                    >  1:3547826/1‑71 (MQ=255)
                                                            cAGTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACg           <  1:1892791/71‑1 (MQ=255)
                                                             aGTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACg           <  1:995138/70‑1 (MQ=255)
                                                             aGTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACg           <  1:610129/70‑1 (MQ=255)
                                                             aGTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACg           <  1:2818797/70‑1 (MQ=255)
                                                             aGTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACg           <  1:138220/70‑1 (MQ=255)
                                                             aGTACACGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGa          <  1:2254248/71‑1 (MQ=255)
                                                                 cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGt        >  1:2642853/1‑69 (MQ=255)
                                                                 cacGCTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTc       >  1:454694/1‑70 (MQ=255)
                                                                     ttGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTCCAttt  >  1:3363270/1‑71 (MQ=255)
                                                                     |                                                                      
TTTGCCTGAAAACGGGGAAAGCGTTATGCCTGCCCGCCACCGATCCGCTCACCACTTATCCAGTACACGTTGAAGGTGGTGACATTTTCATCGACTTACCGGAGGCGCAGCCATGAGCGATCTGCATAACGAGTCCATTT  >  W3110S.gb/2669769‑2669908

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: