Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 218,942 A→G V218A (GTT→GCT)  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_000913218,9420AG93.5% 94.6 / ‑2.8 31V218A (GTT→GCT) tsaAtRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (17/12);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

GAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTC  >  NC_000913/218795‑219087
                                                                                                                                                   |                                                                                                                                                 
gAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGt                                                                                                                                                >  2:317419/1‑151 (MQ=255)
         aCAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCAt                                                                                                                                       >  2:196325/1‑151 (MQ=255)
                  aGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGcc                                                                                                                                 >  1:112615/1‑148 (MQ=255)
                     cGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACgg                                                                                                                           >  1:109929/1‑151 (MQ=255)
                                cAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTg                                                                                                                >  2:346332/1‑151 (MQ=255)
                                       aaGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGttt                                                                                                         >  1:600164/1‑151 (MQ=255)
                                       aaGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGttt                                                                                                         >  1:93454/1‑151 (MQ=255)
                                              ttAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCAcc                                                                                                  >  2:518605/1‑151 (MQ=255)
                                                     aCCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACg                                                                                            <  2:459301/150‑1 (MQ=255)
                                                     aCCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGa                                                                                           <  1:122804/151‑1 (MQ=255)
                                                      ccAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGAt                                                                                          >  2:448345/1‑151 (MQ=255)
                                                             aaaaGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCgg                                                                                   >  1:432876/1‑151 (MQ=255)
                                                                aGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGgcg                                                                                <  2:223515/151‑1 (MQ=255)
                                                                    aTGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCggg                                                                            <  1:544021/151‑1 (MQ=255)
                                                                     tGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGt                                                                           >  2:471197/1‑151 (MQ=255)
                                                                        cAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCGCCTTTACGATAGGCCGGGCGCGGGTcc                                                                         <  2:431523/150‑1 (MQ=255)
                                                                         aaaaGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTcc                                                                         >  2:471953/1‑149 (MQ=255)
                                                                           aaGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTgcg                                                                     <  1:363489/151‑1 (MQ=255)
                                                                                aaGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGt                                                                >  1:152363/1‑151 (MQ=255)
                                                                                     cGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTc                                                           <  1:613253/151‑1 (MQ=255)
                                                                                       tGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTcgc                                                         >  1:174303/1‑151 (MQ=255)
                                                                                        gAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTcgc                                                         <  1:158897/150‑1 (MQ=255)
                                                                                          aTTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCgg                                                       >  1:370540/1‑150 (MQ=255)
                                                                                                   ggTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCg                                             >  1:441536/1‑151 (MQ=255)
                                                                                                         aGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGGCGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACt                                       <  1:326634/151‑1 (MQ=255)
                                                                                                                  aCTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCg                                    >  2:430116/1‑145 (MQ=255)
                                                                                                                  aCTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCg                                    <  1:430123/145‑1 (MQ=255)
                                                                                                                  aCTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAAc                              <  2:246581/151‑1 (MQ=255)
                                                                                                                                  cGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCa               <  1:226376/150‑1 (MQ=255)
                                                                                                                                       aCGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAgct         <  2:95017/151‑1 (MQ=255)
                                                                                                                                            ccAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGCTTTGACCGTTTCCTCACCTCTACGATAGGCCGGGCGCGGGTCCTGCGACAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCttt    <  1:539420/151‑1 (MQ=255)
                                                                                                                                              agcgagcgTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTc  >  1:447308/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAAGGTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTC  >  NC_000913/218795‑219087

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