Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,201,969	T→G	L197F (TTA→TTC)	iclR ←	DNA‑binding transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,201,969	0	T	G	96.2%	106.2 / ‑2.6	26	L197F (TTA→TTC)	iclR	DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (12/13);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

AACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCAC  >  REL606/4201896‑4202040                                                                          |                                                                        aaCGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGaa                                                                       >  1:421651/1‑76 (MQ=255)   cGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAAt                                                                     <  1:80005/76‑1 (MQ=255)      aGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGca                                                                  >  1:412162/1‑76 (MQ=255)        ccccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGcaca                                                                <  1:345416/76‑1 (MQ=255)          ccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAgg                                                              <  2:351104/76‑1 (MQ=255)                  tGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCa                                                      >  1:457120/1‑76 (MQ=255)                  tGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGCTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCa                                                      >  1:704209/1‑76 (MQ=255)                       cTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGtt                                                 <  1:267082/76‑1 (MQ=255)                        tCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTg                                                <  2:422482/76‑1 (MQ=255)                          aTCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCTNCGTTGCg                                              <  2:677842/76‑1 (MQ=255)                            cGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTg                                            >  2:416213/1‑76 (MQ=255)                             gTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTgg                                           >  1:669064/1‑76 (MQ=255)                                  aTGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGtata                                      <  1:300442/76‑1 (MQ=255)                                  aTGAATAACCCCGCTTGCGCNTTTGGGCGATATCTTNTTTGAAATGCACAGGAGACACTAGCGTTGCGTGGGtata                                      <  2:837273/76‑1 (MQ=255)                                   tGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAg                                     <  1:322464/76‑1 (MQ=255)                                    gAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAgg                                    >  1:380841/1‑76 (MQ=255)                                       tAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCAt                                 <  1:595569/76‑1 (MQ=255)                                       tAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCAt                                 >  1:609744/1‑76 (MQ=255)                                       tAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCAt                                 >  2:430428/1‑76 (MQ=255)                                        aaCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATg                                <  1:613045/76‑1 (MQ=255)                                             cGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAAcc                           <  2:293578/76‑1 (MQ=255)                                               tttGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCt                         <  2:431497/76‑1 (MQ=255)                                                   cgcgTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGc                     >  2:82095/1‑76 (MQ=255)                                                          gggCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgcag              >  1:464068/1‑76 (MQ=255)                                                                gaTCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCtt        <  2:284678/76‑1 (MQ=255)                                                                      cttTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCAc  >  1:321915/1‑76 (MQ=255)                                                                          |                                                                        AACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCAC  >  REL606/4201896‑4202040

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