Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,398,077 C→T 22.7% A53T (GCA→ACA)  yhdE ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,398,0770CT22.7% 40.3 / 8.9 22A53T (GCA→ACA) yhdEconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (7/10);  new base T (3/2);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGATAGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGC  >  W3110S.gb/3398016‑3398145
                                                             |                                                                    
aCGATAGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGTTTTCTcg                                                               >  1:1483828/1‑69 (MQ=255)
aCGATAGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGTTTTCTcg                                                               >  1:2180114/1‑69 (MQ=255)
aCGATAGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGTTTTATcg                                                               >  1:1412490/1‑69 (MQ=255)
      gTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCgcgcgc                                                          <  1:3004168/68‑1 (MQ=255)
         tCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGa                                                      >  1:1876868/1‑69 (MQ=255)
               cccAGCACCGGGAGATCCTTCGCCGTTTGGGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGcac                                                 >  1:1198140/1‑68 (MQ=255)
               cccAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGcac                                                 >  1:1828098/1‑68 (MQ=255)
                  aGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGTTTTCTCGCGCGCCAGACGCacaaca                                             <  1:2828088/69‑1 (MQ=255)
                  aGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCacaaca                                             <  1:1364237/69‑1 (MQ=255)
                  aGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCacaaca                                             <  1:648791/69‑1 (MQ=255)
                  aGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCacaaca                                             >  1:597064/1‑69 (MQ=255)
                  aGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCacaaca                                             <  1:319363/69‑1 (MQ=255)
                  aGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCacaaca                                             <  1:2900511/69‑1 (MQ=255)
                   gCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGTTTTCTCGCGCGCCAGACGCacaaca                                             <  1:2240630/68‑1 (MQ=255)
                    cACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACAt                                            >  1:2754073/1‑68 (MQ=255)
                            gaTCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCacaa                                               <  1:3085326/57‑1 (MQ=255)
                            gaTCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCacaa                                               <  1:1416034/57‑1 (MQ=255)
                            gaTCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCacaa                                               <  1:884622/57‑1 (MQ=255)
                               ccTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTgcgcgc                                <  1:1556554/69‑1 (MQ=255)
                                 ttCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTgcg                                   >  1:1676964/1‑64 (MQ=255)
                                 ttCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTgcg                                   >  1:1372369/1‑64 (MQ=255)
                                                         cGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGctcct        <  1:2343330/67‑1 (MQ=255)
                                                             cTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGctc          <  1:1388511/61‑1 (MQ=255)
                                                             cTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGc  >  1:2681995/1‑69 (MQ=255)
                                                             cTTTATCGCGCGCCAGACGCACAACATACTGCTGCGCGct                               >  1:3307192/1‑40 (MQ=255)
                                                             |                                                                    
ACGATAGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGC  >  W3110S.gb/3398016‑3398145

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: