Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,559,680	G→A	G77G (GGC→GGT)	ddpC ←	D,D‑dipeptide ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,559,680	0	G	A	100.0%	87.9 / NA	29	G77G (GGC→GGT)	ddpC	D,D‑dipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/16);  total (13/16)

CGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACG  >  NC_000913/1559556‑1559812                                                                                                                             |                                                                                                                                     cGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTgcgc                                                                                                                        >  1:1048107/1‑139 (MQ=255)      aGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAg                                                                                                                   >  2:77016/1‑139 (MQ=255)             cGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAggcgg                                                                                                            >  1:674004/1‑139 (MQ=255)                  aGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCaaaa                                                                                                       >  1:795645/1‑139 (MQ=255)                       aCCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACC‑AACCAGTGCGCCGCAGAAGGCGGCAAAAggcggg                                                                                                 <  1:848667/139‑1 (MQ=255)                                    gCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                                     <  2:1015736/139‑1 (MQ=255)                                    gCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                                     <  2:422174/139‑1 (MQ=255)                                      aaTGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATg                                                                                   >  1:821046/1‑139 (MQ=255)                                         ggCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCg                                                                                <  1:77016/139‑1 (MQ=255)                                          gCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGt                                                                               >  2:398836/1‑139 (MQ=255)                                           cAACCACGACTAATCCGGCGAGAACTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGtt                                                                              <  1:600816/139‑1 (MQ=255)                                             aCCACGACTAATCCGGCGAGAATTGATTGCTGAATGCCGACCCGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACGAGTGCGCCGCAGAAGGCGGCAAAAGGCTGGCGGTTAAATCAATGGCGGTgg                                                                            >  1:750838/1‑139 (MQ=255)                                                    gtAATCCAGCGAGACTTCATTGCTGACGGCCGACCAGTACGCGGCTAAACAGATCGAGACCCACTTCAGCGGTACCGAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAACCAATGGCGTTGGGATCATg                                                                     <  2:750838/138‑1 (MQ=255)                                                          cGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCgc                                                               <  2:674004/139‑1 (MQ=255)                                                                   ttgattgCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACgg                                                      <  2:795645/139‑1 (MQ=255)                                                                             aCTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAatc                                            >  2:1024457/1‑139 (MQ=255)                                                                             aCTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAatc                                            >  2:429656/1‑139 (MQ=255)                                                                              cTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAatca                                           <  2:381225/139‑1 (MQ=255)                                                                                      cAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAAt                                   <  2:821046/139‑1 (MQ=255)                                                                                        gTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAg                                 <  1:1023628/139‑1 (MQ=255)                                                                                           cgcgGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCAt                              <  2:304233/139‑1 (MQ=255)                                                                                           cgcgGCTAAACAGATCGCGTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCAt                              <  2:1048107/139‑1 (MQ=255)                                                                                                          cgcgTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGCGCGGCTACATCAATGGCGTCGGGATCATGCGTCGCCTGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGAc               >  2:290156/1‑139 (MQ=255)                                                                                                             gTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                                     <  2:543195/66‑1 (MQ=255)                                                                                                             gTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCaa                                                                                     >  1:543195/1‑66 (MQ=255)                                                                                                             gTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATa                         >  2:454794/1‑126 (MQ=255)                                                                                                             gTCCCACTTCATCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATa                         <  1:454794/126‑1 (MQ=255)                                                                                                                       aTCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACg  <  2:645982/139‑1 (MQ=255)                                                                                                                       aTCGGTACCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACg  >  1:645982/1‑139 (MQ=255)                                                                                                                             |                                                                                                                                     CGGATAGACATCCGAGTAGCGAACCAATCATACCCGCAATGGCAACCACGACTAATCCGGCGAGAATTGATTGCTGACTGCCGACCAGTACGCGGCTAAACAGATCGCGTCCCACTTCATCGGTGCCAAACCAGTGCGCCGCAGAAGGCGGCAAAAGGCGGGCGGTTAAATCAATGGCGTTGGGATCATGCGTCGCCAGCCACGGTGAAAAAATCATCATCAATAGCATTAATACAATGATGACTGCGCTGGTCACG  >  NC_000913/1559556‑1559812

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