Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,874,154 (GT)5→4 coding (513‑514/936 nt) lysR → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,874,1460G.94.1% 142.3 / 4.6 34coding (505/936 nt)lysRDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base . (23/9);  total (23/11)
Fisher's exact test for biased strand distribution p-value = 9.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
*REL6062,874,1470T.96.9% 136.0 / 2.3 34coding (506/936 nt)lysRDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base . (22/9);  minor base A (0/1);  total (23/11)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CATGATTTAGGACTCACTGAAACGCTCCATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAAGTGTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCGGATGATTTTCAGGGTG  >  REL606/2874071‑2874215
                                                                           ||                                                                    
cATGANTCAGGACTCACTGAAACGCTCCTTACGCCTGCGGGAACAGAACGTACCGAATTACNCTCTTTAGATGAAg                                                                       >  2:1718828/1‑76 (MQ=255)
   gATTTAGGACTCACTGAAACGCTCCATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAAgtgt                                                                    >  1:875340/1‑76 (MQ=255)
     tttAGGACTCACTGAAACGCTCCATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAAgtgtgt                                                                  >  1:1212034/1‑76 (MQ=255)
        aGGACTCACTGAAACGCTCCATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAAgtgtgtgtt                                                               <  1:1519190/76‑2 (MQ=255)
         ggACTCACTGAAACGCTCCATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAAgtgtgtgtta                                                              <  2:1611315/76‑3 (MQ=255)
         ggACTCACTGAAACGCTCCATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAAgtgtgtgtta                                                              <  2:1420575/76‑3 (MQ=255)
                  gAAACGCTCCATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGt                                                   >  2:543304/1‑76 (MQ=255)
                  gAAACGCTCCATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGt                                                   >  2:396908/1‑76 (MQ=255)
                   aaaCGCTCCATACGCCTGCGGGAACAGAACGTACCGAATTGCTCTCTTTAGATGA‑‑CGTGTGTGTTACCTCCGGGTc                                                  >  2:1363265/1‑76 (MQ=255)
                      cGCTCCATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATc                                               >  2:164333/1‑76 (MQ=255)
                       gCTCCATACGCCTGCGGGAANAGAACGTACAGAATTNCTCTCTTTAGATG‑‑GACTGTGTCTTACCTCCGGGTCATcc                                              <  2:1698907/76‑1 (MQ=255)
                       gCTCCATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATcc                                              >  1:502793/1‑76 (MQ=255)
                       gCTCCATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATcc                                              >  1:502796/1‑76 (MQ=255)
                           cATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTg                                          >  2:615336/1‑76 (MQ=255)
                            aTACGNCGGAGGGAACAGAACGTACCGACATACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCNGGGTCATCCGCTgg                                         >  2:1711019/1‑76 (MQ=255)
                            aTACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTgg                                         >  1:476932/1‑76 (MQ=255)
                            aTACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTgg                                         >  2:1231784/1‑76 (MQ=255)
                                 ccTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTa                                    <  1:1295386/76‑1 (MQ=255)
                                         aaCAGAACGTACCGNATTACTCTCGTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGt                            >  1:8268/1‑76 (MQ=255)
                                         aaCAGAACGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGt                            >  1:118230/1‑76 (MQ=255)
                                         aaCAGAACGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGt                            <  2:1293642/76‑1 (MQ=255)
                                           cAGAACGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTAt                          >  1:1524756/1‑76 (MQ=255)
                                              aaCGTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTaa                       <  1:796892/76‑1 (MQ=255)
                                                 gTACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGc                    >  2:124127/1‑76 (MQ=255)
                                                  tACCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGcc                   >  1:75456/1‑76 (MQ=255)
                                                   aCCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCg                  >  2:1029384/1‑76 (MQ=255)
                                                   aCCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCg                  >  1:938966/1‑76 (MQ=255)
                                                   aCCGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCg                  >  1:813241/1‑76 (MQ=255)
                                                    ccGAATTACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCgg                 <  1:277771/76‑1 (MQ=255)
                                                     cGAATTACTCTCTTTAGATGAAGTGTGGGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCGga                <  1:173649/76‑1 (MQ=255)
                                                       aaTTACTCTCTTTANATGAA‑‑GTGTGTGTTACCTCCGGGTCATACGCTGGCGGTAAAAAAGTTAGTAACGCCGgatg              <  1:1466549/76‑1 (MQ=255)
                                                          tACTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCGGATGAtt           >  2:721882/1‑76 (MQ=255)
                                                           aCTCTCTTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGAATTAACGCCGGATGAttt          >  2:52671/1‑76 (MQ=255)
                                                             tctctTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCGGATGATTTTc        <  2:1145979/76‑1 (MQ=255)
                                                                ctTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCGGATGATTTTCAgg     <  1:1095913/76‑1 (MQ=255)
                                                                ctTTAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCGGATGATTTTCAgg     <  1:1518208/76‑1 (MQ=255)
                                                                 tttAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCGGATGATTTTCAggg    >  2:1012426/1‑76 (MQ=255)
                                                                 tttAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCGGATGATTTTCAggg    >  2:641908/1‑76 (MQ=255)
                                                                  ttAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCGGATGATTTTCAGGGt   <  1:156901/76‑1 (MQ=255)
                                                                   tAGATGAA‑‑GTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCGGATGATTTTCAGGGTg  >  2:735510/1‑76 (MQ=255)
                                                                           ||                                                                    
CATGATTTAGGACTCACTGAAACGCTCCATACGCCTGCGGGAACAGAACGTACCGAATTACTCTCTTTAGATGAAGTGTGTGTGTTACCTCCGGGTCATCCGCTGGCGGTAAAAAAGGTATTAACGCCGGATGATTTTCAGGGTG  >  REL606/2874071‑2874215

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