Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,389,988 C→T 32.3% intergenic (+13/+31) tyrR → / ← tpx DNA‑binding transcriptional dual regulator, tyrosine‑binding/lipid hydroperoxide peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,389,9880CT32.3% 40.8 / 26.6 34intergenic (+13/+31)tyrR/tpxDNA‑binding transcriptional dual regulator, tyrosine‑binding/lipid hydroperoxide peroxidase
Reads supporting (aligned to +/- strand):  ref base C (8/15);  new base T (9/2);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 2.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAA  >  W3110S.gb/1389923‑1390052
                                                                 |                                                                
ttGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCTTGa                                                               <  1:2314200/69‑1 (MQ=255)
ttGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCt                                                                 <  1:146625/67‑1 (MQ=255)
  gCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCTTGATg                                                             >  1:3788806/1‑69 (MQ=255)
  gCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCTTGATGGt                                                           >  1:1169144/1‑71 (MQ=255)
      aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTg                                                                <  1:3419913/62‑1 (MQ=255)
      aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTg                                                                <  1:1618379/62‑1 (MQ=255)
      aTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTg                                                                <  1:397702/62‑1 (MQ=255)
       tAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCaa                                                       >  1:561649/1‑70 (MQ=255)
        aaGTTGCGGGAATATGGTCTGAGTCATAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAaca                                                     >  1:3694283/1‑71 (MQ=255)
             gCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCAt                                                 >  1:3315133/1‑70 (MQ=255)
             gCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCAt                                                 >  1:412208/1‑70 (MQ=255)
              cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATc                                                <  1:566893/70‑1 (MQ=255)
              cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATc                                                <  1:2226273/70‑1 (MQ=255)
              cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATc                                                <  1:2104557/70‑1 (MQ=255)
              cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATc                                                <  1:3429665/70‑1 (MQ=255)
              cGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATc                                                <  1:1276897/70‑1 (MQ=255)
                ggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCTTGATGGTGCAACACCATCAgg                                             <  1:821391/71‑1 (MQ=255)
                  aaTATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCTTGATGGTGCAACACCATCAGGCa                                           >  1:2800912/1‑71 (MQ=255)
                      tGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCTTGATGGTGCAACACCATCAGGCatat                                        >  1:3027274/1‑70 (MQ=255)
                      tGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCatat                                        >  1:1875353/1‑70 (MQ=255)
                        gTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCTTGATGGTGCAACACCATCAGGCATATTa                                      >  1:431544/1‑70 (MQ=255)
                             aGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTa                                 <  1:297282/70‑1 (MQ=255)
                                 agaagaagaACGAAGAGTAAGCGCGAATATGCTTGATGGTGCAACACCATCAGGCATATTAAATTATGCtt                            >  1:2607872/1‑71 (MQ=255)
                                   aagaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTa                                 >  1:3599538/1‑64 (MQ=255)
                                    agaagaaCGAAGAGTAAGCGCGAATATGCTTGATGGTGCAACACCATCAGGCATATTAAATTATGCttt                           >  1:3263981/1‑69 (MQ=255)
                                    agaagaaCGAAGAGTAAGCGCGAATATGCTTGATGGTGCAACACCATCAGGCATATTAAATTATGCttt                           >  1:2653508/1‑69 (MQ=255)
                                     gaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCAt                                                 <  1:3269056/46‑1 (MQ=255)
                                     gaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCAt                                                 <  1:2887348/46‑1 (MQ=255)
                                     gaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCAt                                                 <  1:3050264/46‑1 (MQ=255)
                                     gaagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAg                        <  1:2753129/71‑1 (MQ=255)
                                        gaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTa                      >  1:4026774/1‑70 (MQ=255)
                                             aaGAGTAAGCGCGAATATGCTTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGc                  >  1:2256118/1‑69 (MQ=255)
                                             aaGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCa                                           >  1:487339/1‑44 (MQ=255)
                                                           atatGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTaa  <  1:2979439/71‑1 (MQ=255)
                                                                 |                                                                
TTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAA  >  W3110S.gb/1389923‑1390052

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: