Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,232,992	G→T	P81Q (CCG→CAG)	lysC ←	aspartate kinase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,232,992	0	G	T	100.0%	142.6 / NA	41	P81Q (CCG→CAG)	lysC	aspartate kinase III
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/25);  total (16/25)

AACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCGGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACCAGCAGATTAGTGATACCAGCAGAAGCCGAGAGGAC  >  NC_000913/4232857‑4233125                                                                                                                                        |                                                                                                                                      aacaGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTa                                                                                                                                   >  1:261761/1‑140 (MQ=255)  acaGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTaa                                                                                                                                  >  1:164027/1‑140 (MQ=255)    agcagGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACg                                                                                                                                >  2:163934/1‑140 (MQ=255)       agGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAg                                                                                                                             <  2:24760/140‑1 (MQ=255)       agGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAg                                                                                                                             >  1:250213/1‑140 (MQ=255)       agGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAg                                                                                                                             <  2:443401/140‑1 (MQ=255)             gACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTc                                                                                                                       <  1:42706/140‑1 (MQ=255)               cATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCa                                                                                                                     <  2:99669/140‑1 (MQ=255)               cATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCa                                                                                                                     <  1:212425/140‑1 (MQ=255)                    gCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATg                                                                                                                >  2:363168/1‑140 (MQ=255)                     cTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATgg                                                                                                               >  2:430365/1‑140 (MQ=255)                        gCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCaa                                                                                                            <  2:271131/140‑1 (MQ=255)                         ccGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCaaa                                                                                                           <  1:455465/140‑1 (MQ=255)                               cTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGAt                                                                                                     <  2:261761/140‑1 (MQ=255)                                       cTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGa                                                                                             <  2:164027/140‑1 (MQ=255)                                            cTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCg                                                                                        >  1:205462/1‑140 (MQ=255)                                              gTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTc                                                                                      <  1:375733/140‑1 (MQ=255)                                                   cgcCGGAGACGTTGCCAGCGCCGCCGCTTCTGCAAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGtt                                                                                 >  2:206576/1‑140 (MQ=255)                                                     ccGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGtttt                                                                               <  2:236958/140‑1 (MQ=255)                                                            cGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAAtc                                                                        <  1:328159/140‑1 (MQ=255)                                                               tGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAAtcgct                                                                     >  2:257256/1‑140 (MQ=255)                                                               tGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAAtcgct                                                                     <  1:449292/140‑1 (MQ=255)                                                                   aGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGcc                                                                 <  1:163934/140‑1 (MQ=255)                                                                           cgcTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCa                                                         <  2:296070/140‑1 (MQ=255)                                                                                       aaCAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTa                                             >  2:426984/1‑140 (MQ=255)                                                                                             aaTGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGa                                       >  1:67099/1‑140 (MQ=255)                                                                                              aTGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGAc                                      <  2:183405/140‑1 (MQ=255)                                                                                              aTGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGAc                                      <  2:98143/140‑1 (MQ=255)                                                                                                 ttCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACcag                                   >  1:90771/1‑140 (MQ=255)                                                                                                 ttCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACcag                                   <  1:257256/140‑1 (MQ=255)                                                                                                    tcCAGCAGACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACcagcag                                <  1:430365/140‑1 (MQ=255)                                                                                                      cagcagACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACCAGCAGAt                              <  2:205462/140‑1 (MQ=255)                                                                                                      cagcagACGTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGCTCCAGTCCTTCAGCTAAAGCGACCAGCAGAt                              >  1:102860/1‑140 (MQ=255)                                                                                                              gTTCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACCAGCAGATTAGTGATa                      <  2:328682/140‑1 (MQ=255)                                                                                                                tCAATCTCTTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACCAGCAGATTAGTGATAcc                    <  1:363168/140‑1 (MQ=255)                                                                                                                     ctctTCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACCAGCAGATTAGTGATACcagcag               <  2:410021/140‑1 (MQ=255)                                                                                                                        ttCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACCAGCAGATTAGTGATACCAGCAGAAg            >  2:361400/1‑140 (MQ=255)                                                                                                                        ttCACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACCAGCAGATTAGTGATACCAACAGAAg            >  2:229875/1‑140 (MQ=255)                                                                                                                          cACGGATAACGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACCAGCAGATTAGTGATACCAGCAGAAGcc          <  1:361400/140‑1 (MQ=255)                                                                                                                                 aaCGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACCAGCAGATTAGTGATACCAGCAGAAGCCGAGAGGa   <  2:67099/140‑1 (MQ=255)                                                                                                                                  aCGTTCTGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACCAGCAGATTAGTGATACCAGCAGAAGCCGAGAGGAc  >  2:341652/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                      AACAGCAGGGTCGACATCAGCTCGCCGTGGCTGACCAGCTCATCTGTCAGCGCCGGAGACGTTGCCAGCGCCGCCGCTTCTGCCAGAACAGTAATGTTCTCCAGCAGACGTTCAATCTCTTCACGGATAACGTTCGGGTAACGCAGACGTTCCAGAATGGCAAACTGGATGTTGCGGATAGCGTCGAGTTTTTCGAATCGCTCGCCAGGTTCCAGTCCTTCAGCTAAAGCGACCAGCAGATTAGTGATACCAGCAGAAGCCGAGAGGAC  >  NC_000913/4232857‑4233125

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