Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,032,972	C→T	T36T (ACG→ACA)	yedJ ←	putative HD superfamily phosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,032,972	0	C	T	100.0%	186.1 / NA	50	T36T (ACG→ACA)	yedJ	putative HD superfamily phosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/33);  total (17/33)

AGGCGACGCGTCTCTTCTGCCGCCAGGATTGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGT  >  NC_000913/2032829‑2033116                                                                                                                                                |                                                                                                                                                 aGGCGACGCGTCTCTTCTGCCGCCAGGATTGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGccc                                                                                                                                           >  1:533431/1‑151 (MQ=255)          gTCTCTTCTGCCGCCAGGATTGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACgcggc                                                                                                                                  <  1:451279/151‑1 (MQ=255)            ctctTCTGCCGCCAGGATTGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACgcggcg                                                                                                                                 >  2:522909/1‑150 (MQ=255)             tcttctGCCGCCAGGATTGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGa                                                                                                                                >  1:292502/1‑150 (MQ=255)              cttctGCCGCCAGGATTGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACgcggcg                                                                                                                                 <  1:742907/148‑1 (MQ=255)                tctGCCGCCAGGATTGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGaaa                                                                                                                              <  2:458660/149‑1 (MQ=255)                     cgccAGGATTGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGacac                                                                                                                       >  2:651161/1‑151 (MQ=255)                         aGGATTGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACg                                                                                                                    >  1:660426/1‑150 (MQ=255)                         aGGATTGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGt                                                                                                                   <  1:320163/151‑1 (MQ=255)                             ttGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGa                                                                                                                                        <  1:447391/126‑1 (MQ=255)                             ttGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGa                                                                                                                                        >  2:447370/1‑126 (MQ=255)                               gaagaaCGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATggg                                                                                                             <  2:368899/151‑1 (MQ=255)                               gaagaaCGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATggg                                                                                                             <  2:432699/151‑1 (MQ=255)                               gaagaaCGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATgg                                                                                                              <  1:653816/150‑1 (MQ=255)                                   aaCGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCcgc                                                                                                         <  2:757050/151‑1 (MQ=255)                                   aaCGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCcgc                                                                                                         <  1:348242/151‑1 (MQ=255)                                         gCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCt                                                                                                    <  2:484732/150‑1 (MQ=255)                                         gCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCt                                                                                                    <  2:576225/150‑1 (MQ=255)                                               gCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCtg                                                                                             <  2:636740/151‑1 (MQ=255)                                               gCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCt                                                                                              <  2:180614/150‑1 (MQ=255)                                               gCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCt                                                                                              <  2:476458/150‑1 (MQ=255)                                                  gatgatTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCtgatg                                                                                          <  2:50688/151‑1 (MQ=255)                                                  gatgatTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCtgat                                                                                           >  2:727233/1‑150 (MQ=255)                                                     gatTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCcg                                                                                                          <  1:355834/132‑1 (MQ=255)                                                     gatTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCcg                                                                                                          >  2:355819/1‑132 (MQ=255)                                                     gatTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGatgatg                                                                                       >  2:620316/1‑151 (MQ=255)                                                          ttCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGttt                                                                                   <  1:667336/150‑1 (MQ=255)                                                           tCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGtttt                                                                                  <  1:634241/150‑1 (MQ=255)                                                                           atCGTGAAAATAACAGGCGGGTAAAATCACCAGAATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAAACGCCGCGAAAATGACACACGTCATGTGCCGCGTACTGGTGATGATGATGGATTTTCAACCAGTTTTCGAAc                                                                 >  2:762674/1‑151 (MQ=255)                                                                             cGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACt                                                                <  2:388135/150‑1 (MQ=255)                                                                                   aaTAACAGGCGGTTAAAATCACCAGCATATCGAAATCGTCATCTGCCGCCAGCTTCTGTGCTTTAGCCCAAACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGGGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGtgcc                                                           >  2:426072/1‑149 (MQ=255)                                                                                      aaCAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGtgcctgcc                                                       >  1:283847/1‑150 (MQ=255)                                                                                          ggCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAgtgt                                                  >  1:213397/1‑151 (MQ=255)                                                                                          ggCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAgtgt                                                  >  1:119188/1‑151 (MQ=255)                                                                                              tttAAAATCACCAGCATATCTACAGCGTCATCTGCCGCCAGCTTCTGTGCTTTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGCCCTGGTGCTGATGATGGTTTTGCACCCAGTTTACGAACTGTGCCTGCCAGtgttgt                                               <  2:462767/149‑1 (MQ=255)                                                                                              gTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTa                                              <  1:134038/151‑1 (MQ=255)                                                                                                  aaaTCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTc                                          <  2:241104/151‑1 (MQ=255)                                                                                                       aCCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATAc                                     >  2:570914/1‑151 (MQ=255)                                                                                                           gCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACgcgc                                 <  1:745389/151‑1 (MQ=255)                                                                                                            cATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCt                                <  1:319399/151‑1 (MQ=255)                                                                                                                       cGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCaaa                     <  1:207808/151‑1 (MQ=255)                                                                                                                       cGTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCaa                      >  1:29340/1‑150 (MQ=255)                                                                                                                        gTCATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCATTGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCaaa                     <  2:561660/150‑1 (MQ=255)                                                                                                                          cATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAagag                  <  1:481431/151‑1 (MQ=255)                                                                                                                          cATCTGCCGCCAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAag                    >  2:183992/1‑149 (MQ=255)                                                                                                                               gccgccAGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGc              <  2:638030/150‑1 (MQ=255)                                                                                                                                     aGCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTa       <  2:169189/151‑1 (MQ=255)                                                                                                                                      gCTTCTGTGCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAg      <  1:694354/151‑1 (MQ=255)                                                                                                                                            gtgCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAg   <  2:539562/148‑1 (MQ=255)                                                                                                                                            gtgCTGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGt  <  1:337381/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 AGGCGACGCGTCTCTTCTGCCGCCAGGATTGAAGAACGCTGCCTTTGCGGATGATTTTTCGCCAGGCTCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGT  >  NC_000913/2032829‑2033116

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: ‑