Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,723,735	T→C	E221E (GAA→GAG)	glyS ←	glycine tRNA synthetase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,723,735	0	T	C	100.0%	183.3 / NA	48	E221E (GAA→GAG)	glyS	glycine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/23);  total (25/23)

AGGAATTTCTCTTCGAATTTTGCGGTCAGAACGACCGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGTTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGATCGGACTGAATGCCCAGAATGGTTGCCGGAAT  >  NC_000913/3723589‑3723881                                                                                                                                                   |                                                                                                                                                   aGGAATTTCTCTTCGAATTTTGCGGTCAGAACGACCGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGctctt                                                                                                                                                >  2:280230/1‑151 (MQ=255)   gAATTTCTCTTCGAATTTTGCGGTCAGAACGCCCGGCCACTCCCCCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCg                                                                                                                                              <  1:659794/151‑1 (MQ=255)           ctTCGAATTTTGCGGTCAGAACGACCGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCg                                                                                                                                        <  2:117167/149‑1 (MQ=255)             tCGAATTTTGCGGTCAGAACGACCGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCgg                                                                                                                                       <  1:489476/148‑1 (MQ=255)             tCGAATTTTGCGGTCAGAACGACCGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCg                                                                                                                                     <  2:205676/150‑1 (MQ=255)             tCGAATTTTGCGGTCAGAACGACCGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCg                                                                                                                                     <  2:576460/150‑1 (MQ=255)                         gTCAGAACGACCGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCAc                                                                                                                        <  2:184415/151‑1 (MQ=255)                           cAGAACGACCGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGc                                                                                                                      >  2:558085/1‑151 (MQ=255)                              aaCGACCGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTc                                                                                                                    <  1:316043/150‑1 (MQ=255)                                 gACCGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGc                                                                                                                >  2:158656/1‑151 (MQ=255)                                   ccGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAg                                                                                                              <  2:148012/151‑1 (MQ=255)                                    cGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGa                                                                                                             <  2:98741/151‑1 (MQ=255)                                        cACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAAttt                                                                                                         >  1:128685/1‑151 (MQ=255)                                          ccccaccaGCGATGCCACTTCTTCGAGCAGGCGTTCGCTTAATACATCGTTCCCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTc                                                                                                        <  1:146413/148‑1 (MQ=255)                                             caccaGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCgg                                                                                                     >  1:307830/1‑150 (MQ=255)                                                 aGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACt                                                                                                >  2:434528/1‑151 (MQ=255)                                                   cGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTg                                                                                               >  2:447604/1‑150 (MQ=255)                                                        ccACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGc                                                                                         >  2:489520/1‑151 (MQ=255)                                                         cACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCg                                                                                        >  2:221471/1‑151 (MQ=255)                                                             tcttcCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGtt                                                                                      <  2:122247/149‑1 (MQ=255)                                                                      agGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTg                                                                            >  2:458953/1‑150 (MQ=255)                                                                        gCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAAt                                                                         >  2:252854/1‑151 (MQ=255)                                                                           ttCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTc                                                                       <  2:632799/150‑1 (MQ=255)                                                                             cGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCgg                                                                    >  2:631045/1‑151 (MQ=255)                                                                                ttAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTc                                                                 >  2:348636/1‑151 (MQ=255)                                                                                  aaGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGc                                                               <  1:403459/151‑1 (MQ=255)                                                                                      cAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCAt                                                           <  1:443717/151‑1 (MQ=255)                                                                                       aGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATa                                                          <  1:120782/151‑1 (MQ=255)                                                                                           ttACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGc                                                      >  1:682450/1‑151 (MQ=255)                                                                                            tACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGc                                                      <  2:422884/150‑1 (MQ=255)                                                                                              ccGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCGCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCgg                                                    >  1:410085/1‑150 (MQ=255)                                                                                                 ccAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGc                                                >  2:625264/1‑151 (MQ=255)                                                                                                   aaTCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCcg                                              >  2:322766/1‑151 (MQ=255)                                                                                                    aTCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCcgc                                             >  2:16111/1‑151 (MQ=255)                                                                                                     tCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCcgcg                                            >  1:101901/1‑151 (MQ=255)                                                                                                          cgcgCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCa                                       <  2:495165/151‑1 (MQ=255)                                                                                                                gctTCTTCGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGAt                                 <  2:397302/151‑1 (MQ=255)                                                                                                                  ttcctcGGCATCGGCTTTAATCTTCGCCTTAAGCTCCTCATACTCGTCGATGACTTTCCCACGCTCACGCAGATTTACCGGATACTGATCGGCGTTATCGATGGTGAATTCCGCCTCTCCCATAAACCGACGGCCGCGACTCACGCGATc                                >  1:424881/1‑150 (MQ=255)                                                                                                                   tcttcGGCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGATc                                <  1:508579/149‑1 (MQ=255)                                                                                                                        ggCATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGATCGGACTGa                         <  2:385937/151‑1 (MQ=255)                                                                                                                          cATCGGCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGATCGGACTGAAt                       >  1:82961/1‑151 (MQ=255)                                                                                                                               gCTTTAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGATCGGACTGAATGccc                   >  1:250582/1‑150 (MQ=255)                                                                                                                                  ttAATCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGATCGGACTGAATGCCCAGa                >  2:168019/1‑150 (MQ=255)                                                                                                                                    aaTCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGATCGGACTGAATGCCCAGAATg             <  1:416994/151‑1 (MQ=255)                                                                                                                                     aTCTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGACTTTCCGGATACTGATCGCCGTTATCGCTGGTGAATTCCGGCTCGCCCATAAAGCGCTGGCCGCGAATCACGCGATCGGACTGAATGCCCAGaa               >  1:684917/1‑148 (MQ=255)                                                                                                                                       cTTCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGATCGGACTGAATGCCCAGAATGGtt          >  2:35549/1‑151 (MQ=255)                                                                                                                                         tCGCCTTACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGATCGGACTGAATGCCCAGAATGGtt          <  1:321910/149‑1 (MQ=255)                                                                                                                                               tACGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGATCGGACTGAATGCCCAGAATGGTTGCCGGAAt  >  2:276216/1‑151 (MQ=255)                                                                                                                                                aCGCTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGATCGGACTGAATGCCCAGAATGGTTGCCGGAAt  <  2:20192/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AGGAATTTCTCTTCGAATTTTGCGGTCAGAACGACCGGCCACTCCACCAGCGAAGCCACTTCTTCCAGCAGGCTTTCGCTTAAGTCAGCGTTACCGCCAATCTTACGCGCTGCTTCTTCGGCATCGGCTTTAATCTTCGCCTTACGTTCTTCGTAATCGGCGATGACTTTCCCACGCTCACGCAGAATTTCCGGATACTGATCGGCGTTATCGATGGTGAATTCCGGCTCGCCCATAAAGCGGTGGCCGCGAATCACGCGATCGGACTGAATGCCCAGAATGGTTGCCGGAAT  >  NC_000913/3723589‑3723881

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