Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,100,981 A→G 15.6% 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,9810AG15.6% 75.8 / 8.1 32H176H (CAT→CACyggHtRNA (m7G46) methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (5/22);  new base G (3/2);  total (8/24)
Fisher's exact test for biased strand distribution p-value = 8.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

TTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTT  >  W3110S.gb/3100912‑3101051
                                                                     |                                                                      
tttATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCata                                                                        <  1:2530701/70‑1 (MQ=255)
    taACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                                   <  1:592026/71‑1 (MQ=255)
    taACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                                   <  1:2762769/71‑1 (MQ=255)
    taACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                                   <  1:2759548/71‑1 (MQ=255)
    taACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                                   <  1:1786157/71‑1 (MQ=255)
    taACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                                   <  1:1932186/71‑1 (MQ=255)
    taACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                                   <  1:2704041/71‑1 (MQ=255)
     aaCCGTCAATAGAAGACATCACTTCAAGCATATGTT‑‑GCATAAGGTTCCCAGTCGGTCGCCATATGGAATAc                                                                >  1:294549/1‑71 (MQ=255)
                gaagaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGAATa                                                                 >  1:1239163/1‑61 (MQ=255)
                           ttCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                                   <  1:314860/48‑1 (MQ=255)
                           ttCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                                   <  1:2209898/48‑1 (MQ=255)
                                 cATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGaa                                                                   >  1:897718/1‑42 (MQ=255)
                                 cATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGaa                                                                   >  1:867121/1‑42 (MQ=255)
                                   tatGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGaa                                                                   <  1:332627/40‑1 (MQ=255)
                                   tatGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                                     <  1:1325549/38‑1 (MQ=255)
                                   tatGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATa                                                                 >  1:1471554/1‑42 (MQ=255)
                                     tGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAg                                  <  1:428719/71‑1 (MQ=255)
                                     tGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAg                                  >  1:1478181/1‑71 (MQ=255)
                                         ccGCATAAGGTTCCCAGTCGGTCGCCATATGGATTACGCCCCCCAGCTGCAGTTTGCt                                           <  1:2284212/58‑1 (MQ=255)
                                            cATAAGGTTCCCAGTCGGTCGCCATGTGGAATACGCCCCCCAGCTGCAGTTTGCtt                                          <  1:489246/56‑1 (MQ=255)
                                            cATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCa                           <  1:2716958/71‑1 (MQ=255)
                                            cATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCa                           <  1:2608298/71‑1 (MQ=255)
                                              tAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCaa                          <  1:2118744/70‑1 (MQ=255)
                                                 ggTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACgg                      <  1:1448004/71‑1 (MQ=255)
                                                       cAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCAcc                  >  1:1768981/1‑69 (MQ=255)
                                                          tcggtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCt                                           <  1:2449313/41‑1 (MQ=255)
                                                            ggtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGa           <  1:1774140/71‑1 (MQ=255)
                                                            ggtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGa           >  1:2187962/1‑71 (MQ=255)
                                                             gtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTg                <  1:623323/65‑1 (MQ=255)
                                                             gtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAAc             <  1:2630544/68‑1 (MQ=255)
                                                             gtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGa           <  1:2781716/70‑1 (MQ=255)
                                                             gtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGa           <  1:2072235/70‑1 (MQ=255)
                                                                gCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATAc        <  1:1352824/70‑1 (MQ=255)
                                                                    tatGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTg                <  1:2186693/58‑1 (MQ=255)
                                                                     atGGAATACGCCCCGCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGtt  >  1:405628/1‑71 (MQ=255)
                                                                     |                                                                      
TTTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGATACGGCGTT  >  W3110S.gb/3100912‑3101051

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: