Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 219,364 T→C P77P (CCA→CCG tsaA ← tRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913219,3640TC100.0% 93.3 / NA 27P77P (CCA→CCGtsaAtRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/14);  total (13/14)

CCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGG  >  NC_000913/219230‑219461
                                                                                                                                      |                                                                                                 
ccGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGAc                                                                                               >  2:775850/1‑139 (MQ=255)
ccGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGAc                                                                                               >  1:114820/1‑139 (MQ=255)
                  cTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAgc                                                                             >  1:791470/1‑139 (MQ=255)
                               aaaCAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGt                                                                >  1:19201/1‑139 (MQ=255)
                               aaaCAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGt                                                                <  1:312647/139‑1 (MQ=255)
                                 acaacCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCAt                                                                   <  2:653803/134‑1 (MQ=255)
                                 acaacCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCAt                                                                   >  1:653803/1‑134 (MQ=255)
                                               ctctACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGaaa                                                <  2:791470/139‑1 (MQ=255)
                                                 ctACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAg                                              <  2:25338/139‑1 (MQ=255)
                                                   aCCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGa                                            >  1:91859/1‑139 (MQ=255)
                                                              gCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGcac                                                                                           <  1:549182/81‑1 (MQ=255)
                                                              gCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGcac                                                                                           >  2:549182/1‑81 (MQ=255)
                                                                   ttGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGaa                            <  2:91859/139‑1 (MQ=255)
                                                                   ttGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGaa                            <  1:775850/139‑1 (MQ=255)
                                                                                      agagCGTGTTGCGAAAACCCCCATTCTAGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATAGCTGAACGCTTCCAGGCCGCGTACg         >  1:541634/1‑139 (MQ=255)
                                                                                         gCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCg      <  2:62016/139‑1 (MQ=255)
                                                                                            tgtTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGaa           <  1:314115/131‑1 (MQ=255)
                                                                                            tgtTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGaa           >  2:314115/1‑131 (MQ=255)
                                                                                             gtTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCgg  <  2:114820/139‑1 (MQ=255)
                                                                                               tGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGa                                                                                    <  1:586388/55‑1 (MQ=255)
                                                                                               tGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGa                                                                                    >  2:586388/1‑55 (MQ=255)
                                                                                                        ccccATTCTTGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGc                                                                                 >  2:702856/1‑49 (MQ=39)
                                                                                                        ccccATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGc                                                                                 <  1:702856/49‑1 (MQ=255)
                                                                                                                tGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGaa           <  2:307814/111‑1 (MQ=255)
                                                                                                                tGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGaa           >  1:307814/1‑111 (MQ=255)
                                                                                                                       aCCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAg                                              >  1:345156/1‑69 (MQ=255)
                                                                                                                       aCCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAg                                              <  2:345156/69‑1 (MQ=255)
                                                                                                                                   cggcggACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAg                                              <  1:644167/57‑1 (MQ=255)
                                                                                                                                      |                                                                                                 
CCGAGCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGG  >  NC_000913/219230‑219461

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: