Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,202,405 +A coding (155/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,4051.A93.8% 173.3 / 3.8 48coding (155/825 nt)iclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base A (24/21);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCA  >  REL606/4202330‑4202473
                                                                            |                                                                    
ccGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCtt                                                                       <  2:897983/76‑1 (MQ=255)
ccGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCtt                                                                       <  1:583580/76‑1 (MQ=255)
 cGACCTGACGCACGCAACCCTGCTATTGCATNGTGGTTAGCAGGCGNTGGGTNGTGGAATTGGGTAGCCCGGCtta                                                                      <  2:180974/76‑2 (MQ=255)
   aCCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCttagt                                                                    >  2:1121833/1‑73 (MQ=255)
      tGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTgc                                                                 >  1:914381/1‑76 (MQ=255)
      tGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTgc                                                                 >  1:659383/1‑76 (MQ=255)
        aCGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTgcgc                                                               <  1:1229670/76‑1 (MQ=255)
           cACGAAACCCTGCTGTTGCATTCTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAg                                                            >  1:7730/1‑76 (MQ=255)
             cGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTc                                                         >  1:1027169/1‑76 (MQ=255)
               aaaCCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTcc                                                        <  1:747485/76‑1 (MQ=255)
                 aCCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGt                                                      <  1:1067833/76‑1 (MQ=255)
                  cccTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTg                                                     <  2:1004671/76‑1 (MQ=255)
                    ctgctgTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAg                                                   >  2:952819/1‑76 (MQ=255)
                       ctgTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGc                                                >  1:1003720/1‑76 (MQ=255)
                          ttGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCcac                                             <  1:919154/76‑1 (MQ=255)
                          atacATTGTGGTTNGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCcac                                             <  1:181090/73‑1 (MQ=255)
                           tGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCcaca                                            >  2:595418/1‑76 (MQ=255)
                            gCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCacac                                           <  2:845088/76‑1 (MQ=255)
                             cATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACt                                          <  1:702334/76‑1 (MQ=255)
                              aTTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTg                                         <  1:1200194/76‑1 (MQ=255)
                              aTTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTg                                         <  2:417303/76‑1 (MQ=255)
                               ttGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGc                                        >  1:231066/1‑76 (MQ=255)
                                   ggTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCAAtt                                    <  1:1107646/76‑1 (MQ=255)
                                        gCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGAtt                               <  2:1104245/76‑1 (MQ=255)
                                        gCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGAtt                               <  2:581530/76‑1 (MQ=255)
                                          aGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCg                             >  1:726723/1‑76 (MQ=255)
                                          aGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCg                             >  2:430238/1‑76 (MQ=255)
                                          aGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCGAGTTCCGTGAGTGCCACACTGCCATTGGATTCgg                            >  2:133085/1‑76 (MQ=255)
                                            gCGGTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGc                           <  1:1098070/76‑1 (MQ=255)
                                               gTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAAt                        >  1:1019641/1‑76 (MQ=255)
                                               gTGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAAt                        <  1:209274/76‑1 (MQ=255)
                                                tGGGTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATc                       >  2:870677/1‑76 (MQ=255)
                                                  ggTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCa                     >  1:827190/1‑76 (MQ=255)
                                                  ggTCGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGAATTCGGCAATCCa                     <  1:93004/76‑1 (MQ=255)
                                                     cGTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTc                  >  2:977807/1‑76 (MQ=255)
                                                      gTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTcc                 >  2:1064201/1‑76 (MQ=255)
                                                      gTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTcc                 >  1:1259078/1‑76 (MQ=255)
                                                      gTGGAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTCCCATTGGATTNGGCAATCCAccct                 >  1:186058/1‑73 (MQ=255)
                                                       tGGAATTGGGTAACCCGGNTTAGTTGCGCCAGTTCCGTGAGTGCCAGACTGCCATTGGATTCGGCAATCCACTCCa                >  1:1351843/1‑76 (MQ=255)
                                                        ggAATTGGGTAACCCGGCTTAGCTGCGCCAGTTCCGTGAGTGCCACACTGCCATCGGATTCTGCAATCCACTCCAg               >  2:206903/1‑76 (MQ=255)
                                                         gAATTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGt              <  2:765347/76‑1 (MQ=255)
                                                         gAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCNACTCCAGTa             <  1:37858/76‑1 (MQ=255)
                                                          aaTTGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTa             >  2:757765/1‑76 (MQ=255)
                                                            ttGGGTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAAt           >  1:1014564/1‑76 (MQ=255)
                                                              gggTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAAttt         <  1:972959/76‑1 (MQ=255)
                                                                gTAACCCGGCTTAGTTGTGCCNCTTCCNTGAGTGCCACACTGCCATTGNATACGGCAATCTACCCAAGTACTTTCa       >  1:1378157/1‑76 (MQ=255)
                                                                gTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCa       >  2:258818/1‑76 (MQ=255)
                                                                gTAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCa       <  1:755100/76‑1 (MQ=255)
                                                                 tAACCCGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAg      >  2:684391/1‑76 (MQ=255)
                                                                  aaCCCGGCTTAGTTGCGCCAGTTCCGTGAGAGCCACACTGCCATTGGATTCGGCAATCCACTNCAGTAATTTCAgg     >  2:1294574/1‑76 (MQ=255)
                                                                    cccGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGcc   <  2:827190/76‑1 (MQ=255)
                                                                     ccGGCTTAGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCa  <  2:1065492/76‑1 (MQ=255)
                                                                            |                                                                    
CCGACCTGACGCACGAAACCCTGCTGTTGCATTGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCA  >  REL606/4202330‑4202473

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