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%	162.3 / NA	52	*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 (24/28);  total (24/28)

GCCCGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCG  >  NC_000913/3781556‑3781816                                                                                                                                     |                                                                                                                                 gcccGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTg                                                                                                                            >  1:1096086/1‑139 (MQ=255)     gCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAgg                                                                                                                        <  2:1096086/139‑1 (MQ=255)                  ttCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTCGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAgatgat                                                                                                           <  2:1023629/139‑1 (MQ=255)                   tCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGAtt                                                                                                          >  2:1041642/1‑139 (MQ=255)                       aTTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACAt                                                                                                      <  1:56075/139‑1 (MQ=255)                                   ccGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                          <  2:191425/139‑1 (MQ=255)                                     ggACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGa                                                                                                                   >  1:944921/1‑112 (MQ=255)                                     ggACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGGATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCAGa                                                                                                                   <  2:944921/112‑1 (MQ=255)                                          gCCGCAGCATGAGTCTGTCCAATGCGGTGTCGGTAGAGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATcaca                                                                                   <  1:698445/139‑1 (MQ=255)                                           ccGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAg                                                                                  >  1:493320/1‑139 (MQ=255)                                            cGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAgg                                                                                 <  2:241064/139‑1 (MQ=255)                                                 caTGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACAtt                                                                            <  2:568439/139‑1 (MQ=255)                                                       tCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGatcc                                                                      <  2:392283/139‑1 (MQ=255)                                                              aaTGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATcccc                                                               <  1:378393/139‑1 (MQ=255)                                                              aaTGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATcccc                                                               <  1:1041642/139‑1 (MQ=255)                                                                 gCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATa                                                            <  2:1076408/139‑1 (MQ=255)                                                                        cGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAAt                                                     <  2:493320/139‑1 (MQ=255)                                                                              tggtgTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACt                                                                                                                             <  1:691512/61‑1 (MQ=255)                                                                              tggtgTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACt                                                                                                                             >  2:691512/1‑61 (MQ=255)                                                                                     tGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGa                                                                                                                   <  2:231215/64‑1 (MQ=255)                                                                                     tGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGa                                                                                                                   >  1:231215/1‑64 (MQ=255)                                                                                     tGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGatcccatcc                                                                 >  1:108166/1‑114 (MQ=255)                                                                                     tGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGatcccatcc                                                                 <  2:108166/114‑1 (MQ=255)                                                                                        aGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACt                                                                                                                             >  2:744265/1‑51 (MQ=255)                                                                                        aGCCTGGCGGCAGATGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACt                                                                                                                             <  1:744265/51‑1 (MQ=38)                                                                                         gCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGtt                                    >  1:641709/1‑139 (MQ=255)                                                                                          ccTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACt                                                                                                                             <  2:964971/49‑1 (MQ=255)                                                                                          ccTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACt                                                                                                                             >  1:964971/1‑49 (MQ=255)                                                                                          ccTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATcaca                                                                                   >  1:643271/1‑91 (MQ=255)                                                                                          ccTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATcaca                                                                                   <  2:643271/91‑1 (MQ=255)                                                                                             ggcggcAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAgat                                                                                                              <  1:336229/61‑1 (MQ=255)                                                                                             ggcggcAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAgat                                                                                                              >  2:336229/1‑61 (MQ=255)                                                                                             ggcggcAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                          >  1:912273/1‑81 (MQ=255)                                                                                             ggcggcAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                          <  2:261915/81‑1 (MQ=255)                                                                                             ggcggcAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                          <  2:912273/81‑1 (MQ=255)                                                                                             ggcggcAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCAc                                                                                          >  1:261915/1‑81 (MQ=255)                                                                                               cggcAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATg                              <  2:625664/139‑1 (MQ=255)                                                                                                            aTCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATcac                                                                                    >  2:575077/1‑72 (MQ=255)                                                                                                            aTCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATcac                                                                                    <  1:575077/72‑1 (MQ=255)                                                                                                             tCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTaa                                           <  2:495758/112‑1 (MQ=255)                                                                                                             tCCGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTaa                                           >  1:495758/1‑112 (MQ=255)                                                                                                               cGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACAtt                                                                            >  1:801455/1‑77 (MQ=255)                                                                                                               cGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACAtt                                                                            <  2:801455/77‑1 (MQ=255)                                                                                                               cGGGAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCAtt              >  2:748988/1‑139 (MQ=255)                                                                                                                  gAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGAt                                                                                       <  2:1100388/63‑1 (MQ=255)                                                                                                                  gAGCGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGAt                                                                                       >  1:1100388/1‑63 (MQ=255)                                                                                                                     cGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCACATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGc        >  2:282773/1‑139 (MQ=255)                                                                                                                     cGGTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGc        >  1:307363/1‑139 (MQ=255)                                                                                                                       gTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGCtt      >  2:184635/1‑139 (MQ=255)                                                                                                                       gTATTGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGCtt      >  2:627838/1‑139 (MQ=255)                                                                                                                           tGAGAGATTAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCg  <  1:627838/139‑1 (MQ=255)                                                                                                                                  ttAATTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCatat                   >  1:151424/1‑115 (MQ=255)                                                                                                                                     |                                                                                                                                 GCCCGCTGAACAAAAAATTCGCATTCCGATGGTGCCGGACAGCCGCAGCATGAATCTGTCCAATGCGGTGTCGGTAGTGGTGTATGAAGCCTGGCGGCAGTTGGGGTATCCGGGAGCGGTATTGAGAGATTAGTTACTGTAGGCCGGATAGATGATTACATCGCATCCGGCACGATCACAGGACATTAGATCCCATCCCCATACTCAAATGTATGGTTAATACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCG  >  NC_000913/3781556‑3781816

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