Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 251,041 G→A 25.9% A48A (GCG→GCA dinB → DNA polymerase IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb251,0410GA25.9% 45.8 / 13.7 27A48A (GCG→GCAdinBDNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base G (8/12);  new base A (3/4);  total (11/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.74e-01

TCCCTATTGCTATTGGCGGCAGCCGCGAACGTCGGGGGGTGATCAGCACCGCCAATTATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGACAGGGATGGCGCTCAAATTATGCCCACATCT  >  W3110S.gb/250977‑251103
                                                                |                                                              
tCCCTATTGCTATTGGCGGCAGCCGCGAACGTCGGGGGGTGATCAGCACCGCCAATTATCCCGCACGTaa                                                           <  1:1410284/70‑1 (MQ=255)
      ttgctattgGCGGCAGCCGCGAACGTCGGGGGGTGATCAGCACCGCCAATTATCCCGCGCGTAAAt                                                         >  1:558658/1‑66 (MQ=255)
      ttgctattgGCGGCAGCCGCGAACGTCGGGGGGTGATCAGCACCGCCAATTATCCCGCACGTAAATTTg                                                      <  1:624014/69‑1 (MQ=255)
        gctattgGCGGCAGCCGCGAACGTCGGGGGGTGATCAGCACCGCCAATTATCCCGCGCGTAAATTTGGcg                                                   >  1:1520350/1‑70 (MQ=255)
              ggcggcAGCCGCGAACGTCGGGGGGTGATCAGCACCGCCAATTATCCCGCGCGTAAATTTGGCGTACGTAg                                            >  1:481870/1‑71 (MQ=255)
                 ggcAGCCGCGAACGTCGGGGGGTGATCAGCACCGCCAATTATCCCGCACGTAAATTTGGCGTACGTAgcgc                                         <  1:800192/71‑1 (MQ=255)
                             cGTCGGGGGGTGATCAGCACCGCCAATTATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGACAgg                             <  1:1141042/71‑1 (MQ=255)
                               tCGGGGGGTGATCAGCACCGCCAATTATCCCGCGCGTAAATTTgg                                                     <  1:1249130/45‑1 (MQ=255)
                                cGGGGGGTGATCAGCACCGCCAATTATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGACAGGGAt                          >  1:804443/1‑71 (MQ=255)
                                   ggggTGATCAGCACCGCCAATTATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGACa                               <  1:1116052/63‑1 (MQ=255)
                                   ggggTGATCAGCACCGCCAATTATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGACa                               <  1:101868/63‑1 (MQ=255)
                                     ggTGATCAGCACCGCCAATTATCCCGCGCGTAAATTTGGcgtacg                                               <  1:1627697/45‑1 (MQ=255)
                                       tGATCAGCACCGCCAATTATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGACAgg                             <  1:624288/61‑1 (MQ=255)
                                       tGATCAGCACCGCCAATTATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGACAgg                             <  1:1186550/61‑1 (MQ=255)
                                       tGATCAGCACCGCCAATTATCCCGCACGTAAATTTGGCGTACGTAgcgc                                         <  1:1737502/49‑1 (MQ=255)
                                         aTCAGCACCGCCAATTATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGa                                 <  1:143767/55‑1 (MQ=255)
                                         aTCAGCACCGCCAATTATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGa                                 <  1:1734810/55‑1 (MQ=255)
                                         aTCAGCACCGCCAATTATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGa                                 <  1:693129/55‑1 (MQ=255)
                                            aGCACCGCCAATTATCCCGCACGTAAATTTGGCGTACGTAGCGCTATGCCGACAGGGATGGCGCTCaaa                >  1:1553435/1‑69 (MQ=255)
                                                  gCCAATTATCCCGCGCGTAAATTTGGCGTACGTAGCGCTa                                       >  1:594301/1‑40 (MQ=255)
                                                       ttATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGAc                                <  1:273188/42‑1 (MQ=255)
                                                       ttATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGACAGGGATGGCGCTCAAATTa             >  1:1201055/1‑61 (MQ=255)
                                                         aTCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGACAGGGATGGCGCTCAAATTATGCCCACAtct  >  1:1286967/1‑70 (MQ=255)
                                                         aTCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGACAGGGATGGCGCTCAAATTATGCCCACAtct  >  1:1101217/1‑70 (MQ=255)
                                                         aTCCCGCACGTAAATTTGGCGTACGTAGCGCTATGCCGACAGGGATg                         >  1:1356581/1‑47 (MQ=255)
                                                         aTCCCGCACGTAAATTTGGCGTACGTAGCGCTATGCCGACAGGGATg                         >  1:655815/1‑47 (MQ=255)
                                                            ccGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGACAGGGATGGCGCTCAAATTATGCCCACAtc   <  1:1646064/66‑1 (MQ=255)
                                                                |                                                              
TCCCTATTGCTATTGGCGGCAGCCGCGAACGTCGGGGGGTGATCAGCACCGCCAATTATCCCGCGCGTAAATTTGGCGTACGTAGCGCTATGCCGACAGGGATGGCGCTCAAATTATGCCCACATCT  >  W3110S.gb/250977‑251103

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: