Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 330,450 | 0 | A | C | 51.4% | 44.5 / 55.1 | 35 | V266G (GTT→GGT) | cynR | DNA‑binding transcriptional dual regulator |
Reads supporting (aligned to +/- strand): ref base A (14/3); new base C (0/18); total (14/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.93e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT > REL606/330376‑330521 | cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGTCATCTCTTTTTTCGCCG‑CAACAAAAc > 2:2489920/1‑76 (MQ=21) cgcATTTATCCAACGCCATGTGCAAAAAGGCCTCCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAAc < 2:2108078/76‑1 (MQ=255) gcATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTTCGCCG‑CCACaaaaac > 2:831914/1‑74 (MQ=21) cATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAAccg > 1:280455/1‑76 (MQ=255) cATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAAccg > 1:1444991/1‑76 (MQ=255) aTTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CCACAAAAccgc > 2:1717595/1‑76 (MQ=255) tttATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAAccgcc < 1:1542344/76‑1 (MQ=255) tttATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAAccgcc > 1:647504/1‑76 (MQ=255) ttATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGGCAACTATTTTTCCGCAG‑CAACAAAAccgccg > 1:1285755/1‑76 (MQ=25) ttATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTTCGCCG‑CAACAAAAccgccg > 2:300134/1‑76 (MQ=255) tATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGt > 1:2174583/1‑76 (MQ=255) tATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGGTTGCCCGGCATTTTTTCGCCG‑GAAACAAAccgccct > 1:4018268/1‑74 (MQ=12) aTCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTTCGCCG‑CAACAAAACCGCCGtt > 2:5066139/1‑76 (MQ=255) aTCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGtt > 1:3476303/1‑76 (MQ=255) tCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAAACGCCGGtc > 1:3631756/1‑76 (MQ=255) ccAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTtct > 1:837610/1‑76 (MQ=255) ccAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTtct > 1:3059645/1‑76 (MQ=255) aCGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTTCGCCG‑CAACAAAAACACCGTtctctc > 2:3155097/1‑76 (MQ=21) aCGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTtctctc > 1:305359/1‑76 (MQ=255) ccATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCag > 1:1963513/1‑76 (MQ=255) ccATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGGTCTCTcccg > 1:2784601/1‑74 (MQ=255) aTGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCagta > 1:3710554/1‑76 (MQ=255) tgtgCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTTCGCCG‑CAACAAAAACGCCGTTCTCTCCcgtag > 2:5106948/1‑76 (MQ=25) tgtgCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCagtag > 1:4911855/1‑76 (MQ=255) gtgCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCcgtagt > 1:3237885/1‑76 (MQ=255) tcaaaaatTCCCTCCGC‑GGCGGCTGTCTCCCAGCTTTTTTCCCGCCG‑CAACAAACCCGCCGTTCTCTCCAGTAGTg < 1:1385196/73‑1 (MQ=255) tcaaaaaaTCCCTCGGG‑GGCGGCTGTCTGCCAGCTTTTTTCCCGCCG‑CAACAAACCCGCCGTTCTCTCCAGTAGTg < 1:3421377/73‑1 (MQ=255) gcaaaaaaTTCCTTTGCGGGCGTCTCCCTCCGAGTTTTTTTCCCCCCG‑CAACAAACCCGCCGTTCTCTCCagtagt < 2:1861230/73‑1 (MQ=255) aaaaaTGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTTCGCCG‑CAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg > 2:3064163/1‑76 (MQ=255) aaaaaTGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCACTCGTGGCg > 2:1094276/1‑76 (MQ=255) aaatTGCCTGCGG‑GGCGGCTGTCTCCCAGTTTTTTTCCCGCCA‑AAAAAAACCCGCCGTTCTCTCCAGTAGTGGCgg < 1:396737/76‑1 (MQ=255) aaaaTGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTTCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 1:3318118/1‑76 (MQ=255) actGCCTTCGC‑GGGGGCTGTCCGCCACCTTTTTTCCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg < 1:4444512/74‑1 (MQ=255) aaTGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTTCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 2:2140390/1‑76 (MQ=255) aaTGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAAACGCCGTTCTCTCCCGTAGTGGCgggg > 1:1848445/1‑76 (MQ=255) aTGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTTCGCCG‑CAACAAAAACGCCGTTCTCTCCAGTAGTGGCggggg > 2:2200403/1‑75 (MQ=255) tcccTTGGC‑GGCGGCTGTCTCCCATTTTTTTTCCCGCCG‑CAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCa < 1:4086564/74‑1 (MQ=255) gCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggg > 1:3873291/1‑73 (MQ=255) gCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTACAGTAGTGGCggggggg > 1:5049216/1‑73 (MQ=255) tccgc‑gGGGGCTGTCCCCCGTTTTTTTTCCCGCCG‑CAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag < 1:4027062/74‑1 (MQ=255) tCGG‑GGGGGCTGTCTGCCAGTTTTTTTCCCGCCG‑CAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga < 1:2373704/76‑1 (MQ=255) ggc‑ggcgCCTGTCTGCCATCTTTTTTCCCGCCG‑CAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa < 1:3954169/75‑1 (MQ=255) gc‑ggcgGCTGTCGCCCATCTTTTTTCCCCCCG‑CAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa < 1:3113697/76‑1 (MQ=255) c‑ggcggCTTTCGGCCAGTTTTTTTCCCGCCG‑CAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt < 1:1097858/76‑1 (MQ=255) gcgccTGTCTCCCAGTTTTTTCCCCCCCG‑CAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg < 1:2672040/76‑1 (MQ=255) cTGTCGCCCAGTTTTTTCCCCGCCG‑AAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt < 1:2479289/76‑1 (MQ=255) tctctcCCAGCTTTTTCCCCGCCG‑CAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa < 1:2189948/74‑1 (MQ=255) tGTCTGCCACTTTTTTTCCCGCCG‑CAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa < 1:4886243/76‑1 (MQ=255) tctcCCAGCTTTTTTCCCCCCG‑CAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg < 1:4878608/76‑1 (MQ=255) ctcTATTTACCCCCCG‑AAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc < 1:3593727/74‑1 (MQ=255) ttttttcCCGCCG‑CAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg < 1:90326/74‑1 (MQ=255) tttttcCCGCCGAAAACCAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTGAtg < 1:3810758/75‑1 (MQ=255) tttttcCCCCCG‑CAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt < 1:2037334/75‑1 (MQ=255) aTTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGGGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt > 1:530004/1‑76 (MQ=255) ttcccgccg‑AAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt < 1:2356483/73‑1 (MQ=255) cccccgcca‑CAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt < 2:2046182/73‑1 (MQ=255) cggccg‑cAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg < 1:3118291/74‑1 (MQ=255) ccgccg‑CAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg < 1:1744640/76‑1 (MQ=255) cccg‑aAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGc < 1:1952444/75‑1 (MQ=255) ccg‑cAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCa < 1:1050889/76‑1 (MQ=255) cg‑cAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGGGGAAAAGAGTTAAAACCGTCACGTTGTGTgggga > 1:3414133/1‑73 (MQ=255) ccacaaaaaACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCa < 1:1412657/73‑1 (MQ=255) g‑cAACAACCCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt < 1:4585974/76‑1 (MQ=255) aaaacaacCCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt < 2:1723276/74‑1 (MQ=255) caaACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgct < 1:3691699/76‑1 (MQ=255) | CGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGC‑GGCGGCTGTCTGCCAGCTATTTTTCCGCCG‑CAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCT > REL606/330376‑330521 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |