Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,202,462 A→G L33S (TTA→TCA)  iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,202,4620AG94.4% 128.1 / ‑1.0 36L33S (TTA→TCA) iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (15/19);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

TGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC  >  REL606/4202391‑4202534
                                                                       |                                                                        
tGGGTAACCCGGCTTGTTGCGCCAGTGCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCTCTCCAGTAAttt                                                                      <  2:918089/76‑1 (MQ=255)
 gggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTc                                                                     <  2:919877/76‑1 (MQ=255)
 gggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGCGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTc                                                                     >  1:1366244/1‑76 (MQ=255)
 gcgTAACCCGGCTTGTTGCGCCAGTTCCGTGACTGCCACACTGCCATTGGATTCGCCAATCCACTCCAGTGATTTc                                                                     <  2:1601407/74‑1 (MQ=255)
  ggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCa                                                                    <  2:1191951/76‑1 (MQ=255)
   gTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAg                                                                   >  1:1040790/1‑76 (MQ=255)
      aCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGcc                                                                >  1:1109984/1‑76 (MQ=255)
            cTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGt                                                          >  1:849131/1‑76 (MQ=255)
              tgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTa                                                        >  1:1104325/1‑76 (MQ=255)
               gttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTaa                                                       <  1:579930/76‑1 (MQ=255)
                 tgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAg                                                     <  1:1468223/76‑1 (MQ=255)
                  gcgcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGa                                                    <  1:1146750/76‑1 (MQ=255)
                       aGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGaa                                               <  1:15221/76‑1 (MQ=255)
                          tCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCt                                            <  1:1035176/76‑1 (MQ=255)
                           ccGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTg                                           <  2:1189962/76‑1 (MQ=255)
                             gTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTc                                         <  1:1019949/76‑1 (MQ=255)
                               gAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCa                                       <  1:385234/76‑1 (MQ=255)
                                  tGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTc                                    >  2:1714486/1‑76 (MQ=255)
                                      acacTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTg                                <  2:571231/76‑1 (MQ=255)
                                        acTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGt                              >  1:171130/1‑76 (MQ=255)
                                        acTGCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGCCGCTGGt                              >  1:846950/1‑76 (MQ=255)
                                           gCCATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCg                           >  1:204816/1‑76 (MQ=255)
                                            ccATTGGATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCgg                          >  1:1277156/1‑76 (MQ=255)
                                                 ggATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCa                     <  2:1768356/76‑1 (MQ=255)
                                                  gATTCGGCAATCCTCTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCaa                    <  2:1707548/76‑1 (MQ=255)
                                                  gATTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCaa                    <  1:449414/76‑1 (MQ=255)
                                                   aTTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAAc                   >  2:406968/1‑76 (MQ=255)
                                                   aTTCGGCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGCGGCaag                   >  2:250515/1‑75 (MQ=255)
                                                       ggCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAAcggcg               <  1:304478/76‑1 (MQ=255)
                                                       ggCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAAcggcg               <  2:550460/76‑1 (MQ=255)
                                                        gCAATCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAAcggcgg              >  2:1664634/1‑76 (MQ=255)
                                                           aTCCACTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCggggt           >  1:384454/1‑74 (MQ=255)
                                                                cTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTg      <  2:986012/76‑1 (MQ=255)
                                                                cTCCAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTg      <  2:659036/76‑1 (MQ=255)
                                                                  ccAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGcc    >  2:66096/1‑76 (MQ=255)
                                                                  ccAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGcc    >  1:219127/1‑76 (MQ=255)
                                                                   cAGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCcg   <  1:445249/76‑1 (MQ=255)
                                                                    aGTGATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCcgc  <  2:708596/76‑1 (MQ=255)
                                                                       |                                                                        
TGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC  >  REL606/4202391‑4202534

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