Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,215,186 A→G 69.0% V115V (GTT→GTC yehS ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,215,1860AG69.0% 37.9 / 20.8 29V115V (GTT→GTCyehSconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base G (10/10);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

GCCGCATTCGCGGAAGTTTTTATGATCCGGTGCGCGCATCATCGCTGTAATTTCCGGCATCGAAACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTTTTCAGCGAAAACGCGATGC  >  W3110S.gb/2215123‑2215250
                                                               |                                                                
gccgcATTCGCGGAAGTTTTTATGATCCGGTGCGCGCATCATCGCTGTAATTTCCGGCATCGAGACGCgg                                                            <  1:1574317/70‑1 (MQ=255)
      ttCGCGGAAGTTTTTATGATCCGGTGCGCGCATCATCGCTGTAATTTCCGGCATCGAGAc                                                                >  1:1955770/1‑60 (MQ=255)
        cgcgGAAGTTTTTATGATCCGGTGCGCGCATCATCGCTGTAATTTCCGGCATCGAGACGCGGaa                                                          >  1:2191521/1‑64 (MQ=255)
               gTTTTTATGATCCGGTGCGCGCATCATCGCTGTAATTTCCGGCATCGAGACGCgg                                                            <  1:163591/55‑1 (MQ=255)
                tttttATGATCCGGTGCGCGCATCATCGCTGTAATTTCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAg                                           <  1:1247567/71‑1 (MQ=255)
                tttttATGATCCGGTGCGCGCATCATCGCTGTAATTTCCGGCATCGAAACGCGGAACTGCTGTTCGGTGAg                                           <  1:271038/71‑1 (MQ=255)
                tttttATGATCCGGTGCGCGCATCATCGCTGTAATTTCCGGCATCGAAACGCGGAACTGCTGTTCGGTGAg                                           <  1:2675579/71‑1 (MQ=255)
                         tCCGGTGCGCGCATCATCGCTGTAATTTCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGATAGCCAg                                  <  1:2588662/71‑1 (MQ=255)
                           cGGTGCGCGCATCATCGCTGTAATTTCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGATAGCCAGa                                 <  1:2309272/70‑1 (MQ=255)
                                      tcatcGCTGTAATTTCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGAt                                        <  1:91455/52‑1 (MQ=255)
                                      tcatcGCTGTAATTTCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGAt                                        <  1:1143553/52‑1 (MQ=255)
                                        atcGCTGTAATTTCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGATAg                                      >  1:559146/1‑52 (MQ=255)
                                        atcGCTGTAATTTCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGATAg                                      >  1:2575529/1‑52 (MQ=255)
                                        atcGCTGTAATTTCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGtttt                   >  1:1927283/1‑71 (MQ=255)
                                        atcGCTGTAATTTCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGtttt                   >  1:1105059/1‑71 (MQ=255)
                                        atcGCTGTAATTTCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGtttt                   >  1:599776/1‑71 (MQ=255)
                                           gCTGTAATTTCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTTTTCAg                >  1:895998/1‑71 (MQ=255)
                                           gCTGTAATTTCCGGCATCGAAACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTTTTCAg                >  1:2093344/1‑71 (MQ=255)
                                           gCTGTAATTTCCGGCATCGAAACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTTTTCAg                >  1:2703987/1‑71 (MQ=255)
                                           gCTGTAATTTCCGGCATCGAAACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTGTTCAg                >  1:2179112/1‑71 (MQ=255)
                                               tAATTTCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTTTTCAGCGaa            <  1:2638125/71‑1 (MQ=255)
                                               tAATTTCCGGCATCGAAACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTTTTCAGCGaa            <  1:459154/71‑1 (MQ=255)
                                                  tttCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGa                                         <  1:812888/39‑1 (MQ=255)
                                                    tCCGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATCTCATCGGTTTTCAGCGAAAAcg        >  1:1838497/1‑70 (MQ=255)
                                                     ccGGCATCGAGACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTTTTCAGCGAAAAcgc       <  1:2637685/70‑1 (MQ=255)
                                                     ccGGCATCGAAACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTTTTCAGCGAAAAcgc       <  1:1268902/70‑1 (MQ=255)
                                                         aaTCGAGACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTTTTCAGCGAAAACGCGATGc  >  1:1478497/2‑71 (MQ=255)
                                                           tCGAAACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTTTTCAGCGAAAAcgc       <  1:1681225/64‑1 (MQ=255)
                                                             gAAACTCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTTTTCAGCGa             >  1:2401636/1‑56 (MQ=255)
                                                               |                                                                
GCCGCATTCGCGGAAGTTTTTATGATCCGGTGCGCGCATCATCGCTGTAATTTCCGGCATCGAAACGCGGAACTGCTGTTCGGTGAGGATAGCCAGAATGTCATCGGTTTTCAGCGAAAACGCGATGC  >  W3110S.gb/2215123‑2215250

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: