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%	144.0 / NA	41	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 (19/22);  total (19/22)

GCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGG  >  NC_000913/219234‑219474                                                                                                                                   |                                                                                                               gCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGcac                                                                                                        <  2:198197/139‑1 (MQ=255)    tCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGt                                                                                                     <  2:957101/139‑1 (MQ=255)                ggTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAgcc                                                                                         <  2:837152/138‑1 (MQ=255)                  tCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAgccgc                                                                                       >  1:50171/1‑139 (MQ=255)                   cTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAgccgcc                                                                                      <  2:50171/139‑1 (MQ=255)                    tttATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCt                                                                                     >  1:830423/1‑139 (MQ=255)                            aaaCAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAgccg                                                                                        >  2:61903/1‑128 (MQ=255)                            aaaCAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAgccg                                                                                        <  1:61903/128‑1 (MQ=255)                            aaaCAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAg                                                                                           <  2:1006189/125‑1 (MQ=255)                            aaaCAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAg                                                                                           >  1:1006189/1‑125 (MQ=255)                            aaaCAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTcc                                                                                  <  1:732617/134‑1 (MQ=255)                            aaaCAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTcc                                                                                  >  2:732617/1‑134 (MQ=255)                                          aGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGa                                                               >  1:198300/1‑139 (MQ=255)                                             tctACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGaaaa                                                            <  2:198300/139‑1 (MQ=255)                                                      gACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATa                                                   <  1:1023142/139‑1 (MQ=255)                                                         aTGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGc                                                                                              <  1:382294/93‑1 (MQ=255)                                                         aTGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGc                                                                                              >  2:382294/1‑93 (MQ=255)                                                           gCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGa                                                               >  1:141087/1‑122 (MQ=255)                                                           gCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGa                                                               <  2:141087/122‑1 (MQ=255)                                                           gCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATgg                                              >  1:466498/1‑139 (MQ=255)                                                           gCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATgg                                              >  2:914117/1‑139 (MQ=255)                                                           gCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATgg                                              >  1:466027/1‑139 (MQ=255)                                                                ttGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGaa                                         <  1:914117/139‑1 (MQ=255)                                                                             gAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAg                                                           <  1:480208/108‑1 (MQ=255)                                                                             gAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAg                                                           >  2:480208/1‑108 (MQ=255)                                                                                       cGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGt                  <  2:830423/139‑1 (MQ=255)                                                                                           ttGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGc              >  1:386906/1‑139 (MQ=255)                                                                                            tGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGcc             <  2:466498/139‑1 (MQ=255)                                                                                            tGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGcc             <  2:386906/139‑1 (MQ=255)                                                                                                         aTTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGaaaa                                                            >  1:763880/1‑79 (MQ=255)                                                                                                         aTTCTGGCGTTACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGaaaa                                                            <  2:763880/79‑1 (MQ=255)                                                                                                                  ttACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGaa                                                                     <  1:892779/61‑1 (MQ=255)                                                                                                                  ttACCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGaa                                                                     >  2:892779/1‑61 (MQ=255)                                                                                                                    aCCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTa     >  2:407294/1‑123 (MQ=255)                                                                                                                    aCCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTa     <  1:407294/123‑1 (MQ=255)                                                                                                                    aCCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAggg  <  1:660435/126‑1 (MQ=255)                                                                                                                    aCCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAgg   <  2:395430/125‑1 (MQ=255)                                                                                                                    aCCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAgg   >  1:395430/1‑125 (MQ=255)                                                                                                                    aCCGCCGAGGCGCGGCGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGAACCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAggg  >  2:660435/1‑126 (MQ=255)                                                                                                                         cgAGGCGCGGCGGACGCACAGTCGGACGCCAgccgcc                                                                                      >  1:944280/1‑37 (MQ=255)                                                                                                                         cgAGGCGCGGCGGACGCACAGTCGGACGCCAgccgcc                                                                                      <  2:944280/37‑1 (MQ=255)                                                                                                                                   |                                                                                                               GCTTCAGAATCACGCTGTCTTTATGGCAAACAACCTCTTTCAGCTCTACCAGCGACATGCCAATTGGGTTAGGGCGGAAAGTAGAGCGTGTTGCGAAAACCCCCATTCTGGCGTTACCGCCGAGGCGCGGTGGACGCACAGTCGGACGCCAGCCGCCTTCCATCGTTTGATGAAAGACGAAAAGGATCCATAAATGGCTGAACGCTTCCAGGCCGCGAACGGCGTCGGCCTGGTTGTAGGG  >  NC_000913/219234‑219474

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