Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,565	C→T	D81D (GAC→GAT)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,565	0	C	T	100.0%	97.5 / NA	32	D81D (GAC→GAT)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/13);  total (19/13)

ATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGG  >  NC_000913/4617318‑4617807                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   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)                                                                                                                                                                                                                              gTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGACACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGGTCATCGAAACCGGCGAAATGAAAGGAGAAGCGCTGATCCGTAAAGCGTCTg                     >  1:170452/1‑250 (MQ=255)                                                                                                                                                                                                                                          ccACACGGTAACGCTGACATCGAAATCGCGCTGGCAGAAACCAGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGAAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCCAAACCGGCCAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATATACATCa         >  1:55223/1‑250 (MQ=255)                                                                                                                                                                                                                                                ggTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAAAGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGAACATCGAAACCGGCGAACTGAAAGAAGAAGGGCTGATCCGTAAAGCGTCTGAAATCTCCAACAAAGCggg  >  1:176121/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   ATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGACCTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGG  >  NC_000913/4617318‑4617807

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