Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,631,772 | 0 | A | C | 37.8% | 77.9 / 13.5 | 45 | intergenic (‑267/‑141) | rpoE/nadB | RNA polymerase, sigma 24 (sigma E) factor/L‑aspartate oxidase |
Reads supporting (aligned to +/- strand): ref base A (24/4); new base C (0/17); total (24/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.43e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAG > REL606/2631698‑2631845 | aTCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccac < 1:289717/76‑1 (MQ=255) cTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTACCCCAACGGCTTTTTTTTTTCCCCCCCTAACCCCccccca < 2:184343/76‑1 (MQ=255) cTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccacca < 1:1407523/76‑1 (MQ=255) tgtgTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag < 2:551482/76‑1 (MQ=255) tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc < 2:703860/76‑1 (MQ=255) tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc < 1:14231/76‑1 (MQ=255) aGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGt > 1:1405067/1‑76 (MQ=255) aGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGt > 2:1577391/1‑76 (MQ=255) aGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGt > 2:1577376/1‑76 (MQ=255) aaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAAc > 1:883742/1‑76 (MQ=255) tGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCCCCTAATCCTCCCCCAGCCAGTAACtt < 2:425586/76‑1 (MQ=255) gCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCCCCTAATCCTCCCCCAGCCAGTAACTtc < 2:37195/76‑1 (MQ=255) cATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTctct > 1:593883/1‑76 (MQ=255) nttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCttt > 1:92787/2‑76 (MQ=255) ttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtttt > 1:600644/1‑76 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCCCCAGCCAGTAACTTCTCttttt < 2:1071797/76‑1 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCttttt > 2:669004/1‑76 (MQ=255) accTGCCTGAAGAGTAACCCAACGGCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTtc < 1:556871/74‑1 (MQ=255) aaCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtc > 2:892413/1‑76 (MQ=255) gAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCgccgcc > 2:1344257/1‑76 (MQ=255) gAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCgccgcc > 2:1102279/1‑76 (MQ=255) aaGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGccc > 2:1224148/1‑76 (MQ=255) gagTACCCCAACGGCCTTTTTATTTCACCCCCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTg < 2:434976/76‑1 (MQ=255) agTACCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGc < 1:431463/76‑1 (MQ=255) accccACCGGCTTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc < 1:405285/74‑1 (MQ=255) aaCCCAACGGCCTTTTTATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc < 1:226478/76‑1 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc > 2:1082339/1‑76 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc > 2:1283816/1‑76 (MQ=255) cccACCGGCCTTTTTATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAg < 2:1529353/76‑1 (MQ=255) cAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCg > 1:1238670/1‑76 (MQ=255) cGGCTTTTTTATTTCCCCCCCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgtgt < 2:1488297/76‑1 (MQ=255) cGGCCTTTTTATTTCACCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgtgt < 2:530734/76‑1 (MQ=255) gCCTTTTTATTTCCCCACCTAACCCCCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt < 1:434050/76‑1 (MQ=255) gCCTTTTTATTCCACCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTNGCCGCCCTGCGTCAGCGTGttt < 2:234451/76‑1 (MQ=255) ccTTTTTATTTCCCCCCCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 2:1591459/76‑1 (MQ=255) ccTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa < 2:725541/76‑1 (MQ=255) tttATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCa < 1:434755/76‑1 (MQ=255) tttATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCa < 1:1210456/76‑1 (MQ=255) tttATTCCCCCCCCTAACCCCCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCa < 2:1680097/76‑1 (MQ=255) ttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCAc > 1:977077/1‑76 (MQ=255) ttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCAc > 2:1099408/1‑76 (MQ=255) aTTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTg > 2:1255427/1‑76 (MQ=255) tttCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGt < 1:1673024/76‑1 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa > 2:34488/1‑76 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa > 2:1486394/1‑76 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa > 2:1480094/1‑76 (MQ=255) aTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaa > 2:188470/1‑76 (MQ=255) ctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAAt > 1:1068518/1‑76 (MQ=255) tccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATa > 1:269086/1‑76 (MQ=255) ccaccaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAg > 1:422524/1‑76 (MQ=255) | ATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAG > REL606/2631698‑2631845 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |