Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,436	T→C	T38T (ACT→ACC)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,436	0	T	C	100.0%	66.5 / NA	28	T38T (ACT→ACC)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/9);  total (19/9)

CATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGG  >  NC_000913/4617191‑4617671                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            cATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACCCCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGt                                                                                                                                                                                                                                        >  1:136003/1‑251 (MQ=255)                   ttGTACGTAGATATCCTTCATTGTGCCGTATCTCGCAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGCATTTTATTGCCGACAAGCCAGGAGAATGAAACGCCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATCGATGGACCTGACCACCCTGAATGACGCCGACACCGACGAGAAAGTAATTGCT‑CTATGTCATCACGCCCAAACCCCGGTCGGCACTACCGCCGCTAt                                                                                                                                                                                                                      >  2:109881/1‑251 (MQ=255)                    tGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATc                                                                                                                                                                                                                     >  2:23897/1‑251 (MQ=255)                      tAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTg                                                                                                                                                                                                                   >  2:103105/1‑251 (MQ=255)                                          gaaggtaaaTCGAAGTGGGTTGCGGTGAAGAATTGAGAAAAATAAAAAACTCGAAAGGTGATTTTGATGGGCGAAAAGCCAGGAAAATGAAATTACTGATCTAAAAGCAAACAGCCTGCGTGCACTGAAATTGATGGACCTGACCCCCCTAAAGAACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGCCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATcc                                                                                                                                                                                               <  1:185664/243‑1 (MQ=255)                                            aTGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCg                                                                                                                                                                                             >  1:176271/1‑251 (MQ=255)                                                       ggtgtgtTGCGGAGTAGATGTTAAAAAACTAACAAACTCGCAAGTTAAATTTTATTGGCGACAAGCCAGGAGAATAAAATGACTAATCTGAAAGCAAGCAGCCGGCGTGCATTGAAATTGAGGGACCGGCCCACCCTAAATGACGCCGACACCGACGAGAAAGT‑ATTTCCTCTGTTTCATCAGGCCAAACCCCCGTTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCaa                                                                                                                                                                                  <  1:114402/250‑1 (MQ=255)                                                                       gATGTTAGAATACTAACAAACTCGCAAGGTGATTTTTATTGGCGACAAGCCAGGAGCATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGc                                                                                                                                                                    >  2:185852/1‑249 (MQ=255)                                                                                       aaaaCTCGAAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATTACTGATCTGAAAGCAATAAGCCTGCGTGCACTGAAATTGATGGACCTGACCCCCCTGAATGACGACGACCCCGACGAGAAAGTAATTGTT‑TTGTGTCATCGGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATcc                                                                                                                                                  <  1:21014/250‑1 (MQ=255)                                                                                                 aGGTGAATTTTCTTGGCGACACGCCAGGAGAATGAAATGACGGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGCCACCGACGAGAATGTAATTGAT‑CTGTGTCATCAGGCACAAACCCCGGTCGACAATACCGCCGCTATCTGTATCAATACTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGCCAATCCGTAttgatac                                                                                                                                        >  2:197886/1‑245 (MQ=255)                                                                                                 aGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTCAAAGAAAGCAGCCTGCGTGCACTGCAATTGATGGACCTGACCAACCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATAAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAACGAGCAGGGCACCCCGGAAATCCGTATTGCTAc                                                                                                                                        >  1:172401/1‑251 (MQ=255)                                                                                                                     aaGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACg                                                                                                                    >  2:134937/1‑251 (MQ=255)                                                                                                                         cAGGAGAATGAAATGACTGATCTAAAAGCAAGCAGCCGGCGTGCACTAAAATTGATGGACCTGACCCCCCTAAATGACGACGACACCGACGAAAAAGTAATTCCT‑CTGTGTCATCAGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTaa                                                                                                                <  1:23897/251‑1 (MQ=255)                                                                                                                                atgaaatgaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGACCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCCCCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACgatgac                                                                                                         >  1:130476/1‑251 (MQ=255)                                                                                                                                           aatcggaaagcaaaaggcatGCGGGCACTGAAATTGATGGACCTGCCCCCCCTAAAAGACGCCACCCCCGACAAGAAAGTAATGGCT‑CTGTGTCAACAGGCAAAACCCCCGGCCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACATCCCACACGGTAACGATgacatcgaaatc                                                                                                <  2:136003/230‑1 (MQ=255)                                                                                                                                                 aaaGCAACAAGCCTGCGTGCCGGGAAATTGATGGACCTGCCCCCCCTGAATGACAACGACACCGCCGAGAAAGTAATTGCT‑CTGTTTCACCAGCCCAAACCCCCGGCCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGc                                                                                        <  2:15109/251‑1 (MQ=255)                                                                                                                                                        gcagcCTGCGTGCACTGAAATTGCTGGACCTGACCCCCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTCTCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCACAACACTGAAAGCGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGCCATCGCCATCGCGCTTGCAGAAAcc                                                                                 >  1:12427/1‑251 (MQ=255)                                                                                                                                                              tGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGCACCGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCg                                                                           >  1:226463/1‑251 (MQ=255)                                                                                                                                                                    cACTGAAAAGGATGGACCGGCCAACCCGGAATGCCGACGACACCGACTAGAAAGTAATTGCT‑CGGTGCCATCAGGCCAACTCCCGGGGCGCCAAAACCGCCGCTATCTGTATCTATCCGCGCGTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACAGCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                                                                     <  2:12427/251‑1 (MQ=255)                                                                                                                                                                                  gACCTGACCCCCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAACCCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                                                                                   <  1:223821/223‑1 (MQ=255)                                                                                                                                                                                  gACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGaaa                                                                                   >  2:223821/1‑223 (MQ=255)                                                                                                                                                                                  gACCTGACAACCCTGAATGACGACGACACCGAAGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGCCGGCAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGACAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAAACAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGccgg                                                       >  2:163162/1‑248 (MQ=255)                                                                                                                                                                                       gACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCT‑CTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATACGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCcgatga                                                    >  1:41477/1‑249 (MQ=255)                                                                                                                                                                                              cgaaaTGACGACAACACCGACGAAAAAGTAATTGCT‑CTGGGTCATCAGGCAAAACCCCCGGTCGGCAAAACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGt                                           <  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)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            CATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCC‑CTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGG  >  NC_000913/4617191‑4617671

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