Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,824,561	G→A	R16Q (CGG→CAG)	trmH →	tRNA (Gm18) 2'‑O‑methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,824,561	0	G	A	100.0%	134.9 / NA	39	R16Q (CGG→CAG)	trmH	tRNA (Gm18) 2'‑O‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/22);  total (17/22)

GACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACC  >  NC_000913/3824428‑3824691                                                                                                                                      |                                                                                                                                   gACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGcc                                                                                                                              <  2:689646/140‑1 (MQ=255)        tgCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTg                                                                                                                       <  1:573987/140‑1 (MQ=255)             cTGGCGAATATCATGCGCAAAATCTGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGt                                                                                                                  <  1:265104/140‑1 (MQ=255)             cTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGt                                                                                                                  <  1:627224/140‑1 (MQ=255)              tGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTc                                                                                                                 >  1:57632/1‑140 (MQ=255)                  gAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCa                                                                                                             <  1:215312/140‑1 (MQ=255)                  gAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCa                                                                                                             <  2:162098/140‑1 (MQ=255)                          tGCGCAAAATCCGCGTGATGCCAGACGTGATTAGAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAg                                                                                                     >  1:260560/1‑140 (MQ=255)                          tGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGt                                                                                                                  <  2:699916/127‑1 (MQ=255)                          tGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGt                                                                                                                  >  1:699916/1‑127 (MQ=255)                          tGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAg                                                                                                     >  1:124833/1‑140 (MQ=255)                            cgcAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAg                                                                                                        <  2:421302/135‑1 (MQ=255)                            cgcAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAg                                                                                                        >  1:421302/1‑135 (MQ=255)                             gcAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTc                                                                                                  >  1:138906/1‑140 (MQ=255)                                  aTCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACaa                                                                                             <  1:577153/140‑1 (MQ=255)                                      gcgTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCt                                                                                         <  1:81172/140‑1 (MQ=255)                                      gcgTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCt                                                                                         <  1:246511/140‑1 (MQ=255)                                             gCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACg                                                                                  >  1:128553/1‑140 (MQ=255)                                                     tGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCg                                                                          <  2:124833/140‑1 (MQ=255)                                                          aaaGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGattat                                                                     <  2:160177/140‑1 (MQ=255)                                                                aCCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTAc                                                               >  2:691448/1‑140 (MQ=255)                                                                aCCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTAc                                                               >  1:92085/1‑140 (MQ=255)                                                                aCCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTAc                                                               >  1:317462/1‑140 (MQ=255)                                                                  ccgaaaccgaaaTTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCg                                                             <  2:493723/140‑1 (MQ=255)                                                                    gaaaccgaaaTTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCATAAACCTCATAACGTTTCTGCGATTATTCGTACCGCa                                                           <  1:61804/140‑1 (MQ=255)                                                                    gaaaccgaaaTTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCa                                                           <  2:186268/140‑1 (MQ=255)                                                                      aaccgaaaTTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGa                                                         <  2:57632/140‑1 (MQ=255)                                                                                aaTGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGc                                               <  1:319207/140‑1 (MQ=255)                                                                                  tGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGt                                             <  2:317462/140‑1 (MQ=255)                                                                                       tATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATg                                        >  1:296203/1‑140 (MQ=255)                                                                                         tGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGaa                                      >  1:294240/1‑140 (MQ=255)                                                                                                       tATGCACGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGGACATGAATTTCACGCCGTCTg                        >  2:500904/1‑140 (MQ=255)                                                                                                            aCGCATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTg                   <  1:691448/140‑1 (MQ=255)                                                                                                               cATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTa                >  1:263455/1‑140 (MQ=255)                                                                                                               cATCTGCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTa                >  1:504669/1‑140 (MQ=255)                                                                                                                    gCGAAATGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGc           <  2:379094/140‑1 (MQ=255)                                                                                                                         aTGCTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCTGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATgcg      <  1:154233/140‑1 (MQ=255)                                                                                                                            cTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCAc   >  2:244787/1‑140 (MQ=255)                                                                                                                            cTCGCCAGGCAGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCAc   >  1:641364/1‑140 (MQ=255)                                                                                                                             tCGCCAGGTAGCAGCCTGATCTGTCCGTTGGGATGGAGCGAGTCCGCAAACCTGATATCGTTTCTGCGATTATTCGTACCGCAGATGCGGTTGGCGTACATGAAGTTCTCGCCGGCTGGCCTGTTAGCCGCATGCGCGcc  <  2:100202/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   GACCGTGTGCATCTGGCGAATATCATGCGCAAAATCCGCGTGATGCCAGACGTGATTAAAGTCACCCGAAACCGAAATTAATGTTTTATGAACCCAACACGTTATGCACGCATCTGCGAAATGCTCGCCAGGCGGCAGCCTGATCTGACCGTCTGCATGGAGCAGGTCCACAAACCTCATAACGTTTCTGCGATTATTCGTACCGCAGATGCCGTTGGCGTACATGAAGTTCACGCCGTCTGGCCTGGTAGCCGCATGCGCACC  >  NC_000913/3824428‑3824691

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