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,9420AG98.2% 187.0 / ‑6.7 54V218A (GTT→GCT) tsaAtRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (31/22);  minor base T (1/0);  total (32/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 = 1.00e+00

GTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCG  >  NC_000913/218799‑219088
                                                                                                                                               |                                                                                                                                                  
gTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTaaa                                                                                                                                             <  2:610970/151‑1 (MQ=255)
gTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTaaa                                                                                                                                             <  2:60631/151‑1 (MQ=255)
          cTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGc                                                                                                                                   >  2:234373/1‑151 (MQ=255)
             aaGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAAGCATGTTGCCa                                                                                                                                 <  2:626527/150‑1 (MQ=255)
                  gAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGc                                                                                                                           >  1:272758/1‑151 (MQ=255)
                       tGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAgg                                                                                                                      >  1:534754/1‑151 (MQ=255)
                           aCAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGtttt                                                                                                                  <  1:383928/151‑1 (MQ=255)
                            cAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTg                                                                                                                 <  2:307863/151‑1 (MQ=255)
                               aaaaaaGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGccc                                                                                                              >  1:416357/1‑151 (MQ=255)
                               aaaaaaGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGccc                                                                                                              >  1:537277/1‑151 (MQ=255)
                                aaaaaGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCg                                                                                                             >  2:611982/1‑151 (MQ=255)
                                 aaaaGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGt                                                                                                            <  1:70057/151‑1 (MQ=255)
                                 aaaaGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGt                                                                                                            >  2:408430/1‑151 (MQ=255)
                                 aaaaGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGt                                                                                                            >  2:482372/1‑151 (MQ=255)
                                                    aGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAg                                                                                         <  1:662108/151‑1 (MQ=255)
                                                     gTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAgg                                                                                        >  1:529173/1‑151 (MQ=255)
                                                       gCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGcc                                                                                      >  1:643665/1‑151 (MQ=255)
                                                       gCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGc                                                                                       <  2:603440/150‑1 (MQ=255)
                                                        cAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCg                                                                                     >  2:593738/1‑151 (MQ=255)
                                                         aaaaGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCgg                                                                                    >  2:621460/1‑151 (MQ=255)
                                                          aaaGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCggg                                                                                   >  1:515132/1‑151 (MQ=255)
                                                                   tCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTcc                                                                          >  2:613064/1‑151 (MQ=255)
                                                                   tCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTcc                                                                          >  2:626750/1‑151 (MQ=255)
                                                                   tCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTcc                                                                          >  1:42073/1‑151 (MQ=255)
                                                                     aaaaGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTg                                                                        <  1:112824/151‑1 (MQ=255)
                                                                      aaaGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTgc                                                                       >  2:354154/1‑151 (MQ=255)
                                                                      aaaGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTgc                                                                       >  2:146587/1‑151 (MQ=255)
                                                                         gagaAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGcc                                                                    <  2:418513/151‑1 (MQ=255)
                                                                         gagaAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGcc                                                                    >  2:39767/1‑151 (MQ=255)
                                                                             agtgCGAGAATTTAACGCGGCTCCCGCGCAAAGACTTCAAACCCTGCGTCGGTGACGCGTAAGCGATCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATATCCCGGGCGCGGGTCCTgctctggt                                                                 >  1:596198/4‑144 (MQ=255)
                                                                                 cGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTc                                                            >  1:592618/1‑151 (MQ=255)
                                                                                 cGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTc                                                            >  2:394837/1‑151 (MQ=255)
                                                                                    gAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTcgcg                                                         <  1:377442/151‑1 (MQ=255)
                                                                                        ttAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGAt                                                      <  2:151502/150‑1 (MQ=255)
                                                                                        ttAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGAt                                                      <  1:572174/150‑1 (MQ=255)
                                                                                          aaCGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATaa                                                    <  1:291551/150‑1 (MQ=255)
                                                                                            cgcgGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACa                                                 <  2:548137/151‑1 (MQ=255)
                                                                                             gcgGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAg                                                <  1:621530/151‑1 (MQ=255)
                                                                                             gcgGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAg                                                <  2:579793/151‑1 (MQ=255)
                                                                                                  tCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTc                                                            >  2:483986/1‑134 (MQ=255)
                                                                                                  tCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTc                                                            <  1:483984/134‑1 (MQ=255)
                                                                                                     aGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACt                                        >  2:25431/1‑151 (MQ=255)
                                                                                                        gcAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCg                                     >  2:471285/1‑151 (MQ=255)
                                                                                                           aaGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGAt                                  >  1:186534/1‑151 (MQ=255)
                                                                                                            aGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGTCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGTTAAACAGCGTTAACTGCGGAt                                  >  2:242046/1‑150 (MQ=255)
                                                                                                                 tCAAAACCTGCGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCACGGATAAACAGCGTTAACTGCGGATAACGc                             <  2:124335/150‑1 (MQ=255)
                                                                                                                           cGTCGGTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTACCTGCGGATAACGCTTCTCCAGc                    <  1:116913/149‑1 (MQ=255)
                                                                                                                               ggTGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCaaa              <  1:611971/151‑1 (MQ=255)
                                                                                                                                 tGACGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCaaa              >  1:447818/1‑149 (MQ=255)
                                                                                                                                   aCGCGCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAgct          >  2:20012/1‑151 (MQ=255)
                                                                                                                                    cgcgCCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAgctg         <  2:633472/151‑1 (MQ=255)
                                                                                                                                     gcgcCAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAgctgc        >  2:47184/1‑151 (MQ=255)
                                                                                                                                         cAGCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCtttt    >  1:670402/1‑151 (MQ=255)
                                                                                                                                           gCGAGCGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCg  >  2:540153/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GTAGAACAGGCTCAAGGCGAGCCTGGGACAAAAAAAAGTGATTTAGTTTACCAGTGCAAAAGAAATGTCAAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCG  >  NC_000913/218799‑219088

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