Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,201,841 Δ1 bp coding (719/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,8390A.96.2% 116.1 / ‑1.3 26coding (721/825 nt)iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (10/15);  total (10/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

TCCTTCGCCGCTTTAATCACCATTGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCA  >  REL606/4201766‑4201907
                                                                         |                                                                    
tccttcGCCGCTTTAATCACCATTGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGaat                                                                    <  1:927748/76‑2 (MQ=255)
  cttcGCCGCTTTAATCACCATTGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGaatag                                                                  <  2:1147647/76‑4 (MQ=255)
        cgcTTTAATCACCATTGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTg                                                           <  2:379523/76‑1 (MQ=255)
        cgcTTGAATCACCATAGCGCCAAGCTCGGTCACGCTGTCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTg                                                           <  2:379506/76‑1 (MQ=255)
          cTTTAATCACCATTGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCg                                                         >  2:594673/1‑76 (MQ=255)
             tAATCACCATTGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCa                                                      <  1:234753/76‑1 (MQ=255)
              aaTCACCATTGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCaa                                                     <  1:726468/76‑1 (MQ=255)
               aTCACCATTGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCaaa                                                    <  2:622263/76‑1 (MQ=255)
                 cACCATTGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAACg                                                  <  1:471058/76‑1 (MQ=255)
                         cgcCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCACgg                                          <  1:81369/76‑1 (MQ=255)
                         cgcCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCACgg                                          >  1:681844/1‑76 (MQ=255)
                              aaCTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCACGGTGCTc                                     <  1:31797/76‑1 (MQ=255)
                                 tCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCGGCAAACGGTTCACGGTGCTCAt                                   <  1:353197/76‑1 (MQ=255)
                                  cGGTCACGCGGTCATCGGNAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCg                                 <  1:939810/76‑1 (MQ=255)
                                   ggTCACGCGGTCATCGGTAACACGTTAAATCGGTCGGGAATAGCAAGTGCGGCATACGGTTCACGGTGCTCATcta                                >  1:127474/1‑74 (MQ=255)
                                      cACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGaaaa                             <  2:389861/76‑1 (MQ=255)
                                         gcgGTCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATac                          >  2:111935/1‑76 (MQ=255)
                                         gcgGTCATCGGTAATACCTGCAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATac                          <  2:165970/76‑1 (MQ=255)
                                             tCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAATGGTTCACGCTGCTCATCGAAAATACACGc                      >  1:1077142/1‑76 (MQ=255)
                                             tCATCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATACACGc                      >  1:803032/1‑76 (MQ=255)
                                                tCGGTAATACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCTCGGTGCTCATCGAAAATACACGCTGc                   <  1:111935/76‑1 (MQ=255)
                                                   gTAATACGTGAAATCGGTCCAGAATCGAAATTGCGACAAACGGTTCCCGGTGCTCATCGAAAATACACACTGCAAg                >  1:132947/1‑76 (MQ=255)
                                                      aTACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGATAATACACGCTGCAAGGCa             <  2:1080724/76‑1 (MQ=255)
                                                      aTACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATACACGCTGTAAGGCa             <  1:1116572/76‑1 (MQ=255)
                                                       tACGTGAAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATACACGCTGCAAGGCaa            >  2:525601/1‑76 (MQ=255)
                                                           tGAAATCGGTCCGGAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATACACGCTGCAAGGCAACGTa        >  2:711082/1‑76 (MQ=255)
                                                             aaaTCGGTCCGGAATAGAAAATGCGTCAAACGGTTAACGGTGCCCATCGAAAATACACGCTGCAAGGCTACGTacc      >  1:127202/1‑74 (MQ=255)
                                                                 cGGTCCGGAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCa  >  1:735055/1‑76 (MQ=255)
                                                                 cGGTCCGGAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCa  <  1:293795/76‑1 (MQ=255)
                                                                         |                                                                    
TCCTTCGCCGCTTTAATCACCATTGCGCCAAACTCGGTCACGCGGTCATCGGTAATACGTGAAATCGGTCCGGAAATAGAAATTGCGGCAAACGGTTCACGGTGCTCATCGAAAATACACGCTGCAAGGCAACGTAGCCCCA  >  REL606/4201766‑4201907

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