Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,414,325	A→C	K101Q (AAA→CAA)	crp →	DNA‑binding transcriptional dual regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,414,325	0	A	C	100.0%	189.6 / NA	45	K101Q (AAA→CAA)	crp	DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (20/25);  total (20/25)

GTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTC  >  REL606/3414252‑3414397                                                                          |                                                                         gTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACcaa                                                                        >  2:56015/1‑76 (MQ=255)  tttGAAGGGGGCCAGGAACTTAGCGCAGGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACcaaa                                                                       <  1:1489537/76‑1 (MQ=255)   tgaaaGGGGGCCGGGAACGTAGCGCAGGGGTACGGGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACcaaaa                                                                      <  1:1019303/73‑1 (MQ=255)     gAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAAt                                                                    <  2:458972/76‑1 (MQ=255)       agagGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAAttt                                                                  >  2:1064185/1‑76 (MQ=255)          gggCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGc                                                               <  2:466429/76‑1 (MQ=255)               aGGAACGTAGCGCATGGGTACGTGCGGAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAAtt                                                          >  1:1249430/1‑76 (MQ=255)               aGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAAtt                                                          <  1:1755759/76‑1 (MQ=255)               aGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAAtt                                                          >  1:723632/1‑76 (MQ=255)               aGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAAtt                                                          <  2:1082021/76‑1 (MQ=255)                   aCGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGAtt                                                      <  1:1567975/76‑1 (MQ=255)                     gTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCa                                                    <  2:401860/76‑1 (MQ=255)                       aGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAgg                                                  >  1:1112320/1‑76 (MQ=255)                         cgcATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTa                                                <  1:1267295/76‑1 (MQ=255)                          gcATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTaa                                               >  1:1340212/1‑76 (MQ=255)                          gcATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTaa                                               >  1:1331117/1‑76 (MQ=255)                           caggggTACGGGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTaaa                                              <  1:424808/73‑1 (MQ=255)                            gggggTGCGTGCAAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAAc                                             <  1:748537/74‑1 (MQ=255)                            aTGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAAc                                             >  1:1370769/1‑76 (MQ=255)                             tGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAAcc                                            <  2:645359/76‑1 (MQ=255)                              gggTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAAccc                                           <  1:103579/76‑1 (MQ=255)                              gggTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAAccc                                           >  1:1661406/1‑76 (MQ=255)                                  aCGTGCGAAAGCCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGAc                                       >  1:197732/1‑76 (MQ=255)                                  aCGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGAc                                       >  1:91265/1‑76 (MQ=255)                                   cGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACa                                      >  1:826153/1‑76 (MQ=255)                                     tGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACAtt                                    <  2:1468784/76‑1 (MQ=255)                                       ggAAATCCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCt                                  <  1:1768586/75‑1 (MQ=255)                                        gAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTg                                 >  1:1646785/1‑76 (MQ=255)                                         aaaaCCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGa                                >  1:1600680/1‑76 (MQ=255)                                         aaaaCCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGCACATTCTGa                                <  1:146496/76‑1 (MQ=255)                                           aaCCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATg                              <  1:793035/76‑1 (MQ=255)                                           aaCCGCCTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATg                              <  1:1621106/76‑1 (MQ=255)                                                 cTGTGAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTg                        >  2:464540/1‑76 (MQ=255)                                                   gtgAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTc                      >  1:295138/1‑76 (MQ=255)                                                    tgAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCt                     >  1:126405/1‑76 (MQ=255)                                                     gAAGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTg                    <  2:1170269/76‑1 (MQ=255)                                                      aaGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGc                   <  2:1753206/76‑1 (MQ=255)                                                      aaGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGc                   >  1:1392174/1‑76 (MQ=255)                                                       aGTGGCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGca                  <  2:1270284/76‑1 (MQ=255)                                                          ggCTGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAg               <  2:1755011/76‑1 (MQ=255)                                                             tGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATg            <  2:369081/76‑1 (MQ=255)                                                             tGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATg            <  2:544883/76‑1 (MQ=255)                                                             tGAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATg            <  1:1044945/76‑1 (MQ=255)                                                              gAAATTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATgg           >  1:1364151/1‑76 (MQ=255)                                                               aaaTTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGgc          >  2:581498/1‑76 (MQ=255)                                                                 aTTTCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGgcgc        <  2:1213518/76‑1 (MQ=255)                                                                   ttCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGcgt      >  2:716982/1‑76 (MQ=255)                                                                    tCGTACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGcgtc     <  1:1335790/76‑1 (MQ=255)                                                                       tACCAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCgtcgtc  >  2:1073958/1‑76 (MQ=255)                                                                          |                                                                         GTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCCGGACATTCTGATGCGTTTGTCTGCACAGATGGCGCGTCGTC  >  REL606/3414252‑3414397

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