Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	406,492	T→C	F30L (TTT→CTT)	aroL →	shikimate kinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	406,492	0	T	C	100.0%	179.2 / NA	50	F30L (TTT→CTT)	aroL	shikimate kinase II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/23);  total (27/23)

ATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGTTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAACGGCGGCGCTGGAAGCGGTAACTGCGCCATCCAC  >  NC_000913/406358‑406622                                                                                                                                       |                                                                                                                                   aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGt                                                                                                                                >  2:223690/1‑139 (MQ=255) aTGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGt                                                                                                                                >  1:994166/1‑139 (MQ=255)  tGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTc                                                                                                                               >  1:96335/1‑139 (MQ=255)  tGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTc                                                                                                                               >  1:460592/1‑139 (MQ=255)      acCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATa                                                                                                                           >  1:875534/1‑139 (MQ=255)        cGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCg                                                                                                                              <  1:1036252/134‑1 (MQ=255)        cGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCg                                                                                                                              >  2:1036252/1‑134 (MQ=255)        cGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATAcc                                                                                                                         >  2:521313/1‑139 (MQ=255)                         cGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATc                                                                                                        <  2:994166/139‑1 (MQ=255)                                        aaCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGAc                                                                                         <  1:1115802/139‑1 (MQ=255)                                         aCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACg                                                                                        <  2:452922/139‑1 (MQ=255)                                         aCCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACg                                                                                        <  1:223690/139‑1 (MQ=255)                                          cccACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACgg                                                                                       >  2:896135/1‑139 (MQ=255)                                          cccACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACg                                                                                        <  2:779180/138‑1 (MQ=255)                                          cccACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACg                                                                                        >  1:779180/1‑138 (MQ=255)                                                aTGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCgg                                                                                 <  1:521313/139‑1 (MQ=255)                                                    cacaACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGAt                                                                             >  2:121760/1‑139 (MQ=255)                                                       aaCCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGAtcgt                                                                          <  2:496143/139‑1 (MQ=255)                                                                     aTCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGt                                                            >  2:664632/1‑139 (MQ=255)                                                                          gCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCgg                                                                                 >  2:961794/1‑113 (MQ=255)                                                                          gCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCgg                                                                                 <  1:961794/113‑1 (MQ=255)                                                                                gggCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGAttt                                                 <  2:96335/139‑1 (MQ=255)                                                                                      tgGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGcc                                           >  1:492499/1‑139 (MQ=255)                                                                                                acGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGt                                                                                                                                <  1:231262/44‑1 (MQ=255)                                                                                                acGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGt                                                                                                                                >  2:231262/1‑44 (MQ=255)                                                                                                acGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAAt                                                                                                         <  1:216337/67‑1 (MQ=255)                                                                                                acGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAAt                                                                                                         >  2:216337/1‑67 (MQ=255)                                                                                                acGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAAcgg                                 >  2:647739/1‑139 (MQ=255)                                                                                                  ggTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAg                                                                                                    <  2:789025/70‑1 (MQ=255)                                                                                                  ggTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAg                                                                                                    >  1:789025/1‑70 (MQ=255)                                                                                                   gTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAatat                                                                                            <  2:1113082/77‑1 (MQ=255)                                                                                                   gTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAatat                                                                                            >  1:1113082/1‑77 (MQ=255)                                                                                                        aaTGGCCCTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAACGGCGGCGCTg                         <  1:896135/139‑1 (MQ=255)                                                                                                             cccTTGCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAACGGCGGCGCTGGAAGc                    >  1:504645/1‑139 (MQ=255)                                                                                                                  gCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGc                                                                                                            <  2:363643/46‑1 (MQ=255)                                                                                                                  gCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGc                                                                                                            >  1:363643/1‑46 (MQ=255)                                                                                                                  gCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAACGGCGGCGCTGGAAGCGGTaa               >  2:345035/1‑139 (MQ=255)                                                                                                                  gCCGATTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAACGGCGGCGCTGGAAGCGGTaa               <  2:460592/139‑1 (MQ=255)                                                                                                                      aTTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGa                                                                                          <  2:1059579/60‑1 (MQ=255)                                                                                                                      aTTCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGa                                                                                          >  1:1059579/1‑60 (MQ=255)                                                                                                                       ttCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACg                                                                                        <  2:167609/61‑1 (MQ=255)                                                                                                                       ttCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACg                                                                                        >  1:167609/1‑61 (MQ=255)                                                                                                                       ttCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAAcggc                                >  2:376403/1‑117 (MQ=255)                                                                                                                       ttCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAAcggc                                <  1:376403/117‑1 (MQ=255)                                                                                                                       ttCGCTTAACCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAACGGCGGCGCTGGAAGCGGTAACTgcg          <  1:121760/139‑1 (MQ=255)                                                                                                                               aCCGTCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAACGGCGGCGCTGGAAGCGGTAACTGCGCCATCCAc  >  2:230337/1‑139 (MQ=255)                                                                                                                                   tCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAACGGCGGCGCTTGAAg                     <  2:659808/116‑1 (MQ=255)                                                                                                                                   tCGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAACGGCGGCGCTGGAAg                     >  1:659808/1‑116 (MQ=255)                                                                                                                                    cGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGAttt                                                 >  2:395210/1‑87 (MQ=255)                                                                                                                                    cGGCTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGAttt                                                 <  1:395210/87‑1 (MQ=255)                                                                                                                                       |                                                                                                                                   ATGACACCGGCTTTCGCCGCATTGCGACCTATTGGGGAAAACCCACGATGACACAACCTCTTTTTCTGATCGGGCCTCGGGGCTGTGGTAAAACAACGGTCGGAATGGCCCTTGCCGATTCGCTTAACCGTCGGTTTGTCGATACCGATCAGTGGTTGCAATCACAGCTCAATATGACGGTCGCGGAGATCGTCGAAAGGGAAGAGTGGGCGGGATTTCGCGCCAGAGAAACGGCGGCGCTGGAAGCGGTAACTGCGCCATCCAC  >  NC_000913/406358‑406622

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