Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 53.8% | 16.2 / 58.1 | 39 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (14/4); new base A (0/21); total (14/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.03e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-05 |
GCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG > NC_000913/357747‑358037 | gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAAccg > 2:243355/1‑151 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAAccgc < 1:1209719/151‑1 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAAccgcc > 1:1014721/1‑151 (MQ=255) aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGt > 2:892951/1‑151 (MQ=255) cAGGCAATTCCAGCCGCAGACCTGTGCCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTtctc > 1:726604/1‑148 (MQ=255) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg > 1:841571/1‑150 (MQ=255) gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:975387/1‑149 (MQ=255) cagccttgtgTCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTGCAAAAATGCTTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa < 1:51395/144‑1 (MQ=21) cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg > 1:84397/1‑147 (MQ=255) gcccTGTGTCAGCGGCTCCCGTGATTCTTTCCCCCCACCACCCGCGCATTTACCCACCGCCATGCCAAAAAATGCCTGCGCGGCGGCTGTCTGCCA‑GCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag < 2:410939/148‑1 (MQ=255) tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcgcggg > 2:278737/1‑143 (MQ=255) gtgtCAGCGGCACCCGTAATTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCAATGTGAAAAATTGCCTTCGCGGCGGCTGTCTGCCA‑GTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑aa < 1:571519/149‑1 (MQ=255) tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggaaggaaa > 2:832191/1‑138 (MQ=255) ggCTCCCGTGTTTCATTCCCCCAAACACCCGCGCTTTTCCCAACCGCCTTGTGCAAAAATCCCTTGGGGGCGGCTGTCTCCCA‑GTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag < 1:542443/139‑1 (MQ=255) gcctccgttatttCTTTCCCCCCAACACCCCCCCTTTTACCCACCCCATTGTGAAAAAATGCCTCCGGGGGGTTTGTCTCCCA‑GTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTaa < 2:231432/140‑1 (MQ=255) ccGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACCCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑TTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGccc < 1:118006/150‑1 (MQ=25) ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAG‑AAAAAAA‑TTTaaaaccccccctttttg > 1:29861/1‑135 (MQ=25) ccGCCAACACCCGCGCATTTACCCACCCCCATGTGAAAAAATGCTTTCGGGGCGGCTGTCTGCCA‑GCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTtgtgtg < 1:706877/150‑1 (MQ=255) gCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAG‑AAAAAGCTTTTAAAACCG‑CCTGTTGTGTGGcg > 1:1096841/1‑150 (MQ=25) ccGCCCTTTTACCCAACGCATTTTGAAAAAATGCTTCCGGGGCGGCTTTCTGCCA‑GTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgc < 1:1146738/148‑1 (MQ=255) gcgcATTTATCCACCGCCATGTGCAAAAATCCTTTCGCGGCGGCTTTTTCCCA‑GCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctgct < 1:128241/150‑1 (MQ=255) aTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAG‑AAAAAAA‑TTTTAAACCCCCCTGTTGTGTGGGGATCGCTGCTgg > 1:879485/1‑148 (MQ=25) ttactccacccctaTGTGAAAAAAGGCCTTCGCGGCGGTTGTCGGCCA‑GTTTTTTTCCGGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg < 1:258390/133‑1 (MQ=255) ccACCGCCTTGTAAAAAAATGCCCCCGCGGCGTGTGTCGCCCA‑CCTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt < 2:408038/151‑1 (MQ=255) acctccttgtgAAAAAAGGCCTCCGCGGGGGCTGTCTCCCG‑GCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg < 1:122420/143‑1 (MQ=255) gtgCAAAAATTCTTCCGCGGGGGTTGTCTGCCA‑GCTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg < 2:192610/149‑1 (MQ=255) tgaaaaatTCCCTTCGCGGCGGCTGTCGCCCA‑GTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa < 1:1200253/147‑1 (MQ=255) aaaaaTGCCTTCGCGGCGGCTGTCTGCCA‑GCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc < 1:775079/151‑1 (MQ=255) aaaaaTCCCTTCGGGGGGGTTGTCTGCCA‑GTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc < 1:125454/151‑1 (MQ=255) aaaaaTCCCTTCGCGGCGGTTGTCGCCCA‑TTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt < 2:1166190/149‑1 (MQ=255) aaaaAGCCCTCCCGGCGGCCTTCCGCCACGCTTTTTCCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa < 2:660123/147‑1 (MQ=255) ccTCCGGGGCGGCTGTCTGCCA‑GTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt < 1:857688/151‑1 (MQ=255) cgcgGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAAAT‑TTTAATCCCG‑CATTTTTTGGGGGCATAGCTGCTGGGTAAAATGGGGGAAGGGGAgggggggggaacaa > 1:943409/1‑137 (MQ=255) tgccgccgattttttttcCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa < 2:886695/135‑1 (MQ=255) ccA‑GTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg < 1:601815/146‑1 (MQ=255) cTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAA‑GAAAAGA‑GTTTAAACCGCCCCGTTGTGGGGCAATAGCCGCGTGGGACAAAATGGGAAGGGGAGAGGGGCGGATCAGCTCCACAAACgcgcccatt > 1:551112/1‑145 (MQ=255) ttttcCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg < 2:1135941/150‑1 (MQ=255) cccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 2:1090671/148‑1 (MQ=255) ccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc < 2:541791/150‑1 (MQ=255) caacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGGG‑AAAAAAC‑TTTTAACCCCCCATGTTTTGTGGGAATTGCTGCTGGGAAAAAAGTGGGAAGGGGAGAGGGGGGGATCAGCCCCCGAAACGCGCCCATTGAGTTTGCCTCCATg > 1:430697/1‑151 (MQ=255) | GCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCA‑GCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG > NC_000913/357747‑358037 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |