Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,824,905	G→T	A131S (GCC→TCC)	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,905	0	G	T	100.0%	211.2 / NA	60	A131S (GCC→TCC)	trmH	tRNA (Gm18) 2'‑O‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (24/36);  total (24/36)

TCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAA  >  NC_000913/3824769‑3825036                                                                                                                                         |                                                                                                                                    tctcAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTccc                                                                                                                                  <  2:869203/140‑1 (MQ=255)    cAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTgg                                                                                                                               <  2:486759/140‑1 (MQ=255)     aaaGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGc                                                                                                                              >  2:777060/1‑140 (MQ=255)         gCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGAt                                                                                                                          <  2:167532/140‑1 (MQ=255)          ccAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATc                                                                                                                         <  1:432872/140‑1 (MQ=255)          ccAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATc                                                                                                                         >  2:169690/1‑140 (MQ=255)            aGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAg                                                                                                                       >  1:194056/1‑140 (MQ=255)            aGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAg                                                                                                                       <  2:797802/140‑1 (MQ=255)              ggCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGa                                                                                                                     <  2:690882/140‑1 (MQ=255)              ggCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGa                                                                                                                     <  1:213080/140‑1 (MQ=255)              ggCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGa                                                                                                                     <  2:324471/140‑1 (MQ=255)                 aTGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGAcat                                                                                                                  <  1:457747/140‑1 (MQ=255)                   gCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGAcatca                                                                                                                >  1:329163/1‑140 (MQ=255)                   gCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGAcatca                                                                                                                >  2:859862/1‑140 (MQ=255)                   gCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGAcatca                                                                                                                <  1:894449/140‑1 (MQ=255)                         tCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTc                                                                                                          >  2:254401/1‑140 (MQ=255)                         tCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTc                                                                                                          >  1:328445/1‑140 (MQ=255)                           tGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCg                                                                                                        <  2:908346/140‑1 (MQ=255)                                aCCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCgatgat                                                                                                   <  1:212270/140‑1 (MQ=255)                                aCCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCgatgat                                                                                                   <  2:328445/140‑1 (MQ=255)                                    aTCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGc                                                                                               <  1:506805/140‑1 (MQ=255)                                     tCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCa                                                                                              <  2:194056/140‑1 (MQ=255)                                     tCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCa                                                                                              <  2:895966/140‑1 (MQ=255)                                     tCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCa                                                                                              >  1:873268/1‑140 (MQ=255)                                           tGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACTGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGc                                                                                        <  2:558025/140‑1 (MQ=255)                                                 cGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGc                                                                                  <  1:876812/140‑1 (MQ=255)                                                  gCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCt                                                                                 <  1:254401/140‑1 (MQ=255)                                                       cGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATg                                                                            >  2:377092/1‑140 (MQ=255)                                                                           acacTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTc                                                        <  2:667330/140‑1 (MQ=255)                                                                               tCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTc                                                    >  1:111460/1‑140 (MQ=255)                                                                                 gCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCa                                                              <  2:590230/128‑1 (MQ=255)                                                                                 gCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCa                                                              >  1:590230/1‑128 (MQ=255)                                                                                 gCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCtt                                                  >  1:43123/1‑140 (MQ=255)                                                                                 gCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCtt                                                  >  2:324142/1‑140 (MQ=255)                                                                                  cccGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCttt                                                 >  2:650139/1‑140 (MQ=255)                                                                                   ccGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTa                                                >  1:368282/1‑140 (MQ=255)                                                                                    cGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTAc                                               <  2:626308/140‑1 (MQ=255)                                                                                     gACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACg                                              <  1:667658/140‑1 (MQ=255)                                                                                      aCCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGa                                             >  1:224913/1‑140 (MQ=255)                                                                                      aCCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGa                                             >  1:869582/1‑140 (MQ=255)                                                                                           cATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGccc                                        >  1:103033/1‑140 (MQ=255)                                                                                              tttGATGGGACAGGAGAAAACGAGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCAATCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGACTCAGCCCTCATTCTTTACGAAGCCCAGc                                     <  2:714836/140‑1 (MQ=255)                                                                                                     ggACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCg                              >  1:20013/1‑140 (MQ=255)                                                                                                     ggACAGGAGAAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCg                              >  1:771783/1‑140 (MQ=255)                                                                                                             gaAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATg                      <  2:224913/140‑1 (MQ=255)                                                                                                             gaAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATg                      <  1:324142/140‑1 (MQ=255)                                                                                                             gaAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATg                      <  2:368282/140‑1 (MQ=255)                                                                                                             gaAAACGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATg                      <  2:43123/140‑1 (MQ=255)                                                                                                              aaaaCGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGc                     <  1:169690/140‑1 (MQ=255)                                                                                                              aaaaCGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGc                     <  2:869582/140‑1 (MQ=255)                                                                                                              aaaaCGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGc                     <  2:87797/140‑1 (MQ=255)                                                                                                              aaaaCGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGc                     <  2:179701/140‑1 (MQ=255)                                                                                                              aaaaCGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGc                     <  2:111460/140‑1 (MQ=255)                                                                                                               aaaCGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGc                     >  1:229101/1‑139 (MQ=255)                                                                                                               aaaCGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGc                     <  2:229101/139‑1 (MQ=255)                                                                                                                aaCGGGCATCACGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAg                   >  1:502183/1‑140 (MQ=255)                                                                                                                          aCGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCt         >  2:600933/1‑140 (MQ=255)                                                                                                                           cGCAGGAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTg        >  2:645155/1‑140 (MQ=255)                                                                                                                               ggAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTg    <  2:329163/140‑1 (MQ=255)                                                                                                                                gAAGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGa   <  1:377092/140‑1 (MQ=255)                                                                                                                                 aaGCATTGTCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGaa  <  1:179210/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    TCTCAAAGGCCAGGGCATGCAGATTCTGGCAACCCATCTTTCTGATAACGCTGTCGATTTCCGCGAAATTGATTACACTCGCCCGACCTGCATTTTGATGGGACAGGAGAAAACGGGCATCACGCAGGAAGCATTGGCCCTGGCGGATCAGGACATCATCATTCCGATGATCGGCATGGTGCAGTCGCTGAATGTTTCCGTTGCCTCAGCCCTCATTCTTTACGAAGCCCAGCGTCAGCGGCAAAATGCAGGCATGTACCTGCGTGAA  >  NC_000913/3824769‑3825036

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