Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,202,412	C→T	A50T (GCG→ACG)	iclR ←	DNA‑binding transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,202,412	0	C	T	100.0%	139.1 / NA	37	A50T (GCG→ACG)	iclR	DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/17);  total (20/17)

GCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCG  >  REL606/4202340‑4202477                                                                         |                                                                  gCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAg                                                                <  1:1429394/76‑1 (MQ=255)   aCGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGtt                                                              <  2:1205384/76‑1 (MQ=255)    cGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTc                                                             <  1:1019844/76‑1 (MQ=255)    cGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTc                                                             >  1:1219202/1‑76 (MQ=255)    cGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTc                                                             >  1:1401269/1‑76 (MQ=255)    cGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTc                                                             >  2:1409127/1‑76 (MQ=255)      aaaCCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCg                                                           >  2:1353949/1‑76 (MQ=255)      aaaCCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCg                                                           >  2:1031042/1‑76 (MQ=255)        aCCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCNGGCTTGTTGCGTCAGTTCCGTg                                                         <  2:238063/76‑1 (MQ=255)                gTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCcac                                                 >  1:576271/1‑76 (MQ=255)                 ttGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCcaca                                                <  2:540185/76‑1 (MQ=255)                     aTTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGc                                            <  2:970970/76‑1 (MQ=255)                     aTTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGc                                            >  1:1539849/1‑76 (MQ=255)                       tgtgGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCa                                          <  1:1144768/76‑1 (MQ=255)                       tgtgGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCa                                          >  2:947980/1‑76 (MQ=255)                        gtgGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCAt                                         <  2:685749/76‑1 (MQ=255)                         tgGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCAtt                                        <  1:1317735/76‑1 (MQ=255)                          ggTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTg                                       <  2:1714711/76‑1 (MQ=255)                                cAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCg                                 >  1:323075/1‑76 (MQ=255)                                 aGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCgg                                >  2:1125338/1‑76 (MQ=255)                                 aGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCgg                                >  1:290137/1‑76 (MQ=255)                                  ggCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGc                               <  2:55999/76‑1 (MQ=255)                                     ggTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAAt                            >  1:999209/1‑76 (MQ=255)                                       tGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATcc                          >  1:705340/1‑76 (MQ=255)                                            cGTGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTcc                     <  1:171762/76‑1 (MQ=255)                                              tGGAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAg                   >  2:199452/1‑76 (MQ=255)                                               ggAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGt                  >  2:213779/1‑76 (MQ=255)                                                gAATTGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTa                 <  2:891575/76‑1 (MQ=255)                                                    tGGGTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAAttt             >  1:307805/1‑76 (MQ=255)                                                    tGGGCAACCCGGCTTGTTGCGTCAGGTCCGCGAGTGCCCCACTGCCATTGGATCCAGCCATCCACTCCGGTAAttt             >  2:786757/1‑76 (MQ=255)                                                     gggTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTc            <  2:1235348/76‑1 (MQ=255)                                                     gggTAACCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTc            <  2:1212530/76‑1 (MQ=255)                                                         aaCCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGc        >  2:584859/1‑76 (MQ=255)                                                         aaCCCGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGc        >  1:1518359/1‑76 (MQ=255)                                                             cGGCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCAcg    <  1:438514/76‑1 (MQ=255)                                                              ggCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCAcgc   >  1:586821/1‑76 (MQ=255)                                                               gCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCANTCCAGTAATTTCAGGCCAcgcg  <  2:243301/76‑1 (MQ=255)                                                               gCTTGTTGCGTCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCAcgcg  >  1:190700/1‑76 (MQ=255)                                                                         |                                                                  GCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCG  >  REL606/4202340‑4202477

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: ‑