Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,326,561	A→T	I69N (ATC→AAC)	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,561	0	A	T	100.0%	116.8 / NA	33	I69N (ATC→AAC)	mreD	cell wall structural complex MreBCD transmembrane component MreD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (17/16);  total (17/16)

CTGGTATTTCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCA  >  REL606/3326488‑3326631                                                                          |                                                                       cTGGTATTTCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTc                                                                      >  1:181123/1‑76 (MQ=255)    gTATTTCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGa                                                                   <  1:1028623/76‑1 (MQ=255)    gTATTTCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGa                                                                   <  1:374526/76‑1 (MQ=255)       tttCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATcc                                                                <  2:102724/76‑1 (MQ=255)        ttCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCa                                                               >  2:385765/1‑76 (MQ=255)        ttCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGAACGCCAAGCGTCGAGCCGCTGTTCAGATCCa                                                               >  2:996769/1‑76 (MQ=255)         tCAGCGCCACCAGGTAANCAATGATGCTCATCGCCAATACGCGTACGCCAANCGTCGACCCGCCGTTCAGATCCAg                                                              >  2:898448/1‑76 (MQ=255)            acgcCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTNAGATCCTGTAt                                                           <  1:153686/75‑1 (MQ=255)             cgcCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATg                                                          >  2:1527/1‑76 (MQ=255)             cgcCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATg                                                          <  1:98022/76‑1 (MQ=255)              gcCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATgg                                                         <  1:463362/76‑1 (MQ=255)                    aGGTAAGCAATGACGCCCACCGCCAAAACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCa                                                   <  1:260553/76‑1 (MQ=255)                       tAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATca                                                >  1:887478/1‑76 (MQ=255)                           cAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACaaa                                            <  1:545937/76‑1 (MQ=255)                           cAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACaaa                                            >  2:388468/1‑76 (MQ=255)                            aaTGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACaaaa                                           >  2:29604/1‑76 (MQ=255)                              tgatgCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACNAAAcc                                         >  2:876102/1‑76 (MQ=255)                              tgatgCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAAcc                                         >  2:612117/1‑76 (MQ=255)                                     cATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGccc                                  >  2:389313/1‑76 (MQ=255)                                        cGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCcaca                               <  1:898652/76‑1 (MQ=255)                                        cGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCcaca                               >  2:1071770/1‑76 (MQ=255)                                        cGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCcaca                               >  2:603885/1‑76 (MQ=255)                                            aaTACGCGTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTa                           <  1:333232/76‑1 (MQ=255)                                                  cgTGCGCCCACCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGAt                     <  2:12069/76‑1 (MQ=255)                                                  cgTACGCCAAGCGTCGAGCCGCTNTTCAGATNCNGTATGGAACCCCTCACAAAACCTGAGCCCAGATTTATGCGAt                     >  1:1004489/1‑76 (MQ=255)                                                   gTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGTCCACATTTACGCGATg                    >  2:443872/1‑76 (MQ=255)                                                   gTACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATg                    >  1:740953/1‑76 (MQ=255)                                                    tACGCCAAGCGTCGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGa                   <  2:509252/76‑1 (MQ=255)                                                          aaGCGGCGAGCCGCTGTTCAGATCCAGTATGGCTCCCATCACAAAACCTGTGCTCTCATTTACGCGCTGAGGCAAg             <  1:107019/76‑1 (MQ=255)                                                              gtggAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCa         <  2:1029804/73‑1 (MQ=255)                                                                cGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAgg       <  2:311768/76‑1 (MQ=255)                                                                cGAGCCGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAgg       >  1:594006/1‑76 (MQ=255)                                                                     cGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCa  <  2:226320/76‑1 (MQ=255)                                                                     cGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCa  >  2:206051/1‑76 (MQ=255)                                                                     cGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCa  >  2:642455/1‑76 (MQ=255)                                                                     cGCTGTTCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCa  <  1:799157/76‑1 (MQ=255)                                                                          |                                                                       CTGGTATTTCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCA  >  REL606/3326488‑3326631

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