Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 219,050 T→C 9.1% Q182R (CAG→CGG)  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,0500TC 9.1% 134.8 / 2.8 44Q182R (CAG→CGG) tsaAtRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (20/20);  new base C (2/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCGGCAAAGGGGAGATACGGTTTGATATCCAC  >  NC_000913/218904‑219198
                                                                                                                                                  |                                                                                                                                                    
cAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGc                                                                                                                                                    <  1:554788/149‑1 (MQ=255)
cAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGc                                                                                                                                                    <  2:622851/149‑1 (MQ=255)
 aaaGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCg                                                                                                                                                   <  4:491531/149‑1 (MQ=255)
       ttCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATGATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAAc                                                                                                                                             >  3:450906/1‑149 (MQ=255)
              ccTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTtctc                                                                                                                                      <  3:464879/149‑1 (MQ=255)
                                cAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCtt                                                                                                                    <  3:416944/149‑1 (MQ=255)
                                   cGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACg                                                                                                                                            >  4:199632/1‑122 (MQ=255)
                                   cGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACg                                                                                                                                            <  3:199632/122‑1 (MQ=255)
                                        gTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACct                                                                                                            >  1:35979/1‑149 (MQ=255)
                                                 aTGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTaa                                                                                                   >  2:69522/1‑149 (MQ=255)
                                                   gTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGACAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTaaaa                                                                                                 <  4:413880/149‑1 (MQ=255)
                                                    tAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAAc                                                                                                <  1:130755/149‑1 (MQ=255)
                                                                           ccGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGgc                                                                         >  3:365750/1‑149 (MQ=255)
                                                                                  cTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACttt                                                                  <  4:365750/149‑1 (MQ=255)
                                                                                   tCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCAt                                                                                      <  2:263065/128‑1 (MQ=255)
                                                                                   tCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCAt                                                                                      >  1:263065/1‑128 (MQ=255)
                                                                                   tCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTg                                                                           <  4:235311/139‑1 (MQ=255)
                                                                                   tCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTg                                                                           >  3:235311/1‑139 (MQ=255)
                                                                                   tCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTg                                                                 >  2:373035/1‑149 (MQ=255)
                                                                                   tCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTg                                                                 >  1:269060/1‑149 (MQ=255)
                                                                                     acatttACGATAGGCCGGACGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTgcg                                                               >  4:228307/4‑149 (MQ=255)
                                                                                       cTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCa                                                             >  3:116347/1‑149 (MQ=255)
                                                                                                     gggCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGctggcactggc                                               >  4:43040/1‑149 (MQ=255)
                                                                                                       gcgcGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCAt                                             >  1:1455/1‑149 (MQ=255)
                                                                                                             gTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGaa                                       <  4:450906/149‑1 (MQ=255)
                                                                                                               ccTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGc                                     >  3:216716/1‑149 (MQ=255)
                                                                                                               ccTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGc                                     >  4:246255/1‑149 (MQ=255)
                                                                                                                  gcgcCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGAt                                  >  1:209972/1‑149 (MQ=255)
                                                                                                                   cgcCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGAtt                                 <  1:373035/149‑1 (MQ=255)
                                                                                                                    gcCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTc                                <  4:116347/149‑1 (MQ=255)
                                                                                                                     ccAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCg                               <  3:228307/149‑1 (MQ=255)
                                                                                                                      cAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCgg                              >  3:280891/1‑149 (MQ=255)
                                                                                                                       aGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCGGc                             <  3:246255/149‑1 (MQ=255)
                                                                                                                           cTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCGGCAAAg                         <  1:151139/149‑1 (MQ=255)
                                                                                                                            ttCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGaa                                       <  1:413046/134‑1 (MQ=255)
                                                                                                                            ttCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGaa                                       >  2:413046/1‑134 (MQ=255)
                                                                                                                            ttCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCGGCAAAgg                        >  1:90124/1‑149 (MQ=255)
                                                                                                                                cgGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCGGCAAAGGGgag                    <  1:69522/149‑1 (MQ=255)
                                                                                                                                   tggagaaGCGTTATCCGCAGTTAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCgg                              <  1:119337/130‑1 (MQ=255)
                                                                                                                                   tggagaaGCGTTATCCGCAGTTAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCgg                              >  2:119337/7‑136 (MQ=255)
                                                                                                                                      aaCAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCGGCAAAGGGGAGATACgg              <  4:216716/149‑1 (MQ=255)
                                                                                                                                          gCGTTATCCGCAGTTAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGAtt                                 >  3:464997/1‑126 (MQ=255)
                                                                                                                                          gCGTTATCCGCAGTTAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGAtt                                 <  4:464997/126‑1 (MQ=255)
                                                                                                                                             ttAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCGGCAAAGGGGAGATACGGTTTGata       >  4:420825/1‑149 (MQ=255)
                                                                                                                                                  tGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCGGCAAAGGGGAGATACGGTTTGATATCCAc  >  4:482350/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCACTGGCATCGGGAAGCGATTCGGCAAAGGGGAGATACGGTTTGATATCCAC  >  NC_000913/218904‑219198

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