Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,202,203	Δ1 bp	coding (357/825 nt)	iclR ←	DNA‑binding transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,202,202	0	C	.	100.0%	189.3 / NA	39	coding (358/825 nt)	iclR	DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (17/22);  total (17/22)

GCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGG  >  REL606/4202127‑4202272                                                                            |                                                                       gCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCAATATTGAc                                                                        >  2:31029/1‑76 (MQ=255) gCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAACACCGCCATATTGAc                                                                        >  2:31034/1‑76 (MQ=255)    tCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGAcgtt                                                                     >  2:495965/1‑73 (MQ=255)       gATGCGTACACTGTACCGGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCg                                                                 <  2:247853/76‑1 (MQ=255)        aTGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGc                                                                >  1:504894/1‑76 (MQ=255)          gCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCa                                                              >  1:214026/1‑76 (MQ=255)          gCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCa                                                              <  1:951282/76‑1 (MQ=255)            gTACACTGTACCGGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGa                                                            <  2:767333/76‑1 (MQ=255)              acacTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGAct                                                          >  2:353133/1‑76 (MQ=255)              acacTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGAct                                                          <  1:286241/76‑1 (MQ=255)              acacTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGAct                                                          >  2:176059/1‑76 (MQ=255)                 cTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCtt                                                       >  1:941168/1‑76 (MQ=255)                 cTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCtt                                                       >  1:883411/1‑76 (MQ=255)                  tGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTc                                                      >  2:1059625/1‑76 (MQ=255)                   gTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTcc                                                     >  2:382980/1‑76 (MQ=255)                   gTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTcc                                                     >  2:659723/1‑76 (MQ=255)                   gTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTcc                                                     >  2:206060/1‑76 (MQ=255)                     aCCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCAt                                                   <  1:529710/76‑1 (MQ=255)                        tGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTaa                                                <  2:797185/76‑1 (MQ=255)                         ggTCGATAATAATCGCTTCGTGATCGTTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTaaa                                               <  2:766014/76‑1 (MQ=255)                            cGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTg                                            >  2:537521/1‑76 (MQ=255)                              aaaataatCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTGGCCAGACTCTTCCATTAAATTgcg                                          <  1:755956/74‑1 (MQ=255)                               taataaTCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTgcgc                                         <  1:201548/76‑1 (MQ=255)                                 atattCGTTTCTTGTTCGTTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCTTTAAATTGCGCAg                                       <  2:313699/76‑1 (MQ=255)                                 ataatCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAg                                       <  1:398192/76‑1 (MQ=255)                                   aatCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGa                                     >  2:557018/1‑76 (MQ=255)                                    atCTCTTCTTGTTCGTTTTGTTCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGAt                                    <  2:342731/76‑1 (MQ=255)                                            gtgtTCGTTTTGTTCAGGCTCCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGaa                            <  1:284448/76‑1 (MQ=255)                                            gTGATCGCTTGGACCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGaa                            <  2:111866/76‑1 (MQ=255)                                             tGATCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAAc                           <  1:271026/76‑1 (MQ=255)                                               aTCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACaa                         <  1:961017/76‑1 (MQ=255)                                               aTCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACaa                         <  2:1200966/76‑1 (MQ=255)                                               aTCGCTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACaa                         >  1:503256/1‑76 (MQ=255)                                                   cTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGc                     >  1:760744/1‑76 (MQ=255)                                                   cTTTGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGc                     >  1:1027429/1‑76 (MQ=255)                                                    tttGATCAAGCCCCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCt                    <  2:327580/76‑1 (MQ=255)                                                     tggATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCTa                   <  2:13848/74‑1 (MQ=255)                                                      tGATCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCTaa                  <  1:476192/76‑1 (MQ=255)                                                      tGACCAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCTaa                  <  1:476206/76‑1 (MQ=255)                                                       gATCAAGCCCCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCTAAc                 <  2:87487/76‑1 (MQ=255)                                                          cAAGCACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCTAACaaa              >  2:700548/1‑76 (MQ=255)                                                              cACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTAc          <  2:547513/76‑1 (MQ=255)                                                              cACCGCCATATTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTAc          <  2:1100985/76‑1 (MQ=255)                                                                     atatTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTg   >  1:1194595/1‑76 (MQ=255)                                                                      tatTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTgg  <  1:455592/76‑1 (MQ=255)                                                                      tatTGA‑CGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTgg  >  1:1115228/1‑76 (MQ=255)                                                                            |                                                                       GCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAAATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGG  >  REL606/4202127‑4202272

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