Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,202,096 G→A P155L (CCG→CTG)  iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,202,0960GA96.2% 85.6 / ‑3.7 26P155L (CCG→CTG) iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (19/6);  total (20/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAAT  >  REL606/4202025‑4202160
                                                                       |                                                                
gcagcaGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCaa                                                              <  1:382956/76‑1 (MQ=255)
gcagcaGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCaa                                                              <  2:339963/76‑1 (MQ=255)
 cagcagCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCAAt                                                             >  1:213174/1‑76 (MQ=255)
 cagcagCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCAAt                                                             >  2:385857/1‑76 (MQ=255)
    cagcTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTa                                                          <  2:560507/76‑1 (MQ=255)
             cACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATa                                                 >  1:253438/1‑76 (MQ=255)
              aCCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAg                                                >  1:117514/1‑76 (MQ=255)
               ccTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGCTAgg                                               >  2:707860/1‑76 (MQ=255)
               ccTGTTCTTCGCTCAGGTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGGGCATCAGCAATTTACCGCCGATAgg                                               <  2:1122427/76‑1 (MQ=255)
                  gTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGgcgc                                            >  1:689626/1‑76 (MQ=255)
                  gTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGgcgc                                            >  2:624454/1‑76 (MQ=255)
                    tcttcGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGGGCATCAGCAATTTACCGCCGATAGGCGCgg                                          <  2:1200811/76‑1 (MQ=255)
                     cttcGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGa                                         >  1:1539404/1‑76 (MQ=255)
                       tcGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGACa                                       >  2:907534/1‑76 (MQ=255)
                              gtggggCAAAAAAGGCTTTCCCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGACATTCGCAt                                <  2:563653/73‑1 (MQ=255)
                                  ggCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGa                            >  2:883901/1‑76 (MQ=255)
                                     tAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGc                         >  2:1426075/1‑76 (MQ=255)
                                        aaaGGCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTa                      >  2:591519/1‑76 (MQ=255)
                                         aaGGCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTac                     >  1:459984/1‑76 (MQ=255)
                                            gCTTTACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACt                  >  2:1374831/1‑76 (MQ=255)
                                              tttACCTGCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGt                >  2:727235/1‑76 (MQ=255)
                                                    tGCACCGGAAGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGACAATCGCATCAGATGCCGCCACCGCACCTgg          >  2:149674/1‑76 (MQ=255)
                                                    cgCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTgg          <  2:117514/75‑1 (MQ=255)
                                                     gCACCGGAAGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGt         >  2:618795/1‑76 (MQ=255)
                                                     gCACCGGAAGCGTGCATCAGCAATTTAACGCCGATAGGCGCGGACATTCGCATCAGATGCGTACCCTGTACCTGGt         >  2:405548/1‑76 (MQ=255)
                                                         cGGAAGCGTGCATCAGCAATTTCCCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGat     <  2:689625/76‑1 (MQ=255)
                                                            aaGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGataat  >  2:1065523/1‑76 (MQ=255)
                                                            aaGCGTGCATCAGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTAGTCGataaa  >  2:890098/1‑75 (MQ=255)
                                                                       |                                                                
GCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACCTGCACCGGAAGCGTGCATCGGCAATTTACCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAAT  >  REL606/4202025‑4202160

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: