Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,959,045	T→C	V359A (GTA→GCA)	ilvC →	ketol‑acid reductoisomerase (NADP(+))

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,959,045	0	T	C	100.0%	95.0 / NA	29	V359A (GTA→GCA)	ilvC	ketol‑acid reductoisomerase (NADP(+))
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/12);  total (17/12)

ATGGCGGACTGGGCCAACGATGATAAGAAACTGCTGACCTGGCGTGAAGAGACCGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGTACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCATATTATGAATCACTGCACGAGCTGCCGCTGATTGCCAACACCATCGC  >  NC_000913/3958918‑3959174                                                                                                                                |                                                                                                                                  atgGCGGACTGGGCCAACGATGATAAGAAACTGCTGACCTGGCGTGAAGAGACCGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACtga                                                                                                                              >  1:34216/1‑133 (MQ=255)     cGGACTGGGCCAACGATGATAAGAAACTGCTGACCTGGCGTGAAGAGACCGGCAAAACCGCATTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTgatgat                                                                                                                          >  1:31791/1‑133 (MQ=255)       gACTGGGCCAACGATGATAAGAAACTGCTGACCTGGCGTGAAGAGACCGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAAGAGTACTTCGATAAAGGCGCACTGATGATTg                                                                                                                        >  1:513031/1‑133 (MQ=255)                    atgatAAGAAACTGCTGACCTGGCGTGAAGAGACCGGCAAAACCGCATTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGc                                                                                                           <  2:31791/133‑1 (MQ=255)                               ctgctgACCTGGCGTGAAGAGACCGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGGACTTTGATAAAGGCGCACTGATGATTGCGATGGGGAAAGCGGGCGTTGAAc                                                                                                >  1:290637/1‑133 (MQ=255)                                  ctgACCTGGCGTGAAGAGACCGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTgg                                                                                             >  1:313578/1‑133 (MQ=255)                                  ctgACCTGGCGTGAAGAGACCGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTgg                                                                                             <  1:500459/133‑1 (MQ=255)                                          gCGTGAAGAGACCGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGGTGGTGAAAGCGGGCGTAGAACTGGCGTTCGaa                                                                                     >  2:205157/1‑133 (MQ=255)                                          gCGTGAAGAGACCGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGaa                                                                                     >  1:442094/1‑133 (MQ=255)                                                    aCCGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCg                                                                           >  2:77514/1‑133 (MQ=255)                                                    aCCGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCg                                                                           >  2:41345/1‑133 (MQ=255)                                                    aCCGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACCTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCg                                                                           >  2:454269/1‑133 (MQ=255)                                                    aCCGACAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCg                                                                           >  2:281595/1‑133 (MQ=255)                                                      cGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGAt                                                                         <  1:53621/133‑1 (MQ=255)                                                           aaaCCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCgg                                                                    <  1:469072/133‑1 (MQ=255)                                                                gcgTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGcatca                                                               <  2:313578/133‑1 (MQ=255)                                                                       aaaCCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAaga                                                        <  1:77514/133‑1 (MQ=255)                                                                          ccGCGCCGTAGTATGAAGTCAATATCGGCGCGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTc                                                     <  1:161672/133‑1 (MQ=255)                                                                              gccgcAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCa                                                 <  2:395166/133‑1 (MQ=255)                                                                                  cAGTATGAAGTCAAACTCGGCGAGCAGGGGTACTTCGATAATGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCAtatt                                             <  2:214975/133‑1 (MQ=255)                                                                                        gAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCATATTATGAAt                                       <  2:254227/133‑1 (MQ=255)                                                                                             cAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCATATTATGAATCACTg                                  >  1:410453/1‑133 (MQ=255)                                                                                             cAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCATATTATGAATCACTg                                  >  1:387207/1‑133 (MQ=255)                                                                                              aaaaTCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCATATTATGAATCACTGc                                 <  1:2466/133‑1 (MQ=255)                                                                                               aaaTCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCATATTATGAATCACTGCa                                >  2:61881/1‑133 (MQ=255)                                                                                               aaaTCGGCGAGCAGGAGTACTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGCCTGGATATTATGAATCACTGCa                                >  1:172176/1‑133 (MQ=255)                                                                                                            ggAGTCCTTCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCATATTATGAATCACTGCACGAGCTGCCGCTg                   <  2:442094/133‑1 (MQ=255)                                                                                                                    tCGATAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCAGCATATTAGGAAGCCCAGCACGAGCGGCCGAGGATTGccca           >  2:248179/1‑131 (MQ=255)                                                                                                                       aTAAAGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCATATTATGAATCACTGCACGAGCTGCCGCTGATTGCCAacac        >  2:460460/1‑133 (MQ=255)                                                                                                                         aaaGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCATATTATGAATCACTGCACGAGCTGCCGCTGATTGCCATCACCa      <  2:290637/133‑1 (MQ=255)                                                                                                                         aaaGGCGCACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCATATTATGAATCACTGCACGAGCTGCCGCTGATTGCCAACACCa      <  2:387207/133‑1 (MQ=255)                                                                                                                             gcgcACTGCTGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCATATTATGAATCACTGCACGAGCTGCCGCTGATTGCCAACACCATCGc  <  1:205157/133‑1 (MQ=255)                                                                                                                                |                                                                                                                                  ATGGCGGACTGGGCCAACGATGATAAGAAACTGCTGACCTGGCGTGAAGAGACCGGCAAAACCGCGTTTGAAACCGCGCCGCAGTATGAAGGCAAAATCGGCGAGCAGGAGTACTTCGATAAAGGCGTACTGATGATTGCGATGGTGAAAGCGGGCGTTGAACTGGCGTTCGAAACCATGGTCGATTCCGGCATCATTGAAGAGTCTGCATATTATGAATCACTGCACGAGCTGCCGCTGATTGCCAACACCATCGC  >  NC_000913/3958918‑3959174

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