Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,302,608	T→C	K158R (AAG→AGG)	dacD ←	D‑alanyl‑D‑alanine carboxypeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,302,608	0	T	C	100.0%	268.0 / NA	72	K158R (AAG→AGG)	dacD	D‑alanyl‑D‑alanine carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/40);  total (32/40)

GAGCTATGCTGGCCAGGTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCTTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGTCAGCCAG  >  minE/1302542‑1302672                                                                   |                                                                 gAGCTATGCTGGCCAGGTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCtctc                                                              <  1:2597699/71‑1 (MQ=255)  aGCTATGCTGGCCAGGTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCg                                                             <  1:1421496/71‑1 (MQ=255)      aTGCTGGCCAGGTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCag                                                         >  1:1108206/1‑70 (MQ=255)        gCTGGCCAGGTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCAt                                                         <  1:2448311/69‑1 (MQ=255)         cTGGCCAGGTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCg                                                             <  1:1420203/64‑1 (MQ=255)         cTGGCCAGGTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCg                                                             <  1:2494854/64‑1 (MQ=255)            gCCAGGTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAgttgt                                                   <  1:569131/71‑1 (MQ=255)            gCCAGGTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAgttgt                                                   <  1:3092503/71‑1 (MQ=255)               aGGTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTtca                                                >  1:115277/1‑71 (MQ=255)                ggTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTtcat                                               >  1:790468/1‑71 (MQ=255)                ggTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTtcat                                               >  1:1274214/1‑71 (MQ=255)                 gTGCATCCAGACAATGCACTGTTTCAAAATGCGAATCCTTGAGATGCAGCctct                                                               >  1:515322/1‑54 (MQ=255)                        cAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAgttgt                                                   <  1:644045/59‑1 (MQ=255)                        cAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTcaa                                       >  1:801952/1‑71 (MQ=255)                             cATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAAc                                  >  1:3002065/1‑71 (MQ=255)                               tGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAgttgt                                                   <  1:3178225/52‑1 (MQ=255)                               tGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACaaa                                   >  1:979102/1‑68 (MQ=255)                               tGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACaaa                                   >  1:1202767/1‑68 (MQ=255)                                 cACTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGc                               <  1:1389054/70‑1 (MQ=255)                                  aCTGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGc                               <  1:3162376/69‑1 (MQ=255)                                    tGTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAgt                                                      <  1:1189374/44‑1 (MQ=255)                                     gTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTcatcat                                            <  1:1886514/53‑1 (MQ=255)                                     gTTTCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTg                          >  1:532072/1‑71 (MQ=255)                                      tttCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGtt                           <  1:1326835/69‑1 (MQ=255)                                      tttCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGtt                           <  1:2528255/69‑1 (MQ=255)                                      tttCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTg                          <  1:705966/70‑1 (MQ=255)                                      tttCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTg                          <  1:2732844/70‑1 (MQ=255)                                      tttCAAAATGCGTATCCTAGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGtt                           <  1:525028/69‑1 (MQ=255)                                       ttCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCAttt                                          >  1:3149170/1‑53 (MQ=255)                                       ttCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCAttt                                          >  1:2906201/1‑53 (MQ=255)                                       ttCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCAttt                                          >  1:2732824/1‑53 (MQ=255)                                       ttCAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGc                         <  1:1456348/70‑1 (MQ=255)                                         cAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGccc                       <  1:2989406/70‑1 (MQ=255)                                         cAAAATGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGccc                       <  1:2937803/70‑1 (MQ=255)                                          aaaaTGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTtcat                                               >  1:3152307/1‑45 (MQ=255)                                          aaaaTGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTtcat                                               >  1:2192333/1‑45 (MQ=255)                                          aaaaTGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCAc                     >  1:2534333/1‑71 (MQ=255)                                           aaaTGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCAcc                    >  1:1960527/1‑71 (MQ=255)                                           aaaTGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCAcc                    >  1:414555/1‑71 (MQ=255)                                           aaaTGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCAcc                    >  1:1697682/1‑71 (MQ=255)                                           aaaTGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCAcc                    >  1:1805738/1‑71 (MQ=255)                                           aaaTGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCAcc                    >  1:2185985/1‑71 (MQ=255)                                            aaTGCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCg                   >  1:633698/1‑71 (MQ=255)                                               gCGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGc                 >  1:2562047/1‑70 (MQ=255)                                                cGTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCaa               <  1:1300530/71‑1 (MQ=255)                                                 gTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACaaa                                   <  1:2469658/50‑1 (MQ=255)                                                 gTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACaaa                                   <  1:168512/50‑1 (MQ=255)                                                 gTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACaaa                                   <  1:2213156/50‑1 (MQ=255)                                                 gTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACaaa                                   <  1:1494045/50‑1 (MQ=255)                                                 gTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACaaa                                   <  1:1831471/50‑1 (MQ=255)                                                 gTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACaaa                                   <  1:2045255/50‑1 (MQ=255)                                                 gTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACaaa                                   <  1:2017105/50‑1 (MQ=255)                                                 gTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAat              <  1:1124330/71‑1 (MQ=255)                                                 gTATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATATCAACAAACTGCCGTTGCCCACCGGCAat              <  1:883265/71‑1 (MQ=255)                                                  tATCCTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAata             <  1:1128143/71‑1 (MQ=255)                                                      cTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCAc                     >  1:11790/1‑59 (MQ=255)                                                      cTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCAc                     >  1:1479844/1‑59 (MQ=255)                                                      cTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCAc                     >  1:3163872/1‑59 (MQ=255)                                                      cTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCAc                     >  1:1084040/1‑59 (MQ=255)                                                      cTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCAc                     >  1:837501/1‑59 (MQ=255)                                                      cTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCAc                     >  1:1479747/1‑59 (MQ=255)                                                      cTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGt         <  1:420585/71‑1 (MQ=255)                                                      cTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGt         <  1:2311635/71‑1 (MQ=255)                                                      cTTGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGt         <  1:1886441/71‑1 (MQ=255)                                                       ttGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGt         >  1:258648/1‑70 (MQ=255)                                                       ttGAGATGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGt         >  1:2006688/1‑70 (MQ=255)                                                         gagaTGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGTCa       <  1:1345710/70‑1 (MQ=255)                                                            aTGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGTCAGCCa   <  1:1016199/71‑1 (MQ=255)                                                             tGCAGCCTCTCGGCATATTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGTCAGCCag  <  1:1482932/71‑1 (MQ=255)                                                             tGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGTCAGCCag  <  1:3108966/71‑1 (MQ=255)                                                             tGCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGTCAGCCag  <  1:194399/71‑1 (MQ=255)                                                              gCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAatat            >  1:261733/1‑60 (MQ=255)                                                              gCAGCCTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGTCAGCCa   >  1:2609525/1‑69 (MQ=255)                                                                   |                                                                 GAGCTATGCTGGCCAGGTGCATCCAGACCATGCACTGTTTCAAAATGCGTATCCTTGAGATGCAGCTTCTCGGCATAGTTGTTCATCATTTCAACAAACTGCCGTTGCCCACCGGCAATATAGTCAGCCAG  >  minE/1302542‑1302672

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