Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,532	C→T	I70I (ATC→ATT)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,532	0	C	T	100.0%	108.0 / NA	33	I70I (ATC→ATT)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/14);  total (19/14)

AGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGT  >  NC_000913/4617287‑4617778                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       aGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTCAAAGAAAGCAGCCTGCGTGCACTGCAATTGATGGACCTGACCAACCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATAAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAACGAGCAGGGCACCCCGGAAATCCGTATTGCTAc                                                                                                                                                                                                                                                   >  1:172401/1‑251 (MQ=255)                     aaGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACg                                                                                                                                                                                                                               >  2:134937/1‑251 (MQ=255)                         cAGGAGAATGAAATGACTGATCTAAAAGCAAGCAGCCGGCGTGCACTAAAATTGATGGACCTGACCCCCCTAAATGACGACGACACCGACGAAAAAGTAATTCCTCTGTGTCATCAGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTaa                                                                                                                                                                                                                           <  1:23897/251‑1 (MQ=255)                                atgaaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGACCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCCCCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACgatgac                                                                                                                                                                                                                    >  1:130476/1‑251 (MQ=255)                                           aatcggaaagcaaaaggcatGCGGGCACTGAAATTGATGGACCTGCCCCCCCTAAAAGACGCCACCCCCGACAAGAAAGTAATGGCTCTGTGTCAACAGGCAAAACCCCCGGCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACATCCCACACGGTAACGATgacatcgaaatc                                                                                                                                                                                                           <  2:136003/230‑1 (MQ=255)                                                 aaaGCAACAAGCCTGCGTGCCGGGAAATTGATGGACCTGCCCCCCCTGAATGACAACGACACCGCCGAGAAAGTAATTGCTCTGTTTCACCAGCCCAAACCCCCGGCCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGc                                                                                                                                                                                                   <  2:15109/251‑1 (MQ=255)                                                        gcagcCTGCGTGCACTGAAATTGCTGGACCTGACCCCCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTCTCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCACAACACTGAAAGCGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGCCATCGCCATCGCGCTTGCAGAAAcc                                                                                                                                                                                            >  1:12427/1‑251 (MQ=255)                                                              tGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGCACCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCg                                                                                                                                                                                      >  1:226463/1‑251 (MQ=255)                                                                    cACTGAAAAGGATGGACCGGCCAACCCGGAATGCCGACGACACCGACTAGAAAGTAATTGCTCGGTGCCATCAGGCCAACTCCCGGGGCGCCAAAACCGCCGCTATCTGTATCTATCCGCGCGTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACAGCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                                                                                                                                                                                <  2:12427/251‑1 (MQ=255)                                                                                  gACCTGACCCCCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                                                                                                                                                                                              <  1:223821/223‑1 (MQ=255)                                                                                  gACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                                                                                                                                                                                              >  2:223821/1‑223 (MQ=255)                                                                                  gACCTGACAACCCTGAATGACGACGACACCGAAGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGCCGGCAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGACAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAAACAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGccgg                                                                                                                                                                  >  2:163162/1‑248 (MQ=255)                                                                                       gACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATACGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCcgatga                                                                                                                                                               >  1:41477/1‑249 (MQ=255)                                                                                              cgaaaTGACGACAACACCGACGAAAAAGTAATTGCTCTGGGTCATCAGGCAAAACCCCCGGTCGGCAAAACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGt                                                                                                                                                      <  1:130736/248‑1 (MQ=255)                                                                                                                       aagtaattgctcTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATg                                                                                                                             >  1:28562/12‑251 (MQ=255)                                                                                                                            attgctcTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGCAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCggg                                                                                                                        >  2:201688/7‑251 (MQ=255)                                                                                                                                    gtgtCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGcc                                                                                                                                                                    >  1:65456/1‑198 (MQ=255)                                                                                                                                       tcatcaGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGAACGAGCAgg                                                                                                             >  2:267571/1‑251 (MQ=255)                                                                                                                                                     ccccGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTg                                                                                               <  1:134937/249‑1 (MQ=255)                                                                                                                                                                                gTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTg                                                                        <  1:103105/247‑1 (MQ=255)                                                                                                                                                                                         cTCGCTTTATCCCGATTGCTCGCAAAACACTGACAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGGAAAGAAAAGGGTGGTTTTGAACTGGTAAAAGCCTGTTAAGAGGCTTTCGGGGCaacgga                                                            >  2:244635/1‑245 (MQ=255)                                                                                                                                                                                          tCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAAAAGCTTGCGCGGCcgcggagg                                                          >  2:269486/1‑247 (MQ=255)                                                                                                                                                                                            gCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGCATGt                                                         >  1:124166/1‑250 (MQ=255)                                                                                                                                                                                                   cccGATTGCTCGCAAAACACTGAAAGAGCAGGGCCCCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTActgctg                                                  <  1:235736/250‑1 (MQ=255)                                                                                                                                                                                                      gATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACCTCCCACACGGTAACGATGACATCGAAATCGCTCTGGCAGAAACCCGTGCGGCAATCGCCTACGGGGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTTATTGCGGGTAACGAGCAGGTTGGTTTTGACCTGGGGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAg                                              <  1:109881/251‑1 (MQ=255)                                                                                                                                                                                                             cGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGatc                                         <  1:11693/249‑1 (MQ=255)                                                                                                                                                                                                                 aaaCACTGAAAGAGCAGGGCACCCCGGAAATCCGTCTTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGc                                                                                                                                                                     >  2:29344/1‑121 (MQ=255)                                                                                                                                                                                                                 aaaCACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGc                                                                                                                                                                     <  1:29344/121‑1 (MQ=255)                                                                                                                                                                                                                   acactGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAAc                                 <  1:85341/248‑1 (MQ=255)                                                                                                                                                                                                                    cactGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAAcc                                <  1:185852/249‑1 (MQ=255)                                                                                                                                                                                                                                              aTCCGTATTGCTACGGTAACCAACTTCCCAAACGGTAACGATGACCTAGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGACTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGACACGAGCAGGTTGGTCTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCCGCGGATGTACTGCTGAAAGGGCTCATCGAAACCGGCGAACTGAAAAACGAAGCGCTGAt      >  1:237643/1‑251 (MQ=255)                                                                                                                                                                                                                                                ccGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGGTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTAATGGTGAAAGTGATCATAGAAAACGGCGAACTTAAAGACGAAGCGCGGATcc    >  2:222595/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gTATTGCTACGGTAACCAAATTCCCACGCGGTAACGATGGCATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGATCAGGTTGGTTTTTACCTGGGGAAAGCCTGGAAAGGGGATTGGGCGGCCGCGAAGGTACTGATGGAAGTGTTCATAGAAAACGGGGAACTGAAAAAAGAAGCGCTGagtcgt  >  2:96969/1‑246 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       AGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGT  >  NC_000913/4617287‑4617778

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