Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	470,149	G→C	A427P (GCG→CCG)	mdlA →	ABC transporter family protein MdlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	470,149	0	G	C	100.0%	116.3 / NA	33	A427P (GCG→CCG)	mdlA	ABC transporter family protein MdlA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (20/13);  total (20/13)

ACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCC  >  NC_000913/470071‑470217                                                                               |                                                                     aCGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCg                                                                    <  2:103066/81‑1 (MQ=255)   gAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGTCGCCATTCCTTTTTTCTGACACTGTGCCGaa                                                                  >  2:103116/1‑81 (MQ=255)   gAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGaa                                                                  >  2:103109/1‑81 (MQ=255)   gAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGaa                                                                  >  2:103101/1‑81 (MQ=255)       ttACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATaa                                                               >  1:103100/1‑80 (MQ=255)       ttACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAAc                                                              <  1:103069/81‑1 (MQ=255)          cAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATc                                                           <  1:103070/81‑1 (MQ=255)          cAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATc                                                           <  1:103071/81‑1 (MQ=255)            aCTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATcgc                                                         <  2:103072/81‑1 (MQ=255)              tCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCgcgc                                                       >  2:103106/1‑81 (MQ=255)                gATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTg                                                     >  2:103138/1‑81 (MQ=255)                 aTAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTgg                                                    <  1:103073/81‑1 (MQ=255)                      tGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGc                                               >  1:103121/1‑81 (MQ=255)                         cGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCg                                            <  1:103074/81‑1 (MQ=255)                          gTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCg                                            >  1:103113/1‑80 (MQ=255)                             gccgccTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAAt                                         <  1:103076/80‑1 (MQ=255)                             gccgccTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATg                                        <  1:103075/81‑1 (MQ=255)                                gccTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCa                                     >  1:103130/1‑81 (MQ=255)                                gccTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCa                                     >  2:103122/1‑81 (MQ=255)                                  cTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCAcc                                   >  1:103110/1‑81 (MQ=255)                                      cGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAg                                >  1:103129/1‑80 (MQ=255)                                      cGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAg                                >  2:103128/1‑80 (MQ=255)                                           gTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCaa                             >  2:103118/1‑78 (MQ=255)                                              aGCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGAtt                       >  2:103108/1‑81 (MQ=255)                                               gCCAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGAtt                       >  2:103119/1‑80 (MQ=255)                                                ccAGACGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGAtt                       >  2:103120/1‑79 (MQ=255)                                                    aCGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCAt                 <  2:103077/81‑1 (MQ=255)                                                     cGCCATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATg                >  2:103134/1‑81 (MQ=255)                                                        cATTCCTTTTTTCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTcg             <  1:103078/81‑1 (MQ=255)                                                                ttttCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGt       >  2:103147/1‑79 (MQ=255)                                                                 tttCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAg    >  1:103127/1‑81 (MQ=255)                                                                   tCTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGc   <  2:103079/80‑1 (MQ=255)                                                                    cTGACACTGTGCCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGcc  <  2:103080/80‑1 (MQ=255)                                                                               |                                                                     ACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCACCCAGCAAGAGATTGAGCATGTCGCGCGGTTAGCC  >  NC_000913/470071‑470217

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