Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,036,407 A→G 64.3% E72E (GAA→GAG)  hyaD → protein involved in processing of HyaA and HyaB proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,036,4070AG64.3% 9.0 / 9.6 14E72E (GAA→GAG) hyaDprotein involved in processing of HyaA and HyaB proteins
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base G (4/5);  total (7/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

TATGTCGAAAGCGCCAGCCATCTGTTGATTCTCGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCG  >  W3110S.gb/1036351‑1036461
                                                        |                                                      
tATGTCGAAAGCGCCAGCCATCTGTTGATTCTCGATGCCATTGACTACGGGctggaacc                                                      <  1:1105128/59‑1 (MQ=255)
tATGTCGAAAGCGCCAGCCATCTGTTGATTCTCGATGCCATTGACTACGGGctggaacc                                                      <  1:1504675/59‑1 (MQ=255)
  tGTCGAAAGCGCCAGCCATCTGTTGATTCTCGATGCCATTGACTACGGGCTGGAGCCTGGAACGCTGCGa                                         <  1:415457/70‑1 (MQ=255)
    tCGAAAGCGCCAGCCATCTGTTGATTCTCGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAAcc                                      <  1:619498/71‑1 (MQ=255)
       aaaGCGCCAGCCATCTGTTGATTCTCGATGCCATTGACTACGGGCTGGAGCCTGGAACGc                                              <  1:193160/60‑1 (MQ=255)
           cgcCAGCCATCTGTTGATTCTCGATGCCATTGACTACGGGCTGGAGCCTGGAACGCTGCGAACCTATGCCg                               >  1:26589/1‑71 (MQ=255)
             ccagccaTCTGTTGATTCTCGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGga                             <  1:1769478/71‑1 (MQ=255)
                 ccaTCTGTTGATTCTCGATGCCATTGACTACGGGCTGGAGCCTGGAACGCTGCGAACCTAt                                   <  1:1424155/61‑1 (MQ=255)
                 ccaTCTGTTGATTCTCGATGCCATTGACTACGGGCTGGAGCCTGGAACGCTGCGAACCTAt                                   <  1:1240891/61‑1 (MQ=255)
                         tGATTCTCGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCg                 >  1:1592115/1‑71 (MQ=255)
                          gATTCTCGATGCCATTGACTACGGGCTGGAGCCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCgg                >  1:765320/1‑71 (MQ=255)
                            ttCTCGATGCCATTGACTACGGGctggaacctggaa                                                 >  1:702227/1‑36 (MQ=255)
                              ctcGATGCCATTGACTACGGGCTGGAGCCTGGAACGCTGCGAACCTATGCCGGAGaa                          >  1:1018769/1‑57 (MQ=255)
                              ctcGATGCCATTGACTACGGGCTGGAGCCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCtt             <  1:829051/70‑1 (MQ=255)
                               tcGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCg                                          >  1:429482/1‑40 (MQ=255)
                               tcGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTc                   >  1:1291093/1‑63 (MQ=255)
                                          gACTACGGGCTGGAGCCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGcgcg  >  1:626088/1‑69 (MQ=255)
                                                        |                                                      
TATGTCGAAAGCGCCAGCCATCTGTTGATTCTCGATGCCATTGACTACGGGCTGGAACCTGGAACGCTGCGAACCTATGCCGGAGAACGCATTCCGGCTTATCTCAGCGCG  >  W3110S.gb/1036351‑1036461

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: ‑