Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	310	(T)9→8	intergenic (+55/‑27)	thrL → / → thrA	thr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	302	0	T	.	100.0%	219.1 / NA	52	intergenic (+47/‑35)	thrL/thrA	thr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (31/21);  total (31/21)

AGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGG  >  NC_000913/176‑436                                                                                                                               |                                                                                                                                       agTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGGCAGTGCGGGC‑TTTTTTTTCGAcc                                                                                                                           <  1:381812/139‑1 (MQ=255)   tACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCaa                                                                                                                         >  1:899320/1‑139 (MQ=255)          caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAAc                                                                                                                  >  2:734896/1‑139 (MQ=255)             ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAAc                                                                                                                  >  2:395455/1‑136 (MQ=255)             ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAAc                                                                                                                  <  1:395455/136‑1 (MQ=255)             ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAg                                                                                                               >  1:1244602/1‑139 (MQ=255)             ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAg                                                                                                               >  1:517046/1‑139 (MQ=255)             ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAg                                                                                                               >  2:1070810/1‑139 (MQ=255)             ccATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAg                                                                                                               >  2:277173/1‑139 (MQ=255)               aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGt                                                                                                             >  2:320142/1‑139 (MQ=255)                          tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAAc                                                                                                                  <  2:774570/123‑1 (MQ=255)                          tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAAc                                                                                                                  >  1:774570/1‑123 (MQ=255)                                       accaccaccATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGtt                                                                                     <  1:734896/139‑1 (MQ=255)                                        ccaccaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                                                                    <  2:584475/139‑1 (MQ=255)                                          accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                  >  1:1321299/1‑139 (MQ=255)                                             accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcgg                                                                               >  1:346662/1‑139 (MQ=255)                                             accatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcgg                                                                               >  2:1295259/1‑139 (MQ=255)                                                   accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                                                                    >  2:466058/1‑128 (MQ=255)                                                   accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                                                                    <  1:466058/128‑1 (MQ=255)                                                      atTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGt                                                                      <  1:320142/139‑1 (MQ=255)                                                        tACCACAGGTAACGGTGCGGGCTGACGCGTACAGGGAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                                                                    >  2:970563/1‑139 (MQ=255)                                                               ggTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                             >  1:698747/1‑139 (MQ=255)                                                               ggTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                             >  2:845662/1‑139 (MQ=255)                                                                   aCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa                                                         >  2:706907/1‑139 (MQ=255)                                                                     ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGtt                                                                                     <  1:667166/109‑1 (MQ=255)                                                                     ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGtt                                                                                     >  2:667166/1‑109 (MQ=255)                                                                         cGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAAc                                                                                                                  <  2:464863/76‑1 (MQ=255)                                                                         cGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAAc                                                                                                                  >  1:464863/1‑76 (MQ=255)                                                                         cGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa                                                         >  1:451071/1‑133 (MQ=255)                                                                         cGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCAACTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa                                                         <  2:451071/133‑1 (MQ=255)                                                                          gggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                  >  1:556116/1‑139 (MQ=255)                                                                           ggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCCCCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCt                                                 >  1:1384305/1‑139 (MQ=255)                                                                                aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgtg                                            >  2:1213647/1‑139 (MQ=255)                                                                                      aCAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCg                                      <  1:277173/139‑1 (MQ=255)                                                                                              cacaGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                                                                    >  2:842703/1‑85 (MQ=255)                                                                                              cacaGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTc                                                                                    <  1:842703/85‑1 (MQ=255)                                                                                               acaGAAAAAAGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                             <  1:845662/139‑1 (MQ=255)                                                                                                 aaaaaaaaGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGaa                           >  1:1027085/3‑139 (MQ=255)                                                                                                   aaaaaaGCCCGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAg                         <  1:136646/139‑1 (MQ=255)                                                                                                           ccGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGACAAATGc                                                             <  1:736239/95‑1 (MQ=255)                                                                                                           ccGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGACAAATGc                                                             >  2:736239/1‑95 (MQ=255)                                                                                                            cGCACCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcag                <  2:608905/139‑1 (MQ=255)                                                                                                              cACCTGACAGTGC‑GGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGCACGTTTTCTGCGTGTTGc                                        >  1:429155/1‑112 (MQ=255)                                                                                                              cACCTGACAGTGC‑GGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGc                                        <  2:429155/112‑1 (MQ=255)                                                                                                                   gACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg         <  1:1295259/139‑1 (MQ=255)                                                                                                                      aGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                   >  2:384899/1‑126 (MQ=255)                                                                                                                      aGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGc                   <  1:384899/126‑1 (MQ=255)                                                                                                                        tGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                  >  1:1140304/1‑61 (MQ=255)                                                                                                                        tGCGGGC‑TGTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcgg                                                                                  <  2:1140304/61‑1 (MQ=255)                                                                                                                         gCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                        >  2:914929/1‑118 (MQ=255)                                                                                                                         gCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                        <  1:914929/118‑1 (MQ=255)                                                                                                                         gCGGGC‑TTTTTTTTCGACCAAAGGTAACG‑GGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTgg  >  2:712625/1‑139 (MQ=255)                                                                                                                               |                                                                                                                                       AGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGG  >  NC_000913/176‑436

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