Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,202,529 +CT coding (31/825 nt) iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,202,5261.T89.3% 99.7 / 5.8 28coding (34/825 nt)iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base T (0/25);  total (1/27)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01
*REL6064,202,5262.C89.3% 102.7 / 6.6 28coding (34/825 nt)iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (0/25);  total (1/27)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

TCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTT‑‑TCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA  >  REL606/4202456‑4202586
                                                                         ||                                                            
tCCAGTAATTTCAGGCCGCGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAGCGGCGGG‑‑TTT‑‑TCtct                                                         <  1:836885/76‑3 (MQ=255)
tCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTT‑‑TCtct                                                         <  2:8403/76‑3 (MQ=255)
 ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTT‑‑TCtctg                                                        <  2:72382/76‑4 (MQ=255)
    gTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTT‑‑TCTGCcg                                                       >  1:1457135/1‑76 (MQ=255)
          tCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTc                                                 <  1:1547618/76‑1 (MQ=255)
          tCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTc                                                 <  2:1164655/76‑1 (MQ=255)
               ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCggg                                            <  2:1337958/76‑1 (MQ=255)
                  cgcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAAt                                         <  2:454434/76‑1 (MQ=255)
                   gcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATg                                        <  2:1600675/76‑1 (MQ=255)
                     gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑GTT‑‑TCTTCCGCGGTGCGCGGGGGTggggg                                    >  1:309185/1‑74 (MQ=255)
                      ttAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTT‑‑TCTGCCGCGCTTCGCGGGGGTgggg                                     >  1:735403/1‑75 (MQ=255)
                       tAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTT‑‑TCTTCCGCGGTTCGCGGGGATggggg                                    >  2:680597/1‑74 (MQ=255)
                       tAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTg                                    <  1:81267/76‑1 (MQ=255)
                        aaaGACTGACCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGc                                   <  2:435899/76‑1 (MQ=255)
                        aaaGACTGAACCTGTCCAGTCGCTGGTTCGGTGGCGACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGc                                   <  2:299359/76‑1 (MQ=255)
                        aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTT‑‑TCTGCCGCGTTTCGCGGGGATgggggg                                   >  2:1359977/1‑73 (MQ=255)
                         aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCg                                  <  1:1609624/76‑1 (MQ=255)
                          aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGa                                 <  2:335662/76‑1 (MQ=255)
                            aCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGAcc                               <  2:9388/76‑1 (MQ=255)
                              tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTT‑‑TCTGCCGCGGTTCGCGGGGATGGGGGGGGCCAt                             >  2:1644937/1‑76 (MQ=255)
                               gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATg                            <  2:111435/76‑1 (MQ=255)
                                aaCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCTGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGa                           <  2:243988/76‑1 (MQ=255)
                                aaCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGa                           <  2:1155798/76‑1 (MQ=255)
                                 aCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAc                          <  2:1299892/76‑1 (MQ=255)
                                    tGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTT‑‑TCTGCCGCGTTTCGCGGGGATGGGGGGGGCCATGACAgg                       >  2:600957/1‑75 (MQ=255)
                                          gTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCtt                 <  2:1253227/76‑1 (MQ=255)
                                          gTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCtt                 <  1:1457127/76‑1 (MQ=255)
                                            cGCTGGTGCGGTGGCAACGGCGGGGTTTT‑‑TCTGCCGCGGTTCGCGGGGGTGGGGGGGAGCATGAGCGTCTTCtttt               >  2:862721/1‑76 (MQ=255)
                                               tGGTGCGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCt            <  2:1634338/76‑1 (MQ=255)
                                                    cGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATc       <  1:1446263/76‑1 (MQ=255)
                                                    cGGTGGCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATc       <  1:1594726/76‑1 (MQ=255)
                                                        ggCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTgg   <  2:334154/76‑1 (MQ=255)
                                                        ggCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTgg   <  2:1150611/76‑1 (MQ=255)
                                                        ggCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTgg   <  1:284341/76‑1 (MQ=255)
                                                         gCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGa  <  2:1125821/76‑1 (MQ=255)
                                                         gCAACGGCGGG‑‑TTTTCTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGa  <  1:696868/76‑1 (MQ=255)
                                                                         ||                                                            
TCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑‑TTT‑‑TCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA  >  REL606/4202456‑4202586

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