Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,414,644 A→G Y207C (TAC→TGC)  crp → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,414,6440AG97.6% 165.1 / ‑4.3 42Y207C (TAC→TGC) crpDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (27/14);  total (27/15)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGC  >  REL606/3414572‑3414714
                                                                        |                                                                      
ccGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCgg                                                                     >  2:1469687/1‑76 (MQ=255)
  gTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCa                                                                   >  1:51281/1‑76 (MQ=255)
   tGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCAc                                                                  >  1:855302/1‑76 (MQ=255)
     ggACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTc                                                                <  1:111675/76‑1 (MQ=255)
           aTTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAAt                                                          >  1:259/1‑76 (MQ=255)
           aTTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAAt                                                          >  2:455004/1‑76 (MQ=255)
           aTTCTGAAGATGCTGGAAGATCAGAACCTGATCTCAGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAAt                                                          >  2:1422041/1‑76 (MQ=255)
             tCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATcc                                                        <  1:94050/76‑1 (MQ=255)
             tCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGGAAAACCATCGTCGTTTGCGGCACTCGTTAATcc                                                        >  1:148708/1‑76 (MQ=255)
               tGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCg                                                      >  2:1358369/1‑76 (MQ=255)
                 aaGGTGCTGGAAGATCAGATCCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTc                                                    <  2:310480/76‑1 (MQ=255)
                  aGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAGCCATCGTCGTTTACGGCACTCGTTAATCCCGTCg                                                   <  2:1295173/76‑1 (MQ=255)
                  aGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCg                                                   >  1:303558/1‑76 (MQ=255)
                         ggAAGATCAGAACNTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGgc                                            >  1:1080771/1‑76 (MQ=255)
                         ggAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGgc                                            >  2:1527961/1‑76 (MQ=255)
                           aaGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGgcgc                                          >  2:566072/1‑76 (MQ=255)
                           aaGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGgcgc                                          >  2:130159/1‑76 (MQ=255)
                             gATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGt                                        <  1:1447105/76‑1 (MQ=255)
                                cAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTAc                                     >  1:1323562/1‑76 (MQ=255)
                                cAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTAc                                     <  2:565055/76‑1 (MQ=255)
                                   aaCCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTg                                  <  2:303558/76‑1 (MQ=255)
                                   aaCCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTg                                  <  2:543058/76‑1 (MQ=255)
                                   aaCCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTg                                  <  2:1355766/76‑1 (MQ=255)
                                   aaCCTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTg                                  <  1:1259106/76‑1 (MQ=255)
                                     ccTGATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGt                                <  1:1572335/76‑1 (MQ=255)
                                        gATCTCCGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCCCTACCTGGTAgc                             >  1:1297125/1‑76 (MQ=255)
                                             ccGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGcc                        >  2:1503202/1‑76 (MQ=255)
                                             ccGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGcc                        >  2:1550407/1‑76 (MQ=255)
                                             ccGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGGAGCGCGcc                        >  2:987443/1‑76 (MQ=255)
                                              cGCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCa                       >  2:422367/1‑76 (MQ=255)
                                               gCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAt                      >  1:252124/1‑76 (MQ=255)
                                               gCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAt                      >  2:1528429/1‑76 (MQ=255)
                                               gCACACGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAt                      >  2:643545/1‑76 (MQ=255)
                                                cacaCGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAtt                     <  1:668911/76‑1 (MQ=255)
                                                cacaCGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCCCGCCCtt                     >  1:669488/1‑76 (MQ=255)
                                                 acacGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAttt                    >  2:459480/1‑76 (MQ=255)
                                                  cacGGTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCAtttt                   <  1:1065443/76‑1 (MQ=255)
                                                     ggTAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGtt                <  2:1158505/76‑1 (MQ=255)
                                                       tAAAACCATCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGGGGGGCCTTTTGTTTc              >  2:518397/1‑76 (MQ=255)
                                                              aTCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGAt       <  2:1529197/76‑1 (MQ=255)
                                                              aTCGTCGTTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGAt       >  2:507621/1‑76 (MQ=255)
                                                                   cgtTTGCGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGCAGCGCGCCCTTTTGTTTCCCCCGATGGGgc  >  1:1263395/1‑76 (MQ=255)
                                                                        |                                                                      
CCGTGGGACGCATTCTGAAGATGCTGGAAGATCAGAACCTGATCTCCGCACACGGTAAAACCATCGTCGTTTACGGCACTCGTTAATCCCGTCGGAGTGGCGCGTTACCTGGTAGCGCGCCATTTTGTTTCCCCCGATGTGGC  >  REL606/3414572‑3414714

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