Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,326,660	G→A	P36L (CCA→CTA)	mreD ←	cell wall structural complex MreBCD transmembrane component MreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,326,660	0	G	A	100.0%	214.9 / NA	54	P36L (CCA→CTA)	mreD	cell wall structural complex MreBCD transmembrane component MreD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/34);  total (20/34)

AAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGA  >  REL606/3326588‑3326732                                                                         |                                                                         aaaCCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGcc                                                                       >  1:1326246/1‑76 (MQ=255)  aaCCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCg                                                                      <  2:838759/76‑1 (MQ=255)   aCCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCgg                                                                     >  1:880798/1‑76 (MQ=255)     cTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGaa                                                                   >  1:79750/1‑76 (MQ=255)       gtgCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGaaaa                                                                 <  1:652792/76‑1 (MQ=255)       gtgCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGaaaa                                                                 >  1:232772/1‑76 (MQ=255)          cccACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACaa                                                              <  1:1422370/76‑1 (MQ=255)           ccACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAGCACCCAGTTTAGCCGGAAAACAAt                                                             <  1:528239/76‑1 (MQ=255)            cacaTTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATc                                                            <  1:168375/76‑1 (MQ=255)              caTTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAg                                                          <  2:1020761/76‑1 (MQ=255)               aTTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAgg                                                         <  1:1528613/76‑1 (MQ=255)                  tACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTa                                                      <  2:1295952/76‑1 (MQ=255)                   aCGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTAt                                                     <  1:1675140/76‑1 (MQ=255)                     gcgATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATcc                                                   <  2:58925/76‑1 (MQ=255)                      cgATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCg                                                  <  2:1665566/76‑1 (MQ=255)                        aTGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGc                                                >  1:1340430/1‑76 (MQ=255)                        aTGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGc                                                <  1:938897/76‑1 (MQ=255)                         ttaGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGcc                                               >  1:1599128/3‑76 (MQ=255)                         tggggCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGcc                                               <  2:1756003/73‑1 (MQ=255)                         tGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGcc                                               <  2:1037046/76‑1 (MQ=255)                         tGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGcc                                               <  2:1534092/76‑1 (MQ=255)                         tGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGcc                                               <  2:1136872/76‑1 (MQ=255)                          gAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCa                                              <  2:193938/76‑1 (MQ=255)                           aGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAg                                             <  2:340443/76‑1 (MQ=255)                            ggCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAgg                                            >  1:1004228/1‑76 (MQ=255)                                 ggCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATg                                       <  2:916662/76‑1 (MQ=255)                                  gCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGa                                      <  2:1525617/76‑1 (MQ=255)                                  gCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGa                                      <  2:348334/76‑1 (MQ=255)                                   ccAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGAt                                     >  1:984960/1‑76 (MQ=255)                                    cAGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGAtt                                    >  1:854819/1‑76 (MQ=255)                                     aGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGAttt                                   >  2:359309/1‑76 (MQ=255)                                     aGGATCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGAttt                                   >  1:1125779/1‑76 (MQ=255)                                        aTCCAATACAACAAGATGAGTAACCCCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTgca                                <  1:759693/76‑1 (MQ=255)                                         tCCAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTgcag                               >  1:98673/1‑76 (MQ=255)                                           cAATACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTgcagca                             <  2:1121326/76‑1 (MQ=255)                                            aaTACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCaa                            >  1:1642698/1‑76 (MQ=255)                                            aaTACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCaa                            >  1:1136207/1‑76 (MQ=255)                                             aTACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAAc                           <  1:1568255/76‑1 (MQ=255)                                             aTACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAAc                           <  1:1190249/76‑1 (MQ=255)                                              tACAACAAGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACa                          <  1:113530/76‑1 (MQ=255)                                                  acaaGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAgcgc                      >  2:996449/1‑76 (MQ=255)                                                   caaGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCa                     >  1:964288/1‑76 (MQ=255)                                                   caaGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCa                     <  1:940410/76‑1 (MQ=255)                                                   caaGATGAGTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCa                     <  1:183902/76‑1 (MQ=255)                                                           gTAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGaa             <  2:1075951/76‑1 (MQ=255)                                                           atAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGaa             <  1:67255/75‑1 (MQ=255)                                                            tAACACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGaaa            >  1:536795/1‑76 (MQ=255)                                                             aaCACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAag           <  2:1011372/76‑1 (MQ=255)                                                             aaCACCCAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAag           <  1:1064477/76‑1 (MQ=255)                                                                   cAGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGcc     <  2:639183/76‑1 (MQ=255)                                                                    aGTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCa    <  2:884921/76‑1 (MQ=255)                                                                     gTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAg   >  1:500088/1‑76 (MQ=255)                                                                     gTTTAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAg   >  1:224628/1‑76 (MQ=255)                                                                      tttAGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGa  >  1:685100/1‑76 (MQ=255)                                                                         |                                                                         AAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAACAAGATGAGTAACACCCAGTTTGGCCGGAAAACAATCAGGTTATCCGGCCAGGGCATGATTTGCAGCAACAGCGCAATGAGGAAAGAGAGCCAGA  >  REL606/3326588‑3326732

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