Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 44.0% | 23.1 / 10.7 | 25 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (11/3); new base A (0/11); total (11/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.16e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA > NC_000913/357741‑358038 | tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac > 2:219232/1‑150 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa < 2:579069/151‑1 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac < 2:341052/149‑1 (MQ=255) tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACaaa > 2:29825/1‑150 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccg < 1:370891/151‑1 (MQ=255) cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtc > 2:570210/1‑151 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtct > 2:199486/1‑151 (MQ=255) gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 1:169651/1‑150 (MQ=255) ccGCAGCCTTGTGTCAGCGGCTCCCGTGATTCTTTCCCGCCAACACCCGCGC‑TTTTACCCACCGCCTTGTGAAAAAATCCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg < 2:254116/148‑1 (MQ=255) gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg > 2:325109/1‑146 (MQ=255) ccTGTGTCGGCGCCTCCCGTGTTTCATTCCCCCCAACACCCGCCT‑TTTTACCCACCCCCATGTGAAAAAATGCCTGCGGGGCGGCTGTCTGCCAGTTATTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag‑ag < 2:278963/150‑1 (MQ=255) ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAag > 2:215879/1‑148 (MQ=255) gtggcgccgcctccGGGATTCCTTCCCCCCAACAACCGCCCAATTTACCCAACCCCATGCAAAAAAACGCCTGCGGGCCGGCTGTCCCCAACCTATTTGC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAAt < 1:29821/139‑1 (MQ=255) cggCGGCTACCGTTATTCACTCCCCCCAACACCCGCTT‑ATTTACCCACCCCCTTGAGAAAAAACCCCTCCGCGGCGGTTTTCTGCCAGTTTTTTTC‑CCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa < 1:466669/137‑1 (MQ=255) tACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAG‑CGaaaaaaacttaaaa > 1:6814/1‑139 (MQ=25) ccGCCAACACCCGCGT‑TTTTACCCACCCCTTTGTAAAAAATTCCCTCGGCGGCGGTTTTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTtgtgtg < 2:454333/150‑1 (MQ=255) cgcgC‑TTTTACCCACCCCCGTGTCAAAAAATGCCTCCGCGGCGGCTGTCTCCCAGTTATTTCC‑CCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc < 1:53622/150‑1 (MQ=255) cgc‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAG‑AGAAAAAAATTTTAAACCCCCATGTTGTGTGGGGATCgc > 2:355375/1‑145 (MQ=25) ccccATGTAAAAAAAGGCTTTCGGGGCGGCTGTCGGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGaa < 2:214872/149‑1 (MQ=255) gaacAAAATTGCCTCCGCGGGGTCTG‑CCCCCAGTTTTTTTTCCCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTgg < 1:509728/140‑1 (MQ=255) aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAG‑AGAAAAAAATTTTAAACCCCCCTGTTTTGTGGGGATAGCCGC‑TGGGGAAAAAGTGGGAAGGGGAgaggg > 2:37579/1‑147 (MQ=255) ctttcgggcgggcTGTCTCCCTCTTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTgcggc < 2:496042/137‑1 (MQ=255) gcggcgGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc < 2:251938/150‑1 (MQ=255) cggCCGGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc < 2:2382/148‑1 (MQ=255) ccattttttttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc < 2:92946/140‑1 (MQ=255) ttttttc‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt < 1:570091/143‑1 (MQ=255) tATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGGCG‑AGAAATATCTTTAAAACCCCCCTGTTGTGTGGGTATCTCTGCTTGGTACCAA‑GTGGGAAGGGGAGGGAGGGGAGTCAGCGCCCGaaacccccgcatt > 1:286973/1‑138 (MQ=255) t‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGCGGAAGCGGAAAAAAATTAAACCCCCCGATCTT‑TGTGGGAAAAGGCGT‑TGGTGAGAAAGCGGGAAGGGGAGAGCGGGGAATTAACTCCAGAACcgcgggcgct > 1:70982/1‑137 (MQ=255) aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGa < 2:180415/151‑1 (MQ=255) | TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGC‑ATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAG‑AGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGA > NC_000913/357741‑358038 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |