Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 42.0% | 74.5 / 27.2 | 50 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (23/6); new base A (0/21); total (23/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.11e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA > NC_000913/357741‑358034 | tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac > 2:318637/1‑150 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACaaa < 1:968790/151‑1 (MQ=255) tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACaaa > 1:127642/1‑150 (MQ=255) aagcgcttcgcATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAc < 1:724464/143‑1 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAc > 2:1690216/1‑151 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAcc < 1:1409167/150‑1 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAc < 2:1250385/149‑1 (MQ=255) cGTAGCGCATCAGGCCATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgc > 2:1053502/1‑151 (MQ=255) gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt > 1:268393/1‑151 (MQ=255) aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagtagt > 2:1482015/1‑151 (MQ=255) ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg > 1:502851/1‑151 (MQ=255) cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg > 1:342861/1‑151 (MQ=255) cggaagtccTGTGCCGGGGGCTCCCTTGATTCTTTCCCGCAAACACCCGCGCTTTTACCCACCGCCTTGTGAAAAAA‑GCCCTC‑GGGGGCGGCTGTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg‑gg < 2:1005859/144‑1 (MQ=255) ccGCAGCCTTGTGCCGGCGCCTCCCGTGTTTCTTTCCCCCCAACACCCGCCCTTTTACCCACCGCAATTTGAAAAAT‑TCCTTT‑CGCGGCGGCTGTCTCCCATTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg < 2:340143/148‑1 (MQ=255) cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg‑gg > 2:234779/1‑148 (MQ=255) aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg‑ggg > 1:828424/1‑147 (MQ=255) cctggtgtCGGCGCCTCCCGTGTTTCATTCCCCCAAACACCCGCGCTTTTCCCCACCCCCATGTCAAAAAA‑TCCCTC‑CGCGGCGGCTTTCTGCCAGCTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCaa < 2:58109/144‑1 (MQ=255) tgtgTCAGCGGCTCCGGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCAACCCCTTGTAAAAAAA‑TCCCTT‑CGCGGCGGCTGTCTGCAAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAgaga < 2:858750/149‑1 (MQ=255) tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg‑gggggggg > 1:48924/1‑142 (MQ=255) gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAGGAGaaa > 2:1353738/1‑150 (MQ=255) tgtCGGCGGCTCCGTTAATTCTTTCCCCCAAACACCCGCGTATTTATCAACCGCAATGTGAAAAAA‑TGCCTG‑GGCGGCGGTTTTCTGCAAGTTTTTTTT‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAagag < 2:1546326/146‑1 (MQ=255) tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg‑gggcagaaaa > 2:1329351/1‑140 (MQ=255) cgggggCTCCCGTTTTTCTTTCCCCCAAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAA‑TCCCTT‑CGCGGCGGCTGTCTGCCAGTTATTTTC‑CCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAag < 2:636280/137‑1 (MQ=255) cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCACgagagaaa > 2:1504475/1‑144 (MQ=255) ggctccggtttttCATTCCCCCCAACACCCGCGTATTTACCCACCGCCATGTAAAAAAA‑TGCCTT‑CGCGGGGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTaa < 2:467460/140‑1 (MQ=255) cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGAGAGaaaaaaatttt > 1:1000947/1‑140 (MQ=255) ccGTTATTCTTTCCCCCCACCACCCGCGCATTTATCCACCGCC‑TGTTAAAAAACTCCCTC‑CGCGGCGGCTGTCTCCCGGTTATTTTC‑CCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGccc < 2:1408215/150‑1 (MQ=255) gTGATTCATTTCCGCCAAAAACCGCGCATTTATCCAACCCCATGTGCAAAAA‑TCCCTA‑CGCGGCGGCTGTCTGCCAGCTATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTc < 2:1176847/151‑1 (MQ=255) ttccccccAACACCCGCCCTTTTACCCACCGCCATGTAAAAAAT‑TGCCTC‑CGGGGCGGTTGTCTGCCAGTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCAtgt < 2:98028/141‑1 (MQ=255) cccccAACACCCGCATTTTTACCCACCGCTATGGGAAAAAT‑TGCCTG‑CGGGGCGGCTGTCTGCCATTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg < 2:751113/148‑1 (MQ=255) caacaaCCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGAGAGAAAAAACttttaaaccccc > 2:994518/1‑127 (MQ=25) ccccgtttttACCCACCCCTATGTGAAAAAA‑TCCTTC‑CGCGGCGGCTTTCCGCCAGTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt < 2:935348/138‑1 (MQ=255) ccGCGATTTTACCCACCGCCATGTGAAAAAA‑TGCCTC‑CGGGGGGGCTGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 2:1001024/150‑1 (MQ=255) cATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGAGGGGAAAAAATTTTAACCCCCCCTGTTGTGTggggatagccggtgg > 1:42786/1‑136 (MQ=21) tttcccaaccgcattgtgAAAAAA‑TGCCTC‑GGGGGCGGCTTTCTGCCAGCTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgc‑tg < 1:378935/132‑1 (MQ=255) tttACCCAACCCATTGTAAAAAAG‑TGCTTT‑GGGGGCGCTTGTCGGCCAGTTATTTTC‑CCGCCGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑aa < 1:220809/150‑1 (MQ=255) cgaggagaaaaa‑aGCCTTCCGGGGCGGCTGTCGGCCAGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑AACAATGTGGAAAg < 2:895920/144‑1 (MQ=255) cAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGAGAGAAAAAACTTTTAACACCCCCTGTTTTGTGGGAAATGCTGC‑TGGG‑Aaaaa > 2:1533516/1‑135 (MQ=255) aaaaaa‑TGCCCC‑CGGGGCGGCTGTCGCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑AACAATGTGGAAAGGGa < 1:1401790/146‑1 (MQ=255) aaa‑TGCCTC‑GGGGGCGGCT‑TTTGCCCGGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑AACAATGTGGAAAGGGAAGTgcgg < 2:1532354/151‑1 (MQ=255) a‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCGACGGAAAAAAACTTAATGCCGCCCTTTTCTGTGGGGAACGCTGC‑GGGTCGAAAATGTGGAAAg > 2:252001/1‑140 (MQ=255) tt‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGAAGAGAAAAAACCTTAAAACCGCCCTGTTTTGTGGGGATCGCTGC‑TGGT‑GACAAAGTGGAAAGGGGAgggggggggaac > 2:20743/1‑139 (MQ=255) ggcggcTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGAAGAGAAAAAAATTTTAAACCGCCCTGTTTTGTGGGCATTGCTGC‑TGGG‑GAAAAAGTTGGAAGGGGAGGGGGGGGGAACAGCTc > 2:1398953/1‑150 (MQ=255) cggcTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa < 2:214765/150‑1 (MQ=255) cggcTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGAGAGGAAAAAATTTTAACCCCTCCTGTTGTGTGGGGATCACCGC‑TGGG‑GAAAAAGAGGGAAGGGGGGGGGGGGGGATCAGCTcc > 1:634014/1‑149 (MQ=255) tctgcagttttttttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa < 1:120386/133‑1 (MQ=255) cTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGGAGGGAAAAAAGATTTAACCCCGCCAGTTTTTGGGGC‑GTTGCTGCTTGGG‑AAAAATGTGGAAGGGGAGGGGGGGGGAA‑CGACTCCACAACCgcgcacaat > 1:1214422/1‑145 (MQ=255) ttttcttccgcgcaaacaaATCCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt < 2:261810/134‑1 (MQ=255) ttttt‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGGGAGAGAAAAAAATTTTAGACCGCCCCTTTGTGTGGGGATCGCCGC‑TGGG‑GACAATGTGGGAAGGGAAGGGGGGGGGATCCACTCCCGaaaccccccaatt > 1:601692/1‑136 (MQ=255) ttttc‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa < 2:1639110/149‑1 (MQ=255) atttcccgccgataacaaAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGTCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg < 1:214386/136‑1 (MQ=255) ccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 2:783486/145‑1 (MQ=255) ccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCa < 2:867853/151‑1 (MQ=255) | TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAA‑TGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGG‑GGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGT‑AACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA > NC_000913/357741‑358034 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |