Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	875,094	G→C	H65D (CAT→GAT)	deoR ←	DNA‑binding transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	875,094	0	G	C	100.0%	220.8 / NA	58	H65D (CAT→GAT)	deoR	DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (22/36);  total (22/36)

GCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGA  >  REL606/875047‑875141                                                |                                                gCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAAtcgct                                              >  1:4044990/1‑51 (MQ=255)  cGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAAtcgctg                                             >  1:5433766/1‑51 (MQ=255)  cGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAAtcgctg                                             >  1:2968893/1‑51 (MQ=255)   gTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAAtcgctgg                                            >  1:2866592/1‑51 (MQ=255)     ttttCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCa                                          <  1:3237575/51‑1 (MQ=255)     ttttCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCa                                          <  1:5281772/51‑1 (MQ=255)        tcttcCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTg                                       <  1:5835141/51‑1 (MQ=255)        tcttcCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTg                                       <  1:2063300/51‑1 (MQ=255)        tcttcCACCAGGCGGGATTTTTGATCGCTTAACAGGGAAtcgctgg                                            >  1:6737972/1‑46 (MQ=39)           tcCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTgcgc                                    <  1:5383747/51‑1 (MQ=255)            ccaccaGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGcgcg                                   <  1:7258283/51‑1 (MQ=255)            ccaccaGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGcgcg                                   <  1:6701745/51‑1 (MQ=255)              accaGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGt                                 <  1:5450052/51‑1 (MQ=255)              accaGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGt                                 <  1:605457/51‑1 (MQ=255)               ccaGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGtt                                >  1:6976454/1‑51 (MQ=255)               ccaGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGtt                                <  1:7106872/51‑1 (MQ=255)                caGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTc                               >  1:2500909/1‑51 (MQ=255)                 aGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGcgcg                                   <  1:6914952/46‑1 (MQ=255)                 aGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTcc                              <  1:6264347/51‑1 (MQ=255)                 aGGCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTcc                              >  1:2649016/1‑51 (MQ=255)                 aGGCGGGATTTTTGATCGCCTAACAGGTAATCGCTGGCACTGCGCGGTTcc                              <  1:7221434/51‑1 (MQ=255)                  ggCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCa                             >  1:1604950/1‑51 (MQ=255)                   gCGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAg                            <  1:3117558/51‑1 (MQ=255)                    cGGGATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGa                           <  1:3281620/51‑1 (MQ=255)                     gggATTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGaa                          <  1:886797/51‑1 (MQ=255)                        aTTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACaa                       <  1:3205702/51‑1 (MQ=255)                        aTTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACaa                       <  1:6463287/51‑1 (MQ=255)                        aTTTTTGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACaa                       <  1:4030724/51‑1 (MQ=255)                         tttttGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAat                      >  1:4600866/1‑51 (MQ=255)                         tttttGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAat                      >  1:3421928/1‑51 (MQ=255)                          ttttGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAata                     >  1:2169272/1‑51 (MQ=255)                           tttGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAatat                    >  1:5936384/1‑51 (MQ=255)                            ttGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAAtata                   <  1:6096750/51‑1 (MQ=255)                             tGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAg                  <  1:289601/51‑1 (MQ=255)                             tGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAg                  <  1:2352872/51‑1 (MQ=255)                             tGATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAg                  >  1:4762946/1‑51 (MQ=255)                              gATCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAgc                 <  1:6463478/51‑1 (MQ=255)                               aTCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAgcc                <  1:3420726/51‑1 (MQ=255)                               aTCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAgcc                >  1:5299279/1‑51 (MQ=255)                                tCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAgccg               <  1:3814454/51‑1 (MQ=255)                                tCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAgccg               >  1:5607894/1‑51 (MQ=255)                                tCGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAgccg               <  1:4356753/51‑1 (MQ=255)                                 cGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAgccgc              >  1:7096720/1‑51 (MQ=255)                                 cGCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAgccg               <  1:449540/50‑1 (MQ=255)                                  gCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAgccgcc             <  1:4994945/51‑1 (MQ=255)                                  gCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAgccgcc             >  1:6989268/1‑51 (MQ=255)                                  gCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAgccgcc             <  1:1642696/51‑1 (MQ=255)                                  gCTTAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAgccgcc             >  1:2628950/1‑51 (MQ=255)                                    ttAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGa           >  1:3391661/1‑51 (MQ=255)                                    ttAACAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGa           <  1:2840769/51‑1 (MQ=255)                                      aaCAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGc         <  1:5752463/51‑1 (MQ=255)                                      aaCAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGc         >  1:6685788/1‑51 (MQ=255)                                      aaCAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGc         >  1:2956990/1‑51 (MQ=255)                                       aCAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCa        >  1:2176591/1‑51 (MQ=255)                                        cAGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCaa       <  1:3878212/51‑1 (MQ=255)                                         aGGTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCaaa      >  1:1260824/1‑51 (MQ=255)                                          ggTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCaaaa     <  1:4651162/51‑1 (MQ=255)                                           gTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAac    <  1:4635726/51‑1 (MQ=255)                                           gTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAac    >  1:6230467/1‑51 (MQ=255)                                           gTAATCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAac    <  1:3160640/51‑1 (MQ=255)                                             aaTCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAacga  >  1:6318129/1‑51 (MQ=255)                                             aaTCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAacga  <  1:1986052/51‑1 (MQ=255)                                             aaTCGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAacga  >  1:2741620/1‑51 (MQ=255)                                                |                                                GCGTTTTTCTTCCACCAGGCGGGATTTTTGATCGCTTAACAGGTAATGGCTGGCACTGCGCGGTTCCAGAACAATATAGCCGCCGAGCAAAACGA  >  REL606/875047‑875141

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