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%	157.8 / NA	49	*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 (25/24);  total (25/24)

CTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGC  >  NC_000913/3956743‑3957003                                                                                                                                     |                                                                                                                                 cTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGaaaa                                                                                                                            <  2:1119766/139‑1 (MQ=255)   aCGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAAt                                                                                                                          >  2:778296/1‑139 (MQ=255)        cgcgTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGaaaaa                                                                                                                           >  1:1028081/1‑133 (MQ=255)        cgcgTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGaaaaa                                                                                                                           <  2:1028081/133‑1 (MQ=255)          cgTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGAt                                                                                                                   >  1:173716/1‑139 (MQ=255)          cgTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGAt                                                                                                                   >  1:778114/1‑139 (MQ=255)           gTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATa                                                                                                                  >  1:1185288/1‑139 (MQ=255)            tACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGcc                                                                                                                       >  2:420100/1‑133 (MQ=255)            tACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGcc                                                                                                                       <  1:420100/133‑1 (MQ=255)              cTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAgcg                                                                                                               <  2:1185288/139‑1 (MQ=255)              cTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAgcg                                                                                                               <  2:685437/139‑1 (MQ=255)               tGGCGTCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAGCCCTGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAgcgc                                                                                                              <  2:629405/139‑1 (MQ=255)                   ggcgTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGcgcttcg                                                                                                          <  1:778296/139‑1 (MQ=255)                      gTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAgcg                                                                                                               >  2:199616/1‑131 (MQ=255)                      gTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAgcg                                                                                                               <  1:199616/131‑1 (MQ=255)                            aCTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATAAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAgcg                                                                                                               >  1:561675/1‑125 (MQ=255)                            aCTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATAAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAgcg                                                                                                               <  2:561675/125‑1 (MQ=255)                                          aaCCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGAc                                                                                   >  2:823578/1‑139 (MQ=255)                                                  ttCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGcgcttc                                                                                                           >  2:711404/1‑107 (MQ=255)                                                  ttCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGcgcttc                                                                                                           <  1:711404/107‑1 (MQ=255)                                                         cccGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGt                                                                    <  2:1213424/139‑1 (MQ=255)                                                             gCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGcgcttcg                                                                                                          >  1:965616/1‑97 (MQ=255)                                                             gCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGcgcttcg                                                                                                          <  2:965616/97‑1 (MQ=255)                                                                  aTGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCaaaa                                                           <  2:1173349/139‑1 (MQ=255)                                                                       cTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTATGGCGGGTTAAGGAAAAATGCCTTATAGCGCTTCGCTTAGCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTc                                                      <  1:823578/139‑1 (MQ=255)                                                                           gCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGa                                                  >  1:805867/1‑139 (MQ=255)                                                                           gCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGa                                                  >  2:743256/1‑139 (MQ=255)                                                                                    gCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCaa                                         >  1:694805/1‑139 (MQ=255)                                                                                    gCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCaa                                         >  2:674430/1‑139 (MQ=255)                                                                                     ccACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTTATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCa                                                              >  2:323725/1‑117 (MQ=255)                                                                                     ccACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCa                                                              >  2:998873/1‑117 (MQ=255)                                                                                     ccACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCa                                                              <  1:998873/117‑1 (MQ=255)                                                                                     ccACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCa                                                              <  1:323725/117‑1 (MQ=255)                                                                                      cACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTTGGCAAAATCTTCCAGAATGCCTCAAtt                                       <  2:778114/139‑1 (MQ=255)                                                                                      cACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAAtt                                       <  2:805867/139‑1 (MQ=255)                                                                                       acgacgAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATGTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTa                                      >  1:811958/1‑139 (MQ=255)                                                                                          acgaAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCg                                   <  1:674430/139‑1 (MQ=255)                                                                                               aCCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCa                              >  2:1133170/1‑139 (MQ=255)                                                                                                     aaCCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAAtt                                       <  1:19877/124‑1 (MQ=255)                                                                                                     aaCCCGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAAtt                                       >  2:19877/1‑124 (MQ=255)                                                                                                       cccGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTa                                                                                                      >  1:739259/1‑59 (MQ=255)                                                                                                       cccGGCGTTCAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTa                                                                                                      <  2:739259/59‑1 (MQ=255)                                                                                                                cAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTg             <  2:694805/139‑1 (MQ=255)                                                                                                                cAGGTTCTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTg             <  2:811958/139‑1 (MQ=255)                                                                                                                      cTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGcgcttcg                                                                                                          >  1:969819/1‑40 (MQ=255)                                                                                                                      cTTTTTGGCGGGTTAAGGAAAAATGCCTGATAGcgcttcg                                                                                                          <  2:969819/40‑1 (MQ=255)                                                                                                                         tttGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCacacac    >  1:918513/1‑139 (MQ=255)                                                                                                                         tttGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCacacac    >  2:1170327/1‑139 (MQ=255)                                                                                                                           tGGCGGGTTAAGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGc  >  1:258503/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                 CTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGC  >  NC_000913/3956743‑3957003

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