Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 1,286,699 | C→A | 93.1% | intergenic (+23/‑160) | narI → / → ECB_01206 | nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 1,286,699 | 0 | C→A | 93.1% | 3.6 ‑2.0 ‑0.2 | 49 | intergenic (+23/‑160) | narI/ECB_01206 | nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein |
Reads supporting (aligned to +/- strand): new base (15/30): ref base (3/1): total (18/31) | ||||||||||
Fisher's exact test strand distribution p-value = 1.34e-01 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.32e-01 |
CTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGT > REL606/1286666‑1286731
|
ctagTCACTAAGCGAATTTTAGTTCACATAGACACt < 1:2089961/33‑1 (MQ=38)
cTCGTCACTAAGCGAATTTTAGTTCACATAGACACt > 1:3848927/1‑36 (MQ=255)
cTCGTCACTAAGCGAATTTTAGTTCACATAGACACt < 1:1746794/36‑1 (MQ=255)
gTCACTAAGCGAATTTTAGTTCACATAGACACTGCt < 1:5040305/36‑1 (MQ=255)
gTCACTAAGCGAATTTTAGTTCACATAGACACTGCt < 1:837865/36‑1 (MQ=255)
gTCACTAAGCGAATTTTAGTTCACATAGACACTGCt < 1:2971519/36‑1 (MQ=255)
cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc < 1:5892558/36‑1 (MQ=255)
cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc < 1:6743314/36‑1 (MQ=255)
cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc < 1:3136583/36‑1 (MQ=255)
cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc < 1:1693917/36‑1 (MQ=255)
aCTAAGTGAATTTTAGTTCACATAGACACTGCTTcg < 1:6869881/36‑1 (MQ=38)
cTAAGCGAATTTTAGTTCACATAGACACTGCTTcgg < 1:3689718/36‑1 (MQ=255)
tAAGCGAATTTTAGTTCACATAGACACTGCTTcggc > 1:4259560/1‑36 (MQ=255)
tAAGCGAATTTTAGTTCACATAGACACTGCTTcggc < 1:7922416/36‑1 (MQ=255)
aaGCGAATTTTAGTTCACATAGACACTGCTTcggcg < 1:1817074/36‑1 (MQ=255)
aaGCGAATTTTAGTTCACATAGACACTGCTTcggcg > 1:6208847/1‑36 (MQ=255)
aGCGAATTTTAGTTCACATAGACCCTGCTTcggcgg > 1:4259022/1‑36 (MQ=255)
gCGAATTTTAGTTCACATAGACACTGCTTCGGCggg > 1:1468870/1‑36 (MQ=255)
gCGAATTTTAGTTCACATAGACACTGCTTCGGCggg < 1:6530717/36‑1 (MQ=255)
cGAATTTTAGTTCACATAGACCCTGCTTCGGCgggg < 1:4898156/36‑1 (MQ=255)
gAATTTTAGTTCACATAGACACTGCTTCGGCGGGGt < 1:2517595/36‑1 (MQ=255)
aaTTTTAGTTCACATAGACACTGCTTCGGCGGGGtt < 1:7665839/36‑1 (MQ=255)
aaTTTTAGTTCACATAGACACTGCTTCGGCGGGGtt < 1:4755305/36‑1 (MQ=255)
aTTTTAGTTCACATAGACACTGCTTCGGCGGGGttt < 1:2952910/36‑1 (MQ=255)
aTTTTAGTTCACATAGACACTGCTTCGGCGGGGttt > 1:1082475/1‑36 (MQ=255)
tttAGTTCACATAGACCCTGCTTCGGCGGGGttttt > 1:3178623/1‑36 (MQ=255)
tttAGTTCACATAGACACTGCTTCGGCGGGGttttt < 1:7409515/36‑1 (MQ=255)
ttAGTTCACATAGACCCTGCTTCGGCGGGGtttttt > 1:3950724/1‑36 (MQ=255)
ttAGTTCACATAGACACTGCTTCGGCGGGGtttttt > 1:7448567/1‑36 (MQ=255)
tAGTTCACATAGACACTGCTTCGGCGGGGttttttt < 1:8157667/36‑1 (MQ=255)
tAGTTCACATAGACACTGCTTCGGCGGGGttttttt < 1:1068933/36‑1 (MQ=255)
ttCACATAGACACTGCTTCGGCGGGGTTTTTTTAtt > 1:1378409/1‑35 (MQ=255)
ttCACATAGACACTGCTTCGGCGGGGTTTTTTTAtt > 1:4642612/1‑35 (MQ=255)
ttCACATAGACACTGCTGCGGCGGGGTTTTTTTATg < 1:4299734/36‑1 (MQ=37)
tCACATAGACACTGCTTCGGCGGGGTTTTTTTATgg < 1:2489385/36‑1 (MQ=255)
cacaTAGACACTGCTTCGGCGGGGTTTTTTTATggg < 1:960377/36‑1 (MQ=255)
acaTAGACACTGCTTCGGCGGGGTTTTTTTATTGGc > 1:1421669/1‑36 (MQ=25)
acaTAGACACTGCTTCGGCGGGGTTTTTTTATGGGc < 1:4966621/36‑1 (MQ=255)
caTAGACACTGCTTCGGCGGGGTTTTTTTATTGGCa > 1:7308201/1‑36 (MQ=25)
caTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCa > 1:481629/1‑36 (MQ=255)
caTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCa < 1:7870263/36‑1 (MQ=255)
aTAGACACTGCTTCGGCGGGGTTTTTTTATTTGCAc > 1:8169907/1‑36 (MQ=18)
aTAGACACTGCTTCGGCGGGGTTTTTTTATTGGAAc > 1:2613685/1‑36 (MQ=18)
tAGACACTGCTTCGGCGGGGTTTTTTTTTGGGCagg > 1:2482278/1‑34 (MQ=21)
tAGACACTGCTTCGGCGGGGTTTTTTTATTGGCACg > 1:7961363/1‑36 (MQ=25)
gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt < 1:350876/36‑1 (MQ=255)
gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt < 1:6212519/36‑1 (MQ=255)
gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt < 1:4831082/36‑1 (MQ=255)
gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt < 1:5604073/36‑1 (MQ=255)
|
CTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGT > REL606/1286666‑1286731
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 17 ≤
ATCG
< 26 ≤
ATCG
< 32 ≤
ATCG
< 37 ≤
ATCG
|