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%	99.6 / NA	35	T38T (ACT→ACC)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/9);  total (26/9)

TTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGT  >  NC_000913/4617193‑4617672                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              ttACAGTGATGCAAACTTGTAAGTAGATTTCCTTCATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCg                                                                                                                                                                                                                                       >  1:106400/1‑251 (MQ=255)                              tCCTTAATTGTGATGTCTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGAAAGGTGAA‑TTTTATTGGAGACAAGACAGGAGAATGCAATGACTGATCTGAAAGCAAGCAGCCTGCGTGAACTGAATTTGATGGACCTGACCA‑CCATG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCACAACCACGG‑ACGGCAATACCGCCGCTATCTGTATCTAGCCt                                                                                                                                                                                                          >  2:74447/1‑251 (MQ=255)                               ccTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTc                                                                                                                                                                                                         >  1:263504/1‑251 (MQ=255)                                  tAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCt                                                                                                                                                                                                      >  1:358988/1‑251 (MQ=255)                                    ggctcttccgtcttCGAGGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGAAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCC‑CCCTG‑AATGACGACGACACCGACGAAAAAGTAATTGCTCTGTGTCATCAGGCCAAACCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCttt                                                                                                                                                                                                    <  1:83516/238‑1 (MQ=255)                                      tgtgATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGACAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTACTCTGTGTCATCAGGCCAAAACCACGG‑TCGGCAATACCGCAGCTATCGGTATCAATCCTCGATTTa                                                                                                                                                                                                   >  2:16745/1‑250 (MQ=255)                                                  cGAAGTGTGTTGCGGAGTAGATGTTAGCATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCttt                                                                                                                                                                                                    >  2:83516/1‑237 (MQ=255)                                                             gCGGATTAGATGTTAAAATAATAACAAACTCG‑AAAGTGAAATTTTATGGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCC‑CCCGG‑AATGACGACGACCCCGACGAGAAAGTAATTGCTCTGTGACATCAGGCAAAACCCCCGG‑CCGGCAATTCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGa                                                                                                                                                                           <  1:207878/251‑1 (MQ=255)                                                               ggcgTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAAt                                                                                                                                                                                                                                  >  1:57535/4‑194 (MQ=255)                                                               ggAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCC‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAAt                                                                                                                                                                                                                                  <  2:57536/194‑1 (MQ=255)                                                               ggAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAAt                                                                                                                                                                                                                                  >  1:57536/1‑194 (MQ=255)                                                                  gTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATAACTGATCTGAAAGCACGCAGCCTCCGTGCACTGAAATTGCTAGACCT‑ACCAACCCTG‑AATGCCGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑CCGGCAATACCGCCGCCATCTGAATCTATCCTCGCTTTATCCCGATTGCTCCCAAACCCCTGaaaagg                                                                                                                                                                      >  1:352849/1‑248 (MQ=255)                                                                   tAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGc                                                                                                                                                                     >  1:220891/1‑251 (MQ=255)                                                                       tGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGATTGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGCCGACACCGACGAGAAAGTAATTGCCCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGAAggg                                                                                                                                                                 >  1:326711/1‑251 (MQ=255)                                                                            gAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCAccc                                                                                                                                                            >  1:322283/1‑251 (MQ=255)                                                                              aTACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCC‑CCCTG‑AATGACGACGACACCGACGAGATAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCACCGCTATCTGTATCTATCCTCGCTTTATCACGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCg                                                                                                                                                          >  1:372026/1‑251 (MQ=255)                                                                              aTACTAACAAACTCGCAAGGTGAA‑TTTTATTGGAGACAAGCCAGGAGAATGAAATGCCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATAGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCTACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAAACCCGG‑TCGGCAGTACCGCCGCTATCTGTATCTATCCTAGCTTGATCCCGATTGCTCGAAAAACACTGAAAGGGCAGGGCAccc                                                                                                                                                            >  2:211004/1‑249 (MQ=255)                                                                              aTACTAACAAAATCGCAAGGTGAA‑TTTTATTGGAGACAAGCCAGGAGAAGGAAATGACTGATCTGAAAGAAAGCAGCCTGCGAGAAATGAAATTGAAGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGACATCAGGCCAAAACCCAGG‑TCGGCAATAACGCCGCTATCTGTATCAATCCTAGCTATATCACGATTGCTCGCAACATACTGAAAGAGCAGGGAAccc                                                                                                                                                            >  1:297125/1‑249 (MQ=255)                                                                               tACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCACAACCCCGG‑ACGGCAATACCGCCGCTATATGTCTCTATCCTCGCTTTATCCCGGTTGCTCGCAAAACACTGAAAGAGCCGGGCAcccc                                                                                                                                                           >  1:262782/1‑249 (MQ=255)                                                                                   aaCAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGtatctat                                                                                                                                                                                                             >  1:119540/1‑195 (MQ=255)                                                                                      aaaCTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCg                                                                                                                                                  >  1:14706/1‑251 (MQ=255)                                                                                                       ttttATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAAc                                                                                                                                   >  1:65004/1‑250 (MQ=255)                                                                                                             gggcaaCAAGCCAGGAAAATGAAATGACTGATCTGAAAGAAAGCAGCCTGCGTGCACTAAAATTGATGGACCTGCCCA‑CCCGG‑AATGACGACGACACCGACGAAAAAGTAATTGCTCTGTGTCATCAGGCAAAACCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTc                                                                                                                            <  1:278295/246‑1 (MQ=255)                                                                                                                   aaaGCAAGGAGAATGAATTATCTGATCTAAAAGCAAGAAGCTTGCGGGCACTGAAATTGATGGACCTGACC‑‑CCCTGTAATGACGACGACCCCGACGAGAAAGTAATTGCTCGGTGTCACCAGCCCAAACCCCCGG‑GCGACAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCacac                                                                                                                      <  1:272010/250‑1 (MQ=255)                                                                                                                            agaATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAacgat                                                                                                             >  1:120141/1‑250 (MQ=255)                                                                                                                                tgaaatgaCTGATCTGAAAGCAAGCAGACTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTGCCCACACGGTAacgatgcc                                                                                                          >  1:82832/1‑246 (MQ=255)                                                                                                                                      gaCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGCCCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCaaaa                                                                                                                                                                                 >  1:117924/1‑173 (MQ=255)                                                                                                                                             cTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCgcgc                                                                                             >  1:366178/1‑250 (MQ=255)                                                                                                                                                               cGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGc                                                                          >  2:112977/1‑251 (MQ=255)                                                                                                                                                                   cACTGAAATTGAGGGCCCTGCCCC‑CCCTG‑AAGGACGACGACACCGACGAGAATGTAATTGCTTTGTTTCTCCGGGCCAAAACCCCGG‑TCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                                                                      <  1:221741/251‑1 (MQ=255)                                                                                                                                                                                   ccatgcccccccTG‑AATGACGACGACCCCGACGAAAAAGTAATTGCTCTGTGTCATCGGGCCAAACCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCcgat                                                       <  1:112977/242‑1 (MQ=255)                                                                                                                                                                                                              cccGACGATAAAGTAATTGCGCTGTGACATCAGCCCAAACCCCCGG‑TCGCAAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACGTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCTTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACcg                             <  1:152271/250‑1 (MQ=255)                                                                                                                                                                                                                cgacgaGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACcgcg                           >  1:271164/1‑251 (MQ=255)                                                                                                                                                                                                                   cgaGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCgcgcgc                        >  1:344025/1‑251 (MQ=255)                                                                                                                                                                                                                                          caccagtcAAAACCCCGGGTCGGCAATTCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggt  <  1:262823/244‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              TTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGT  >  NC_000913/4617193‑4617672

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