Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,225,448 G→A 57.9% P59P (CCG→CCA ycgI → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,225,4480GA57.9% 1.7 / 20.6 19P59P (CCG→CCAycgIhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base A (3/8);  total (6/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAG  >  W3110S.gb/1225386‑1225515
                                                              |                                                                   
agcagTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAAcgc                                                             <  1:451879/71‑1 (MQ=255)
agcagTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCATGGAAcgc                                                             <  1:741729/71‑1 (MQ=255)
agcagTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCATGGAAcgc                                                             <  1:1853449/71‑1 (MQ=255)
    gTACGTTTAGCGGACATCTCGTGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCATgg                                                                  <  1:1551194/62‑1 (MQ=255)
    gTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCATgg                                                                  <  1:732128/62‑1 (MQ=255)
         tttAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTgg                                                    <  1:2480878/71‑1 (MQ=255)
            aGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAAcgc                                                             <  1:451878/59‑1 (MQ=255)
                          ggtggggtATATCCATGGTGCCGGTGTGGAATCCCCATGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCt                                   <  1:739780/71‑1 (MQ=255)
                              gtgtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGTTCGTTCTGAc                                >  1:1571500/3‑70 (MQ=255)
                                   taGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTg                                  <  1:2356273/63‑1 (MQ=40)
                                     gCCATGGTGCCGGTGTGGAATCCCCATGGAACGCGGTGGCTGGTGTGaa                                              >  1:1538253/1‑49 (MQ=255)
                                         tGGTGCCGGTGTGGAATCCCCATGGAACGCGGTGGCTGGTTTGAACTGGTCGTTCTGACCATCAACGAt                      <  1:1951001/69‑1 (MQ=14)
                                         tGGTGCCGGTGTGGAATCCCCATGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGAt                      <  1:1702171/69‑1 (MQ=255)
                                         tGGTGCCGGTCTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTa                    <  1:2282493/71‑1 (MQ=255)
                                                    tgGAATCCCCATGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGAtt                     >  1:988520/1‑59 (MQ=255)
                                                    tgGAATCCCCATGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTa                    <  1:1889332/60‑1 (MQ=255)
                                                      gAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAAc                 >  1:2516219/1‑61 (MQ=255)
                                                         tCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGa                                 >  1:750219/1‑42 (MQ=38)
                                                           cccATGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAg  >  1:954992/1‑71 (MQ=255)
                                                              |                                                                   
AGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAG  >  W3110S.gb/1225386‑1225515

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: