Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,100,981 A→G 15.1% 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.1% 130.5 / 15.7 53H176H (CAT→CACyggHtRNA (m7G46) methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (15/30);  new base G (5/3);  total (20/33)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGA  >  W3110S.gb/3100913‑3101042
                                                                    |                                                             
ttATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCatat                                                              <  1:2372139/70‑1 (MQ=255)
 tataACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATg                                                             <  1:3586919/70‑1 (MQ=255)
 tataACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATg                                                             <  1:1655768/70‑1 (MQ=255)
 tataACCGTCAATAGAAGACATCACTTAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                            <  1:4438110/70‑1 (MQ=255)
  ataACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTgg                                                            <  1:4465611/70‑1 (MQ=255)
  ataACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                            <  1:4113034/70‑1 (MQ=255)
  ataACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                            <  1:1194481/70‑1 (MQ=255)
  ataACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                            <  1:4130071/70‑1 (MQ=255)
   taACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                          <  1:3868592/71‑1 (MQ=255)
        gTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                            <  1:1526442/64‑1 (MQ=255)
         tCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                            <  1:1633505/63‑1 (MQ=255)
         tCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATg                                                             <  1:3822450/62‑1 (MQ=255)
         tCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGa                                                           >  1:2784556/1‑64 (MQ=255)
           aaTAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGcccc                                                  >  1:896594/1‑71 (MQ=255)
           aaTAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGcccc                                                  >  1:196106/1‑71 (MQ=255)
               gaagaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATa                                                        >  1:3459576/1‑61 (MQ=255)
               gaagaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCa                                               >  1:1957175/1‑70 (MQ=255)
                 agaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAg                                              >  1:2209907/1‑69 (MQ=255)
                  gaCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTg                                           <  1:1028174/71‑1 (MQ=255)
                   aCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                            <  1:4434772/53‑1 (MQ=255)
                   aCATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                            <  1:1942163/53‑1 (MQ=255)
                    cATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATgg                                                            >  1:2134012/1‑52 (MQ=255)
                     aTCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCGAGCTGCAg                                        <  1:1295796/71‑1 (MQ=255)
                     aTCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAg                                        <  1:3823489/71‑1 (MQ=255)
                          ttCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGAATa                                                        >  1:3430020/1‑50 (MQ=255)
                          ttCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGAATa                                                        >  1:370726/1‑50 (MQ=255)
                          ttCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGAATACGc                                                     >  1:56252/1‑53 (MQ=255)
                           tCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCt                                  >  1:2823013/1‑71 (MQ=255)
                             aaGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATGTGGAATACGCCCCCCAGCTGCAGTTTGCtt                                 <  1:1610054/70‑1 (MQ=255)
                              aGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGGtttt                               <  1:4494626/71‑1 (MQ=255)
                              aGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCtttt                               <  1:3444027/71‑1 (MQ=255)
                              aGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCtttt                               <  1:2201667/71‑1 (MQ=255)
                               gCATATGTTCCGCATAAGGTTCCCAGTCGGGCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCttttt                              <  1:4103853/71‑1 (MQ=255)
                                  tatGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATa                                                        <  1:1252334/42‑1 (MQ=255)
                                  tatGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCttttt                              <  1:2443504/68‑1 (MQ=255)
                                  tatGTTCCGCATAAGGTTCCCAGTCGGGCGCCATATGGaa                                                          <  1:630372/40‑1 (MQ=255)
                                    tGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATTGAATACGCCCCCCAGCTGCAGTTTGc                                   >  1:978578/1‑61 (MQ=255)
                                    tGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                          >  1:3130218/1‑38 (MQ=255)
                                    tGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGaa                                                          >  1:2432237/1‑38 (MQ=255)
                                        ccGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTc                      >  1:3340186/1‑70 (MQ=255)
                                        ccGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCg                     >  1:1026731/1‑71 (MQ=255)
                                        ccGCATAAGCTTCCCAGTCGGTCGCCATGTGGAATa                                                        >  1:4524356/1‑36 (MQ=37)
                                           cATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCa                  <  1:1030012/71‑1 (MQ=255)
                                                     ccAGTCCGGCGCCATGTGGAATACGCCCCCCAGCTGCAGTTTg                                    <  1:4483750/43‑1 (MQ=37)
                                                       aGTCGGTCGCCATGTTGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTg       >  1:3003371/1‑70 (MQ=255)
                                                       aGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGc                   <  1:4381890/58‑1 (MQ=255)
                                                        gtcggtcgCCATATGGAATACGCCCCACAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGaa     >  1:672699/1‑71 (MQ=255)
                                                         tcggtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCtt                                 <  1:721124/42‑1 (MQ=255)
                                                           ggtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGa  <  1:3672944/71‑1 (MQ=255)
                                                           ggtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGa  >  1:1057507/1‑71 (MQ=255)
                                                           ggtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGa  <  1:449974/71‑1 (MQ=255)
                                                           ggtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGa  <  1:3329166/71‑1 (MQ=255)
                                                           ggtcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGa  <  1:3925331/71‑1 (MQ=255)
                                                             tcgCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGa  <  1:1024117/69‑1 (MQ=255)
                                                                  atatGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCa                  <  1:4533274/48‑1 (MQ=255)
                                                                    |                                                             
TTATAACCGTCAATAGAAGACATCACTTCAAGCATATGTTCCGCATAAGGTTCCCAGTCGGTCGCCATATGGAATACGCCCCCCAGCTGCAGTTTGCTTTTTACCAGTTCGGCAAACGGCACCTGAACGA  >  W3110S.gb/3100913‑3101042

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: