Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,393,446	G→T	E74* (GAA→TAA)	mppA →	murein tripeptide ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,393,446	0	G	T	100.0%	82.5 / NA	25	E74* (GAA→TAA)	mppA	murein tripeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (10/15);  total (10/15)

GATCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGC  >  NC_000913/1393371‑1393523                                                                            |                                                                              gATCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTaaaa                                                                           <  2:247672/80‑1 (MQ=255)     ccGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAggg                                                                       <  1:247673/80‑1 (MQ=255)       gCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGgag                                                                    <  2:247674/81‑1 (MQ=255)             gCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTc                                                              <  2:247676/81‑1 (MQ=255)                gTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTcccc                                                           <  1:247677/81‑1 (MQ=255)                  gggCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCgg                                                         >  2:247715/1‑81 (MQ=255)                   ggCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGc                                                        <  1:247678/81‑1 (MQ=255)                     ccTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGt                                                      <  2:247679/81‑1 (MQ=255)                         ccAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGc                                                   <  1:247680/80‑1 (MQ=255)                         ccAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGc                                                   <  2:247681/80‑1 (MQ=255)                          cAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGa                                                 >  1:247728/1‑81 (MQ=255)                                 tCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTg                                          >  2:247714/1‑81 (MQ=255)                                   aGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTgg                                         <  2:247682/80‑1 (MQ=255)                                     gTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGaaa                                      <  1:247683/81‑1 (MQ=255)                                     gTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGaaa                                      <  2:247684/81‑1 (MQ=255)                                            gcgATCTGTTTGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATg                               <  2:247685/81‑1 (MQ=255)                                                    tttGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGt                       <  1:247686/81‑1 (MQ=255)                                                     ttGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTa                      >  1:247707/1‑81 (MQ=255)                                                     ttGAAGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTa                      >  1:247711/1‑81 (MQ=255)                                                        aaGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCt                   <  2:247687/81‑1 (MQ=255)                                                         aGGTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTg                  >  1:247712/1‑81 (MQ=255)                                                           gTCTGGTGAATCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGa                >  2:247713/1‑81 (MQ=255)                                                                    aTCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTAccc       >  2:247710/1‑81 (MQ=255)                                                                    aTCAGAACTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTAccc       >  2:247729/1‑81 (MQ=255)                                                                         aaCTAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTgcgc  >  1:247724/1‑81 (MQ=255)                                                                            |                                                                              GATCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGC  >  NC_000913/1393371‑1393523

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