Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,956,875	G→A	*515* (TAG→TAA)	ilvA →	l‑threonine dehydratase, biosynthetic, also known as threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,956,875	0	G	A	100.0%	95.0 / NA	30	*515* (TAG→TAA)	ilvA	l‑threonine dehydratase, biosynthetic, also known as threonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/17);  total (13/17)

ACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACG  >  NC_000913/3956742‑3957002                                                                                                                                      |                                                                                                                                aCTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGaaa                                                                                                                            >  2:253268/1‑139 (MQ=255)  cTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGaaaa                                                                                                                           <  1:554745/139‑1 (MQ=255)          gcgTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGa                                                                                                                   <  1:797365/139‑1 (MQ=255)           cgTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGAt                                                                                                                  >  1:1037671/1‑139 (MQ=255)             tACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAg                                                                                                                >  1:776050/1‑139 (MQ=255)                        ttCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTa                                                                                                     <  2:247768/139‑1 (MQ=255)                         tCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTAt                                                                                                    <  1:880107/139‑1 (MQ=255)                          cGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGcgctt                                                                                                           >  2:550559/1‑131 (MQ=255)                          cGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGcgctt                                                                                                           <  1:550559/131‑1 (MQ=255)                                       cATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCcgcg                                                                                      <  1:950729/139‑1 (MQ=255)                                         tGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCcgcgcg                                                                                    >  2:649322/1‑139 (MQ=255)                                           aaCCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGAc                                                                                  <  2:692718/139‑1 (MQ=255)                                                gATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGt                                                                             >  2:911998/1‑139 (MQ=255)                                                           ccGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGtt                                                                  >  1:259993/1‑139 (MQ=255)                                                                  aaTGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCaaa                                                           <  2:568625/139‑1 (MQ=255)                                                                       gCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCtt                                                      >  2:452661/1‑139 (MQ=255)                                                                            gCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGa                                                 <  2:452946/139‑1 (MQ=255)                                                                                     gCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGACAAATGCCTGAAAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGAGGTTCGGCAAAATCTTCCAGAATGCCTCaa                                        >  2:922591/1‑139 (MQ=255)                                                                                       cACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAAtt                                      <  1:911998/139‑1 (MQ=255)                                                                                        acgacgAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTa                                     <  1:649322/139‑1 (MQ=255)                                                                                        acgacgAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTa                                     <  1:452661/139‑1 (MQ=255)                                                                                        acgacgAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTa                                     <  1:233853/139‑1 (MQ=255)                                                                                                    aTAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTa                         <  1:253268/139‑1 (MQ=255)                                                                                                      aaCCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGc                       >  1:369093/1‑139 (MQ=255)                                                                                                             cGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCtttt                >  1:910517/1‑139 (MQ=255)                                                                                                                 cAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTg            <  2:369093/139‑1 (MQ=255)                                                                                                                             ggCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGGTTTTTCTGCGCACACACg  <  2:586554/137‑1 (MQ=255)                                                                                                                             ggCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACg  >  1:586554/1‑137 (MQ=255)                                                                                                                              gCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGAc                                                                                  <  1:620583/56‑1 (MQ=255)                                                                                                                              gCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGAc                                                                                  >  2:620583/1‑56 (MQ=255)                                                                                                                                      |                                                                                                                                ACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACG  >  NC_000913/3956742‑3957002

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