Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,201,945 Δ1 bp coding (615/825 nt) iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,201,9450G.96.4% 258.9 / 4.8 55coding (615/825 nt)iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (21/32);  total (22/33)
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.64e-01

AAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT  >  REL606/4201877‑4202019
                                                                    |                                                                          
aaaaTACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGACTAACCCCGTTTCGCGTTTg                                                                    >  2:772802/1‑76 (MQ=255)
aaaaTACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTg                                                                    >  2:1806062/1‑76 (MQ=255)
 aaaTACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTgg                                                                   >  2:1820924/1‑76 (MQ=255)
  aaTACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTggg                                                                  <  2:1274565/76‑1 (MQ=255)
   aTACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGc                                                                 <  2:620398/76‑1 (MQ=255)
     acacGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCga                                                               <  1:628559/76‑1 (MQ=255)
       acGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCgaga                                                             >  2:1663850/1‑76 (MQ=255)
        cGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGGTTGGGCGAGAt                                                            >  1:874355/1‑76 (MQ=255)
         gCTGCAAGGCAACGTAGCCCCAGCGCATGTNCCTCATCGTCAAATGAATNACNCCGTTTCGCGTTTGNGCGAGAtc                                                           >  1:1911718/1‑76 (MQ=255)
           tGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGAtctt                                                         <  2:291528/76‑1 (MQ=255)
            gCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGAtcttc                                                        >  2:603842/1‑76 (MQ=255)
             cAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGAtcttct                                                       <  1:651124/76‑1 (MQ=255)
              aaGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATcttctt                                                      <  1:2012583/76‑1 (MQ=255)
               aGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCttt                                                     >  2:1952748/1‑76 (MQ=255)
               aGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCttt                                                     <  2:551027/76‑1 (MQ=255)
                ggCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCtttt                                                    <  2:755249/76‑1 (MQ=255)
                 gCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTa                                                   <  1:1916432/76‑1 (MQ=255)
                 gCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTa                                                   <  1:742411/76‑1 (MQ=255)
                     cGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAAt                                                <  2:1199419/76‑1 (MQ=255)
                     cGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATg                                               <  1:278121/76‑1 (MQ=255)
                       tAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGca                                             <  2:930548/76‑1 (MQ=255)
                       tAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGca                                             >  1:95570/1‑76 (MQ=255)
                         gCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGcaca                                           <  2:270643/76‑1 (MQ=255)
                          ccccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGcacan                                          <  1:216891/76‑2 (MQ=255)
                          ccccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAg                                          >  2:1652053/1‑76 (MQ=255)
                           cccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAgg                                         <  2:402059/76‑1 (MQ=255)
                           cccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAgg                                         <  2:1385163/76‑1 (MQ=255)
                             cAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAggcg                                       >  1:900683/1‑74 (MQ=255)
                             cAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGgag                                       <  1:978456/76‑1 (MQ=255)
                                cgcATGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGaca                                    <  2:1397487/76‑1 (MQ=255)
                                   aTGTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCa                                 <  1:1370097/76‑1 (MQ=255)
                                     gTTCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGc                               <  1:1730890/76‑1 (MQ=255)
                                      ttCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCg                              <  2:1106039/76‑1 (MQ=255)
                                      ttCCTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCg                              >  2:1737184/1‑76 (MQ=255)
                                        ccTCATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGtt                            <  1:234599/76‑1 (MQ=255)
                                           cATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCg                         <  2:246452/76‑1 (MQ=255)
                                           cATCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCg                         >  2:703231/1‑76 (MQ=255)
                                            aTCGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGt                        >  2:120182/1‑76 (MQ=255)
                                              cGTCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTgg                      >  1:2021413/1‑76 (MQ=255)
                                                tCAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGt                    <  2:865460/76‑1 (MQ=255)
                                                 cAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTgggca                   >  2:648531/1‑74 (MQ=255)
                                                 cAAATGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGta                   <  2:52485/76‑1 (MQ=255)
                                                  aaaTGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGtat                  <  2:1795650/76‑1 (MQ=255)
                                                  aaaTGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGtat                  <  2:1600932/76‑1 (MQ=255)
                                                  aaaTGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGtat                  <  2:1696265/76‑1 (MQ=255)
                                                  aaaTGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGgat                  >  1:1025908/1‑76 (MQ=255)
                                                   aaTGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGtata                 <  2:1902794/76‑1 (MQ=255)
                                                   aaTGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGtata                 >  2:508425/1‑76 (MQ=255)
                                                     tGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAgg               >  1:287492/1‑76 (MQ=255)
                                                     tGAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAgg               <  1:1592159/76‑1 (MQ=255)
                                                      gAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTTAATGCACAGGAGACACCAGCGTTGCGTGGGGATAGGc              >  1:995369/1‑76 (MQ=255)
                                                      gAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGc              >  1:1908088/1‑76 (MQ=255)
                                                      gAATAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGc              >  1:161596/1‑76 (MQ=255)
                                                      aaaTAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGc              <  1:867558/75‑1 (MQ=255)
                                                         tAACCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATg           <  1:1637951/76‑1 (MQ=255)
                                                           aCCCCGTTTCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTa         <  2:968249/76‑1 (MQ=255)
                                                                   tGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCttt  >  1:1658305/1‑76 (MQ=255)
                                                                    |                                                                          
AAAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT  >  REL606/4201877‑4202019

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