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_000913	219,364	0	T	C	100.0%	117.0 / NA	33	P77P (CCA→CCG)	tsaA	tRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/17);  total (16/17)

ACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTA  >  NC_000913/219245‑219471                                                                                                                        |                                                                                                            aCGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGTGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGCGGCGCGGCGGACGCACAGTCGGACGcc                                                                                          >  2:108058/1‑139 (MQ=255)  cGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCa                                                                                         >  1:794086/1‑139 (MQ=255)                 aaaCAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGt                                                                          >  2:371490/1‑139 (MQ=255)                   acaacCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGttt                                                                        >  1:274036/1‑139 (MQ=255)                   acaacCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGttt                                                                        >  1:531262/1‑139 (MQ=255)                        ctctTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTtgatga                                                                   <  2:531262/139‑1 (MQ=255)                                 ctctACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGaaa                                                          >  2:994873/1‑139 (MQ=255)                                  tctACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTtgatga                                                                   >  2:1040578/1‑129 (MQ=255)                                  tctACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTtgatga                                                                   <  1:1040578/129‑1 (MQ=255)                                   ctACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAg                                                        <  2:347606/139‑1 (MQ=255)                                   ctACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAg                                                        <  1:994873/139‑1 (MQ=255)                                    tACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCg                                                                                                >  1:50104/1‑98 (MQ=255)                                    tACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCg                                                                                                <  2:50104/98‑1 (MQ=255)                                    tACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGa                                                                                              >  1:108458/1‑100 (MQ=255)                                    tACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGa                                                                                              <  2:108458/100‑1 (MQ=255)                                     aCCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGa                                                      <  1:853325/139‑1 (MQ=255)                                             cATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATaaa                                              <  2:794086/139‑1 (MQ=255)                                              aTGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGa                                                               >  2:40353/1‑121 (MQ=255)                                              aTGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGa                                                               <  1:40353/121‑1 (MQ=255)                                                     ttGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGaa                                      <  2:216231/139‑1 (MQ=255)                                                     ttGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGaa                                      <  1:465053/139‑1 (MQ=255)                                                                  gAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAgc                                                                                       >  1:838176/1‑77 (MQ=255)                                                                  gAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAgc                                                                                       <  2:838176/77‑1 (MQ=255)                                                                   aaaGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCcgc                        >  2:422856/1‑139 (MQ=255)                                                                              tgtTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGaa                     >  2:107203/1‑131 (MQ=255)                                                                              tgtTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGaa                     <  1:107203/131‑1 (MQ=255)                                                                                 tGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGcc          <  2:274036/139‑1 (MQ=255)                                                                                 tGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGcc          >  1:147893/1‑139 (MQ=255)                                                                                   cGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGcc          >  1:678835/1‑137 (MQ=255)                                                                                   cGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGcc          <  2:678835/137‑1 (MQ=255)                                                                                         cccccATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTa  <  2:334549/139‑1 (MQ=255)                                                                                                         aCCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAg                                                        <  2:553097/69‑1 (MQ=255)                                                                                                         aCCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAg                                                        >  1:553097/1‑69 (MQ=255)                                                                                                                        |                                                                                                            ACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTA  >  NC_000913/219245‑219471

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