Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,781,563	G→A	E117K (GAA→AAA)	trmL →	tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,781,563	0	G	A	100.0%	136.9 / NA	35	E117K (GAA→AAA)	trmL	tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/18);  total (17/18)

GAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTGTAGGCCGGATAGATG  >  NC_000913/3781416‑3781709                                                                                                                                                    |                                                                                                                                                   gAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAAc                                                                                                                                                 <  2:150997/151‑1 (MQ=255)                cAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCAtt                                                                                                                                   >  1:267081/1‑150 (MQ=255)                      cTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGAt                                                                                                                              <  1:163211/149‑1 (MQ=255)                        gTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGt                                                                                                                           <  1:52476/150‑1 (MQ=255)                          tCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGc                                                                                                                         <  2:105061/150‑1 (MQ=255)                                       gAAAGGTACACCTGCTCACAGCGCCGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCgcagc                                                                                                           <  2:34489/151‑1 (MQ=255)                                         aaGGTACACCTGCTCACAGCGCCGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCAt                                                                                                         <  2:505428/151‑1 (MQ=255)                                            gTACACCTGCTCACAGCGCCGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGa                                                                                                       >  2:572803/1‑150 (MQ=255)                                             tACACCTGCTCACAGCGCCGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAAt                                                                                                     >  2:663842/1‑151 (MQ=255)                                                         caGCGCCGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCg                                                                                         <  2:122575/151‑1 (MQ=255)                                                              ccGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTc                                                                                    <  1:17414/151‑1 (MQ=255)                                                                             aTGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATg                                                                      >  1:436782/1‑150 (MQ=255)                                                                              tGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGa                                                                     >  2:127269/1‑150 (MQ=255)                                                                               ggCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAg                                                                   <  1:295863/151‑1 (MQ=255)                                                                                gCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGc                                                                  <  2:275311/151‑1 (MQ=255)                                                                                           aTGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGt                                                       <  2:596984/151‑1 (MQ=255)                                                                                           aTGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGt                                                       >  1:248375/1‑151 (MQ=255)                                                                                                   cccGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTAt                                               >  1:103486/1‑151 (MQ=255)                                                                                                     cGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATcc                                             >  1:656491/1‑151 (MQ=255)                                                                                                     cGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATcc                                             <  2:606512/151‑1 (MQ=255)                                                                                                      ggAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCg                                            >  1:174491/1‑151 (MQ=255)                                                                                                          acacGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAg                                        <  2:417383/151‑1 (MQ=255)                                                                                                           cacGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGc                                       >  2:666989/1‑151 (MQ=255)                                                                                                              gcgGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGt                                    <  2:107437/151‑1 (MQ=255)                                                                                                               cgGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAAGTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTa                                   <  1:439335/151‑1 (MQ=255)                                                                                                                ggCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTAt                                  >  2:71726/1‑151 (MQ=255)                                                                                                                   cTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTg                                >  2:437374/1‑150 (MQ=255)                                                                                                                        aGCGAGCATTCTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGAt                          <  2:436787/151‑1 (MQ=255)                                                                                                                                  cTTGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTGt                >  2:493947/1‑151 (MQ=255)                                                                                                                                   ttGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTGt                >  2:232537/1‑150 (MQ=255)                                                                                                                                   ttGATGCCCTGCCCGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTGt                >  2:155237/1‑150 (MQ=255)                                                                                                                                           ctgcccGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTg                 >  2:175291/1‑141 (MQ=255)                                                                                                                                           ctgcccGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTg                 <  1:175290/141‑1 (MQ=255)                                                                                                                                             gcccGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTACAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTGTAGGCCGGAt       >  1:159844/1‑149 (MQ=255)                                                                                                                                                cGCTAAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTGTAGGCCGGATAgatg  <  2:206188/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GAAGCAGAAAATCCCCAGCGCCTGTTCGCCCTCACCACGAAAGGTACACCTGCTCACAGCGCCGTAAGCTATCAGGATGGCGATTATCTGATGTTCGGCCCGGAAACACGCGGCCTGCCAGCGAGCATTCTTGATGCCCTGCCCGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTGTAGGCCGGATAGATG  >  NC_000913/3781416‑3781709

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