Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 380,012 | 1 | . | G | 66.7% | 45.5 / 21.5 | 21 | intergenic (‑131/+57) | frmR/yaiO | regulator protein that represses frmRAB operon/outer membrane protein |
Reads supporting (aligned to +/- strand): ref base . (5/2); new base G (9/5); total (14/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01 |
GGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGA‑GAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTC > NC_000913/379868‑380153 | ggAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGA‑GAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑g‑ggggg > 2:30390/1‑151 (MQ=255) gTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGA‑GAATAAAATTTATCCGGTGAATGTGGTCGGAAAAAAAAGAGGAAA‑‑‑g‑ggagat > 2:599990/1‑145 (MQ=255) tACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGA‑GAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑g‑gg < 1:600026/144‑1 (MQ=255) gcatgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGA‑GAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGg‑ccgcc < 2:399347/147‑1 (MQ=255) gcagacaATATGCTAAGGGGGGTTTTTCATTGTTTCAATGCATCCCCCCATAAGTTAAAAAACA‑TGGAGGAAAAAAATTTACCCGGTGAAGGTGGTGGGAAAACAAAGAGGAAAGGGG‑GGGGGGGGG‑‑CTAATCGGCAGGGAAgg < 2:755938/139‑1 (MQ=255) atatACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGA‑GAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGG‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCCCCCCGGATAgcg > 2:751430/1‑151 (MQ=255) aGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGA‑GAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGa > 1:167649/1‑150 (MQ=255) aaTGCATACCCCCCTATAGTATATAACACTGGA‑GAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGAAAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTc > 2:90174/1‑151 (MQ=255) aTGCATACCCCCCTATAGTATATAACACTGGA‑GAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑G‑GGGGGGGGGGCCTATTCGGCGGGGAAGG‑CCGCCCCGGAAAGCGGGGGGCAGAAGGAATCAGAATTTCCgggc > 1:373213/1‑149 (MQ=255) atataACACTGGG‑GAATAAAAGTTATCCGGTGAATGTGGGCGGAAAACAAGGAGGAAAG‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCa < 2:319208/143‑1 (MQ=255) cacTGGA‑GAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑G‑GGGGGGGGGT‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGAAGACCGtt > 2:532481/1‑151 (MQ=255) tGGA‑GAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTaa > 2:756543/1‑150 (MQ=255) tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑G‑GGGGGGGGG‑‑CTAAACGGCAGGGGAGGCCCGCCCGGGATAGGGGGGGGAAGAAGAAAT‑AGGATTTCCAGGGCAAGAGGGGTTG‑AAGGTTCAAACCGTTAAAaaca > 1:347880/1‑148 (MQ=255) aaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCCGATAGCGGGGGGGAGAAGGAAACAGAATTTCCAGGGC‑AGGCGGGCCGGAAGGTTCAGACCGGTAAAATCATCgggt > 2:551891/1‑149 (MQ=255) tATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg < 2:682850/150‑1 (MQ=255) ggTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCg < 1:586103/150‑1 (MQ=255) gAATGTGGTCGGAAAACAAAGAGGAAAGG‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGt > 1:245428/1‑150 (MQ=255) gTCGGAAAACAAAGAGGAAAG‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGCAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACtt > 1:356537/1‑151 (MQ=255) gAAAACAAAGAGGAAA‑‑‑GAGGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAcc > 2:100410/1‑151 (MQ=255) aaaCAAAGAGGAAAG‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGGGGGGGGCAGAAGGAATCAAAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGGGTGGTACCACACTTTAc > 2:234013/1‑148 (MQ=255) aaCAAAGAGGAAAG‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGc > 1:159023/1‑151 (MQ=255) aaaGAGGAAAG‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCg < 2:698982/150‑1 (MQ=255) aaaGAGAAAA‑‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTCCCACACTTTACCTGCCGt < 2:504583/150‑1 (MQ=255) aaGAGGAAA‑‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTc > 1:215547/1‑150 (MQ=255) | GGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGA‑GAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑G‑GGGGGGGGG‑‑CTAATCGGCAGGGAAGG‑CCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTC‑AGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTC > NC_000913/379868‑380153 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |