Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,781,688	G→A	*158* (TAG→TAA)	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,688	0	G	A	100.0%	198.9 / NA	61	*158* (TAG→TAA)	trmL	tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/34);  total (27/34)

CTGCCCGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTG  >  NC_000913/3781554‑3781823                                                                                                                                       |                                                                                                                                        ctgcccGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTAc                                                                                                                                     >  1:1601007/1‑139 (MQ=255)   gcccGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTg                                                                                                                                   >  2:575079/1‑139 (MQ=255)   gcccGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTg                                                                                                                                   <  2:515838/139‑1 (MQ=255)   gcccGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTg                                                                                                                                   >  1:1845726/1‑139 (MQ=255)    cccGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGt                                                                                                                                  <  1:1286943/139‑1 (MQ=255)     ccGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTa                                                                                                                                 <  2:1007411/139‑1 (MQ=255)       gCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAgg                                                                                                                               >  2:813672/1‑139 (MQ=255)          gAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGAGTAATTACTGTAGGCCg                                                                                                                            <  2:396494/139‑1 (MQ=255)            aCAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGTATGAATCTGTCCAATGCGGTGTCGGTCGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAACTACTGTAGGCCGGa                                                                                                                          <  2:1141115/139‑1 (MQ=255)            aCAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGa                                                                                                                          <  1:858382/139‑1 (MQ=255)              aaaaaaTTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATa                                                                                                                        <  2:1134427/139‑1 (MQ=255)                     tCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGAtt                                                                                                                 >  2:912355/1‑139 (MQ=255)                        cATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTg                                                                                                                                   <  1:52324/118‑1 (MQ=255)                        cATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTg                                                                                                                                   >  2:52324/1‑118 (MQ=255)                              gATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGa                                                                                                                          >  2:245824/1‑121 (MQ=255)                              gATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGa                                                                                                                          <  1:245824/121‑1 (MQ=255)                               aTGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATc                                                                                                       >  1:69615/1‑139 (MQ=255)                                 ggTGCCGGACAGCCGCAGCATGAATCTGTCCTATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCg                                                                                                     >  1:656107/1‑139 (MQ=255)                                  gTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCgg                                                                                                    <  2:573811/139‑1 (MQ=255)                                     ccGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                                 >  1:618908/1‑139 (MQ=255)                                     ccGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                                 <  2:1252595/139‑1 (MQ=255)                                     ccGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                                 <  2:69615/139‑1 (MQ=255)                                             ccGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAg                                                                                         <  1:813672/139‑1 (MQ=255)                                             ccGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAg                                                                                         >  2:428355/1‑139 (MQ=255)                                                 agcaTGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACa                                                                                     >  1:508870/1‑139 (MQ=255)                                                    aTGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTa                                                                                  <  1:575079/139‑1 (MQ=255)                                                    aTGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTa                                                                                  <  1:912355/139‑1 (MQ=255)                                                              ccAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGatcccatcc                                                                        <  2:508870/139‑1 (MQ=255)                                                               cAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAG‑TGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATcccc                                                                      <  2:1601007/139‑1 (MQ=255)                                                                     ggTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACt                                                                 <  2:656107/139‑1 (MQ=255)                                                                     ggTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACt                                                                 <  1:1790517/139‑1 (MQ=255)                                                                     ggTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACt                                                                 <  2:936259/139‑1 (MQ=255)                                                                       tgtCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCa                                                               >  2:837770/1‑139 (MQ=255)                                                                        gtCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCaa                                                              <  1:428355/139‑1 (MQ=255)                                                                                     tATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATa                                                                                                                        >  2:1709803/1‑68 (MQ=255)                                                                                     tATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATa                                                                                                                        <  1:1709803/68‑1 (MQ=255)                                                                                     tATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAg                                                                                         <  2:1462086/99‑1 (MQ=255)                                                                                     tATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAg                                                                                         >  1:1462086/1‑99 (MQ=255)                                                                                          aGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                                 >  1:973357/1‑86 (MQ=255)                                                                                          aGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                                 <  2:973357/86‑1 (MQ=255)                                                                                            ccTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAg                                                                                         <  1:178245/92‑1 (MQ=255)                                                                                            ccTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAg                                                                                         >  2:178245/1‑92 (MQ=255)                                                                                            ccTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTg                                          <  1:837770/139‑1 (MQ=255)                                                                                               ggcggcAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                                 >  1:1514565/1‑81 (MQ=255)                                                                                               ggcggcAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                                 <  2:1514565/81‑1 (MQ=255)                                                                                               ggcggcAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGaaa                                       <  2:1845726/139‑1 (MQ=255)                                                                                                   gcAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCt                                   <  1:1186203/139‑1 (MQ=255)                                                                                                                 cGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCAtt                     >  2:13632/1‑139 (MQ=255)                                                                                                                     aGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACt                                                                 <  2:969038/91‑1 (MQ=255)                                                                                                                     aGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACt                                                                 >  1:969038/1‑91 (MQ=255)                                                                                                                         gTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGc                                                                                                   >  1:319518/1‑53 (MQ=255)                                                                                                                         gTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGc                                                                                                   >  1:557834/1‑53 (MQ=255)                                                                                                                         gTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGc                                                                                                   <  2:557834/53‑1 (MQ=255)                                                                                                                         gTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGc                                                                                                   <  2:319518/53‑1 (MQ=255)                                                                                                                         gTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                                 >  1:419683/1‑55 (MQ=255)                                                                                                                         gTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                                 <  2:419683/55‑1 (MQ=255)                                                                                                                         gTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGCtt             >  2:63363/1‑139 (MQ=255)                                                                                                                                   attaattaCTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCt   >  1:213760/1‑139 (MQ=255)                                                                                                                                    ttAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACt                                                                 <  2:702081/76‑1 (MQ=255)                                                                                                                                    ttAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACt                                                                 >  1:702081/1‑76 (MQ=255)                                                                                                                                    ttAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTg  >  1:1668972/1‑139 (MQ=255)                                                                                                                                    ttAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTg  >  1:504141/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                        CTGCCCGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTG  >  NC_000913/3781554‑3781823

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