Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,202,382 T→G T60P (ACC→CCC)  iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,202,3820TG96.6% 238.8 / ‑0.8 59T60P (ACC→CCC) iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (26/31);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

ATTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCA  >  REL606/4202309‑4202454
                                                                         |                                                                        
aTTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCg                                                                        >  2:1744415/1‑76 (MQ=255)
   gCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTgg                                                                     >  2:367954/1‑76 (MQ=255)
    cccAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCggggg                                                                    >  1:155344/1‑75 (MQ=255)
      cAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAAt                                                                  <  2:1103375/76‑1 (MQ=255)
      cAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAAt                                                                  <  1:244771/76‑1 (MQ=255)
       aaTGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAAtt                                                                 <  1:1209205/76‑1 (MQ=255)
       aaTGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAAtt                                                                 <  1:676213/76‑1 (MQ=255)
          gTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTggg                                                              >  1:1105980/1‑76 (MQ=255)
           tCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGt                                                             <  2:1018908/76‑1 (MQ=255)
           tCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGt                                                             >  2:1822361/1‑76 (MQ=255)
            cccAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTa                                                            <  1:1506126/76‑1 (MQ=255)
            cccAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTa                                                            <  1:1347358/76‑1 (MQ=255)
            cccAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTa                                                            <  2:742177/76‑1 (MQ=255)
            cccAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTa                                                            <  2:454548/76‑1 (MQ=255)
             ccAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGCGGGGCGTGGAATTGGGTaa                                                           <  2:1638986/76‑1 (MQ=255)
                 ttCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCg                                                       <  2:1820091/76‑1 (MQ=255)
                 ttCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCg                                                       <  1:1259472/76‑1 (MQ=255)
                  tCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCgg                                                      <  2:1296772/76‑1 (MQ=255)
                  tCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCgg                                                      >  1:1441038/1‑76 (MQ=255)
                   cGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGc                                                     <  1:1717468/76‑1 (MQ=255)
                     ccGACGTGACGCACTAAGCCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCtt                                                   <  2:1370157/76‑1 (MQ=255)
                           tGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTgcg                                             <  1:1847065/76‑1 (MQ=255)
                               gCACGAAGCCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAg                                         <  1:1503937/76‑1 (MQ=255)
                               gCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAg                                         <  1:1715767/76‑1 (MQ=255)
                                cACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGt                                        >  2:1977783/1‑76 (MQ=255)
                                 aCGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGtt                                       <  1:1068066/76‑1 (MQ=255)
                                   gAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTcc                                     >  1:784237/1‑76 (MQ=255)
                                    aaaCCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCg                                    >  1:446579/1‑76 (MQ=255)
                                       cccTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGa                                 >  1:981535/1‑76 (MQ=255)
                                        ccTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAg                                <  2:63555/76‑1 (MQ=255)
                                            ctgTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGcc                            >  2:1789920/1‑76 (MQ=255)
                                            ctgTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGcc                            >  2:745060/1‑76 (MQ=255)
                                             tgTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCca                           >  1:478626/1‑76 (MQ=255)
                                             tgTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGCGCca                           >  2:1713918/1‑76 (MQ=255)
                                              gTTGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCcac                          >  1:1795165/1‑76 (MQ=255)
                                               ttGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCcaca                         >  2:95191/1‑76 (MQ=255)
                                                tGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCacac                        >  2:68502/1‑76 (MQ=255)
                                                tGCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCacac                        >  1:1364428/1‑76 (MQ=255)
                                                 gCATTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACt                       <  2:1650495/76‑1 (MQ=255)
                                                   aTTGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGc                     >  1:663251/1‑76 (MQ=255)
                                                    ttGTGGTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGcc                    >  2:523845/1‑76 (MQ=255)
                                                      gtgGTTCGCAGGCGGTGGGGCGTGTAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCAt                  <  2:435361/76‑1 (MQ=255)
                                                      gaggTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCAt                  >  1:1447046/3‑76 (MQ=255)
                                                         gTTAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTgg               >  2:131825/1‑76 (MQ=255)
                                                           tAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGAt             >  2:860156/1‑76 (MQ=255)
                                                           tAGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGAt             >  1:1906141/1‑76 (MQ=255)
                                                             gCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTc           <  2:370194/76‑1 (MQ=255)
                                                             gCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTc           <  1:1443752/76‑1 (MQ=255)
                                                              cAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCg          >  1:1679036/1‑76 (MQ=255)
                                                              cAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCg          <  1:1912956/76‑1 (MQ=255)
                                                              cAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCg          <  1:177158/76‑1 (MQ=255)
                                                              cAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCg          <  1:1768098/76‑1 (MQ=255)
                                                              cAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCg          <  2:927619/76‑1 (MQ=255)
                                                               aGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCgg         >  2:1664740/1‑76 (MQ=255)
                                                                ggCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGc        <  1:1154813/76‑1 (MQ=255)
                                                                ggCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGc        >  2:388328/1‑76 (MQ=255)
                                                                ggCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGc        <  1:913745/76‑1 (MQ=255)
                                                                ggCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGc        <  2:868540/76‑1 (MQ=255)
                                                                 gCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCa       >  2:104352/1‑76 (MQ=255)
                                                                      ggggCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCa  >  2:129549/1‑76 (MQ=255)
                                                                      ggggCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCa  <  2:1441038/76‑1 (MQ=255)
                                                                         |                                                                        
ATTGCCCAATGTCCCAGTTCGCCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCA  >  REL606/4202309‑4202454

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