Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,780	T→C	T349T (ACT→ACC)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,780	0	T	C	100.0%	91.7 / NA	30	T349T (ACT→ACC)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/17);  total (13/17)

AACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCT  >  NC_000913/4533‑5025                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         aaCGTTTTATTGCTGCGACCAACGTGAACGATACCGTGTCACGTTTCCTGCACGGCGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTGATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTTTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCg                                                                                                                                                                                                                                                       <  1:244211/251‑1 (MQ=255)         aTTGCTGCGACCAACGTGAACGATACCGTGCCCCGTTTCCTGCACGACGGTCAGTGGTCCCCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCa                                                                                                                                                                                                                                               <  1:312730/251‑1 (MQ=255)                           ttccgagcggtgcccctttcccagCACGACGGGCAGTGTTCACCCAAAGCGCCTCAGGCGACGTTACCCAACGCGATGGACGGGAGTCAGCCAAACAACTGGCGGCGTGTGGAAGATTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATc                                                                                                                                                                                                                              <  2:394919/227‑1 (MQ=255)                                    gTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCACAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATc                                                                                                                                                                                                                              >  1:394919/1‑241 (MQ=255)                                             ttCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑Gtt                                                                                                                                                                                                          >  1:350103/1‑251 (MQ=255)                                                 tGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAAt                                                                                                                                                                                                      >  1:437112/1‑251 (MQ=255)                                                                    cccaaaaGCGACTCAGGGAACGTTACCCAAGGGGTGGGACGTGAGTCAGCCGAACAACTGGCCGGGTGGGAAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGt                                                                                                                                                                                   <  2:484471/247‑1 (MQ=255)                                                                      ccAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGCACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACAACTCGGAGCAGCACGCTGCCGTAGCGTATCGTGCGCTGCG‑GAACAAGTTGAATCCAGGCGAATATGGCTTGTTc                                                                                                                                                                                 >  1:381167/1‑251 (MQ=255)                                                                                                aaCGCGATGGACGTGATTCAGCCGCACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGCCAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCCGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGGAAtt                                                                                                                                                       >  1:450085/1‑251 (MQ=255)                                                                                                               aGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGAATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCCCGCCTCCGGCGAAATTT‑‑AAAGAGAGCggggg                                                                                                                                      >  2:360378/1‑247 (MQ=255)                                                                                                                gTCAGCCGAACAACGTGCCGGGGGTGGAAGGGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGaa                                                                                                                                     <  2:250062/251‑1 (MQ=255)                                                                                                                 tCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAACTCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAAACAAGCAACAGACAATGAGGGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCAGCCGTAGCGAATCGGGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATGTGGATTGTTCCTCAGCACCTCACATCCGGCGACATTTCAAAagagagagaa                                                                                                                                        >  2:78303/1‑245 (MQ=255)                                                                                                                   cgccgaccactgtgCCGCGTGTGTAAGAGTGGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATGT‑‑AAAGAGAGCGTGGAAGCg                                                                                                                                  <  2:441408/238‑1 (MQ=255)                                                                                                                         acaccTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATtctc                                                                                                                            <  1:469407/247‑1 (MQ=255)                                                                                                                         acaatgggccgggggtgGAAGTGTGGTTCCGCCGCAAAAACTGTCTACTGAAAGAGCTGGGTTATGCAGCCGTTGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTATACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATtctc                                                                                                                            <  2:484006/237‑1 (MQ=255)                                                                                                                             cTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGt                                                                                                                                                                                                                            >  2:200152/1‑154 (MQ=255)                                                                                                                             cTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGt                                                                                                                                                                                                                            <  1:200152/154‑1 (MQ=255)                                                                                                                                gCCGCGTGGGAAAGGGTGGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCGGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATTCTCGGTGaaa                                                                                                                     <  2:6450/251‑1 (MQ=255)                                                                                                                                 ccGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATTCTCGGTGAAAc                                                                                                                    <  1:182889/251‑1 (MQ=255)                                                                                                                                                    tCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAag                                                                                                   >  1:83681/1‑249 (MQ=255)                                                                                                                                                                   ggaaaCTGAAAGAGCGTGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGGTAAAAGAACTGGGCGACACCTCGGGGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCGTGTTCCTCGGCACCGCGCATCCGGCGAAATTG‑‑AAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGc                                                                                   <  2:255114/247‑1 (MQ=255)                                                                                                                                                                             agagCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGGGCTGATTTACCCt                                                                        >  1:267811/1‑251 (MQ=255)                                                                                                                                                                                                  ggATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGccc                                                    <  1:75739/250‑1 (MQ=255)                                                                                                                                                                                                   gatgatGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACCCCTCGGAGCCGCACGCTGCCGTAGCGCATAGTGCGCTGCGTGCCCA‑GTTGAATCCAGCCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGACGAAATTT‑‑AAAGAAAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGACCGTGCTGCTTTAACCCTGCTTTCACATAATCTGCccgc                                                  >  1:484427/1‑251 (MQ=255)                                                                                                                                                                                                                             gCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCCCGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtg                         <  1:344557/250‑1 (MQ=255)                                                                                                                                                                                                                                ggattaaaaaaaACTGGCCTACTCCTCGGGGCCGCACGCGGCCGTAGCGTTTCGTGCGCTGCGTGACCA‑GTTGATTCCAGGCGAATATGGCTTGTTCCTCGGCCCCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGTTTCTCGGTGAACCGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtgatga                     <  1:367367/245‑1 (MQ=255)                                                                                                                                                                                                                                         aGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatcatc            <  1:150570/251‑1 (MQ=255)                                                                                                                                                                                                                                              ggggCTACCCCTCGGAGCCGCACGCTGCGGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCCCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTCCCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTaa       <  1:492819/250‑1 (MQ=255)                                                                                                                                                                                                                                                  cTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATc   >  1:79790/1‑251 (MQ=255)                                                                                                                                                                                                                                                    acacCTCGGGGCCGGACGATGCAGTAGCGGATCGTGCGCTGAGTGACCC‑GTTGAATCCAGGCGAATATGGCTTGGTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGAGTCTCGGTGAAACATTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTCATTTACCCTTGCTTTCGCATAGTCTGCCCGCCGGTTTTGCTGCGTTGTGGAAATTGCTTAAGAATCATCAGTAAActct  >  2:90595/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         AACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCA‑GTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTT‑‑AAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCT  >  NC_000913/4533‑5025

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