Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,202,378 T→C H61R (CAC→CGC)  iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,202,3780TC97.7% 160.8 / NA 43H61R (CAC→CGC) iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/21);  minor base A (1/0);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAA  >  REL606/4202308‑4202450
                                                                      |                                                                        
gATTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTc                                                                     >  1:1310028/1‑76 (MQ=255)
gATTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTc                                                                     <  2:1295213/76‑1 (MQ=255)
   tGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTg                                                                  <  2:226154/76‑1 (MQ=255)
   tGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTg                                                                  >  1:1622955/1‑76 (MQ=255)
         aTGTCCCANTTCGCCGACCTGACGCACTAAACCGCGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTg                                                            <  2:281597/76‑1 (MQ=255)
          tGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTgg                                                           >  2:1551260/1‑76 (MQ=255)
            tCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGt                                                         <  2:639841/76‑1 (MQ=255)
             cccAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTa                                                        <  2:49202/76‑1 (MQ=255)
               cAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAAc                                                      >  1:1357880/1‑76 (MQ=255)
               cAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAANTGGGTAAc                                                      <  1:1639505/76‑1 (MQ=255)
                aGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAAcc                                                     >  2:215187/1‑76 (MQ=255)
                  ttCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCg                                                   >  1:493361/1‑76 (MQ=255)
                   tCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCgg                                                  >  2:1491511/1‑76 (MQ=255)
                           cTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTgc                                          >  1:1596549/1‑76 (MQ=255)
                           cTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTgc                                          <  2:1246093/76‑1 (MQ=255)
                           cTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTgc                                          >  2:1272033/1‑76 (MQ=255)
                            tGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTgcg                                         <  2:236300/76‑1 (MQ=255)
                                  aCGAAACCCTGCTGTTGCCTTNTGNTTAGCAGACGGCGGGTCGTGGNATTGGTTNANCCGGCTTGTTGCGCCAGtt                                   >  2:242988/1‑76 (MQ=255)
                                  aCGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGtt                                   >  2:760806/1‑76 (MQ=255)
                                   cGAACCCCTGCGGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTc                                  <  2:1609145/76‑1 (MQ=255)
                                     aaaCCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCg                                >  2:581847/1‑76 (MQ=255)
                                     aaaCCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCg                                >  2:823290/1‑76 (MQ=255)
                                     aaaCCCTGCTGTTGCATTGTGGGTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCAGGCTTGTTGCGCCAGTTCCg                                <  2:143662/76‑1 (MQ=255)
                                         ccTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAg                            <  2:1527695/76‑1 (MQ=255)
                                            gctgTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGc                         <  1:482276/76‑1 (MQ=255)
                                             ctgTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGcc                        >  2:514990/1‑76 (MQ=255)
                                             ctgTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGcc                        >  2:819510/1‑76 (MQ=255)
                                             ctgTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGcc                        >  2:881088/1‑76 (MQ=255)
                                              tgTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCca                       <  1:1348388/76‑1 (MQ=255)
                                              tgTTGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCca                       <  1:1318122/76‑1 (MQ=255)
                                                ttGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCcaca                     <  2:871218/76‑1 (MQ=255)
                                                ttGCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCcaca                     <  2:352137/76‑1 (MQ=255)
                                                  gCATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACt                   <  1:1033212/76‑1 (MQ=255)
                                                   cATTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTg                  >  1:653375/1‑76 (MQ=255)
                                                    aTTGTGGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGc                 >  1:1147081/1‑76 (MQ=255)
                                                       gtgGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCAt              <  2:1536617/76‑1 (MQ=255)
                                                        tgGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCAtt             >  2:184019/1‑76 (MQ=255)
                                                        tgGTTAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCAtt             <  2:877827/76‑1 (MQ=255)
                                                           tttgCAAGCGGAGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGa          >  1:1572320/4‑76 (MQ=255)
                                                           ttAGCAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGa          >  2:656392/1‑76 (MQ=255)
                                                               cAGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCg      <  2:361021/76‑1 (MQ=255)
                                                                aGGCGGCGGGTCGTGGAATTGGTTCACCCGGCTTGTTGCGCCAGTTCCGTGAGAGCGATACTGCCATTGGATTCgg     >  1:85314/1‑76 (MQ=255)
                                                                aGGCGGCGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCgg     <  1:227944/76‑1 (MQ=255)
                                                                   cggcggGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCaa  >  1:769599/1‑76 (MQ=255)
                                                                      |                                                                        
GATTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAA  >  REL606/4202308‑4202450

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