Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,561,758	G→A	S993S (AGC→AGT)	hsdR ←	endonuclease R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,561,758	0	G	A	100.0%	184.0 / NA	48	S993S (AGC→AGT)	hsdR	endonuclease R
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/28);  total (20/28)

GTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACGCTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCGGTTTTC  >  REL606/4561688‑4561824                                                                       |                                                                   gTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTcac                                                               >  1:1631538/1‑76 (MQ=255)  ttGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTcacc                                                              >  2:1490073/1‑76 (MQ=255)    gCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTcaccac                                                            >  2:443185/1‑76 (MQ=255)     cACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACt                                                           >  2:382322/1‑76 (MQ=255)      aCCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACtt                                                          <  1:761897/76‑1 (MQ=255)      aCCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACtt                                                          <  1:1295998/76‑1 (MQ=255)      aCCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACtt                                                          >  1:1885938/1‑76 (MQ=255)           cGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCAtcg                                                     >  1:1551449/1‑76 (MQ=255)           cGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCAtcg                                                     >  2:118001/1‑76 (MQ=255)            gAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCAtcgt                                                    >  2:180505/1‑76 (MQ=255)               tCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGa                                                 >  2:439977/1‑76 (MQ=255)                cAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGat                                                <  1:1490073/76‑1 (MQ=255)                cAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGat                                                <  2:46546/76‑1 (MQ=255)                     gCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGa                                           <  1:1153880/76‑1 (MQ=255)                      cTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAg                                          <  2:1631538/76‑1 (MQ=255)                      cTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAg                                          >  2:237948/1‑76 (MQ=255)                      cTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAg                                          >  2:1290024/1‑76 (MQ=255)                      cTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAg                                          >  2:483617/1‑76 (MQ=255)                       ttCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAgg                                         <  2:1848603/76‑1 (MQ=255)                        tCGGGGAAATCCTGCGGCGTGACGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGa                                        <  1:1046089/76‑1 (MQ=255)                        tCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGa                                        <  2:1191174/76‑1 (MQ=255)                         cGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGaa                                       >  1:926833/1‑76 (MQ=255)                          gAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAg                                      <  1:1748232/76‑1 (MQ=255)                          gAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAg                                      >  2:474775/1‑76 (MQ=255)                             gAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATc                                   <  1:641975/76‑1 (MQ=255)                                 tCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGgcg                               <  1:1824223/76‑1 (MQ=255)                                       ggcgTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGcatcat                         <  2:1400916/76‑1 (MQ=255)                                         cgTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTc                       <  2:906583/76‑1 (MQ=255)                                         cgTGTCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTc                       <  1:1182151/76‑1 (MQ=255)                                           tgtCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAg                     <  2:737911/76‑1 (MQ=255)                                           tgtCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAg                     <  1:1594009/76‑1 (MQ=255)                                           tgtCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAg                     <  2:1845162/76‑1 (MQ=255)                                            gtCGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAgg                    <  1:1794037/76‑1 (MQ=255)                                              cGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgtt                  >  1:688787/1‑76 (MQ=255)                                              cGTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgtt                  >  1:1697496/1‑76 (MQ=255)                                               gTAATCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgttg                 >  1:204364/1‑76 (MQ=255)                                                 aaTCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGgttgtt               >  2:1104826/1‑76 (MQ=255)                                                  aTCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGttgttg              <  1:57437/76‑1 (MQ=255)                                                  aTCACCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGttgttg              <  1:763226/76‑1 (MQ=255)                                                     aCCGTACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATg           >  1:153588/1‑76 (MQ=255)                                                         tACAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCgg       <  1:1757044/76‑1 (MQ=255)                                                          aCAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCGGt      <  2:601566/76‑1 (MQ=255)                                                           cAGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCGGtt     >  1:1814866/1‑76 (MQ=255)                                                            aGCGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCGGttt    <  2:82124/76‑1 (MQ=255)                                                              cGATTTTACACTCACCTCTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCGGTTTTc  <  2:1671856/76‑1 (MQ=255)                                                              cGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCGGTTTTc  <  2:334052/76‑1 (MQ=255)                                                              cGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCGGTTTTc  <  1:443185/76‑1 (MQ=255)                                                              cGATTTTACACTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCGGTTTTc  <  2:1147284/76‑1 (MQ=255)                                                                       |                                                                   GTTGCACCAGCGAGTCAAAGGCTTCGAGGAAATCCTGCGGCGTGTCGTAATCACCGTACAGCGATTTTACGCTCACCACTTCATCGTCGATATCGAGGAAGATCGGCGCATCATTCAGGTTGTTGATGTCGGTTTTC  >  REL606/4561688‑4561824

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