Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,100,981 A→G 27.8% H176H (CAT→CAC yggH ← tRNA (m7G46) methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,100,9810AG27.8% 51.2 / 24.4 36H176H (CAT→CACyggHtRNA (m7G46) methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (13/13);  new base G (7/3);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGC  >  W3110S.gb/3100911‑3101048
                                                                      |                                                                   
ttttATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCata                                                                     >  1:929591/1‑71 (MQ=255)
ttttATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCata                                                                     >  1:3882947/1‑71 (MQ=255)
ttttATAACCGTCAATAGAAGACAACACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCata                                                                     >  1:2373503/1‑71 (MQ=255)
    ataACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTgg                                                                  <  1:2733294/70‑1 (MQ=255)
    ataACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                                  <  1:1992487/70‑1 (MQ=255)
     taACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCata                                                                     >  1:1618086/1‑66 (MQ=255)
     taACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                                  <  1:1429355/69‑1 (MQ=255)
     taACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                                <  1:1618134/71‑1 (MQ=255)
     taACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                                <  1:1660282/71‑1 (MQ=255)
           tCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCata                                                                     <  1:2084566/60‑1 (MQ=255)
           tCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCata                                                                     <  1:681814/60‑1 (MQ=255)
             aaTAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGccc                                                         >  1:3196197/1‑70 (MQ=255)
             aaTAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGccc                                                         >  1:304128/1‑70 (MQ=255)
             aaTAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGccc                                                         >  1:2969672/1‑70 (MQ=255)
                 gaagaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGaa                                                                >  1:2885066/1‑59 (MQ=255)
                 gaagaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGaa                                                                >  1:2736054/1‑59 (MQ=255)
                    gaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                                <  1:3352342/56‑1 (MQ=255)
                        tCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCa                                               <  1:4058236/69‑1 (MQ=255)
                        tCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCACCCCAGCTGCa                                               <  1:3990705/69‑1 (MQ=255)
                          aCTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATg                                                                   >  1:996124/1‑47 (MQ=255)
                          aCTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGttt                                           >  1:2538186/1‑71 (MQ=255)
                          aCTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGttt                                           >  1:736977/1‑71 (MQ=255)
                              cAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGaa                                                                >  1:1211595/1‑46 (MQ=255)
                              cAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGaa                                                                >  1:1685457/1‑46 (MQ=255)
                              cAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGaa                                                                >  1:2553721/1‑46 (MQ=255)
                              cAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGaa                                                                >  1:2469189/1‑46 (MQ=255)
                              cAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGc                                         <  1:28055/69‑1 (MQ=255)
                              cAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCtt                                       >  1:209440/1‑71 (MQ=255)
                                aGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCtttt                                     <  1:856415/71‑1 (MQ=255)
                                aGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCtttt                                     <  1:71109/71‑1 (MQ=255)
                                    tatGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                                  <  1:6262/38‑1 (MQ=255)
                                      tGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGAATa                                                              >  1:1061368/1‑40 (MQ=255)
                                      tGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAg                               >  1:1361766/1‑71 (MQ=255)
                                      tGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAg                               >  1:1356978/1‑71 (MQ=255)
                                          ccGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCg                           >  1:3438567/1‑71 (MQ=255)
                                                 aGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAAc                     <  1:3520073/70‑1 (MQ=255)
                                                   gTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGc                  <  1:1052707/71‑1 (MQ=255)
                                                        cAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCttttt                                    >  1:3384085/1‑48 (MQ=255)
                                                             ggtcgCCATGTGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGa        <  1:237818/71‑1 (MQ=255)
                                                             ggtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAg                               >  1:569585/1‑48 (MQ=255)
                                                                 gCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACg    >  1:1489647/1‑71 (MQ=255)
                                                                   cATGTGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGc  <  1:1875705/71‑1 (MQ=255)
                                                                      |                                                                   
TTTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGC  >  W3110S.gb/3100911‑3101048

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: