Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,855,731 G→T 17.8% intergenic (‑269/‑18) ygbA ← / → mutS hypothetical protein/methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,855,7310GT17.8% 85.4 / 5.0 28intergenic (‑269/‑18)ygbA/mutShypothetical protein/methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (10/13);  new base T (2/3);  total (12/16)
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.85e-01

CATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCT  >  W3110S.gb/2855671‑2855801
                                                            |                                                                      
tctcgatctggtgcagaagaaaggcggaagttgcctctggcgagcataaTTTAATATCAGTAAACCGGACa                                                              >  1:791875‑M2/50‑71 (MQ=255)
  tgtgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATa                                                            <  1:3135383/71‑1 (MQ=255)
   gtgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATaa                                                           >  1:2319088/1‑71 (MQ=255)
      cgcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCC                                                        >  1:1059870‑M1/1‑71 (MQ=255)
        ccTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGAc                                                               <  1:57209/62‑1 (MQ=255)
         cTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCAt                                                     <  1:2964988/71‑1 (MQ=255)
          gtgcagaagaaaggcggaagttgcctctggcgagcgtaaTTTAATATCAGTAAACCGGACATAAcc                                                         <  1:1471675‑M2/27‑1 (MQ=255)
           tATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCC                                                        <  1:2904546‑M1/66‑1 (MQ=255)
              agaagaaaggcggaagttgcctctggcgagcataaTTTAATATCAGTAAACCGGACATAAc                                                          <  1:957376‑M2/26‑1 (MQ=255)
                 aTTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                             <  1:1140891/71‑1 (MQ=255)
                   tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                             >  1:2200929/1‑69 (MQ=255)
                   tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCaa                                             >  1:1553681/1‑69 (MQ=255)
                      aaCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAAc                                                                    <  1:3071333/43‑1 (MQ=255)
                         gAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGaaa                                       >  1:50747/1‑69 (MQ=255)
                            aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACAt                                                             <  1:953858/44‑1 (MQ=255)
                             aTAAAAACCATCACACCCCATTTAATATCAGGGAAc                                                                    <  1:1369009/36‑1 (MQ=255)
                             aTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCAt                                                     <  1:1329413/51‑1 (MQ=255)
                                         acacCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgc                     <  1:1161956/71‑1 (MQ=255)
                                              ccATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgcccat                 >  1:2280436/1‑70 (MQ=255)
                                              ccATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGacgcccatacgcccat                 >  1:238663/1‑70 (MQ=255)
                                               cATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatga               <  1:315788/71‑1 (MQ=255)
                                                   tAATATCAGTAAACCGGACATAACCCatgaagaaaaatcgcgcttttttgaagtgggcagggggcaagta            <  1:1014154‑M1/70‑45 (MQ=255)
                                                   tAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatgatg             >  1:499699/1‑69 (MQ=255)
                                                      tatCAGTAAACCGGACATAACCCatgaagaaaaatc                                           >  1:159597‑M1/1‑23 (MQ=255)
                                                      tatCAGGGAACCGGACATAACCCCATGAGTGCAATAGaa                                        <  1:1000298/39‑1 (MQ=255)
                                                        tCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCatgatg             >  1:1649716/1‑64 (MQ=255)
                                                          aGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTc                                  <  1:356251/41‑1 (MQ=255)
                                                          aGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGa                                >  1:1738932/1‑43 (MQ=255)
                                                            ggAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTAtct  <  1:2713093/71‑1 (MQ=255)
                                                            |                                                                      
CATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCT  >  W3110S.gb/2855671‑2855801

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: