Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,667,558 T→C 35.5% intergenic (‑6/‑130) hcaR ← / → hcaE DNA‑binding transcriptional regulator/3‑phenylpropionate dioxygenase, large (alpha) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,667,5580TC35.5% 24.7 / 27.8 31intergenic (‑6/‑130)hcaR/hcaEDNA‑binding transcriptional regulator/3‑phenylpropionate dioxygenase, large (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base T (15/5);  new base C (9/2);  total (24/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAAC  >  W3110S.gb/2667521‑2667623
                                     |                                                                 
gCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTgg                                    >  1:999606/1‑69 (MQ=255)
gCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTgg                                    >  1:557622/1‑69 (MQ=255)
gCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTgg                                    >  1:499460/1‑69 (MQ=255)
gCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTgg                                    >  1:2902720/1‑69 (MQ=255)
           gCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTa                       >  1:1429895/1‑71 (MQ=255)
           gCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTa                       >  1:770986/1‑71 (MQ=255)
           gCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTa                       >  1:3173737/1‑71 (MQ=255)
           gCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTa                       >  1:3102082/1‑71 (MQ=255)
           gCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTa                       >  1:2978219/1‑71 (MQ=255)
           gCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTa                       >  1:283942/1‑71 (MQ=255)
           gCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTa                       >  1:2738448/1‑71 (MQ=255)
           gCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTa                       >  1:2149295/1‑71 (MQ=255)
           gCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTa                       >  1:2511653/1‑71 (MQ=255)
           gCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTa                       >  1:2064567/1‑71 (MQ=255)
           gCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTa                       >  1:1550961/1‑71 (MQ=255)
              tAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAttt                    <  1:2638686/71‑1 (MQ=255)
              tAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAttt                    <  1:2765212/71‑1 (MQ=255)
              tAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAttt                    <  1:2855405/71‑1 (MQ=255)
              tAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAttt                    <  1:300377/71‑1 (MQ=255)
              tAAATGCCGTAGTTCCATCACCTTACCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAttt                    <  1:150420/71‑1 (MQ=255)
                aaTGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTc                   >  1:2030463/1‑70 (MQ=255)
                  tGCCGTAGTTCCATCACCTCCCCCTTGTTATGGAAAAAACGTCTCAGCTGGTAGTAATTAACTATTTcaca                <  1:2634774/71‑1 (MQ=255)
                  tGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTcaca                <  1:1614241/71‑1 (MQ=255)
                                 aCCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAttt                    >  1:787229/1‑52 (MQ=255)
                                 aCCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGcaaccaac  >  1:1050881/1‑70 (MQ=255)
                                 aCCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGcaaccaac  >  1:853200/1‑70 (MQ=255)
                                 aCCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGcaaccaac  >  1:163248/1‑70 (MQ=255)
                                 aCCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGcaaccaac  >  1:1706014/1‑70 (MQ=255)
                                 aCCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGcaaccaac  >  1:2017733/1‑70 (MQ=255)
                                 aCCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGcaaccaac  >  1:2877533/1‑70 (MQ=255)
                                 aCCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGcaaccaac  >  1:2771799/1‑70 (MQ=255)
                                     |                                                                 
GCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAAC  >  W3110S.gb/2667521‑2667623

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: