Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 62.7% | ‑4.6 / 45.6 | 51 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (9/10); new base A (0/32); total (9/42) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.04e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-03 |
AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG > NC_000913/357745‑358037 | aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑c < 1:527466/151‑1 (MQ=255) tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTtctctc > 2:1002227/1‑151 (MQ=255) gcaaaTCCCAGCCGAGGCCTTGTGTCAGCGGCTCCCGTGA‑TTCATTCCCGCAAA‑CACCCGCGCTTTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTC‑CGCGGCGGCTGTCTGCAAGTTATTTTTCCGCCGCAAAAAAA‑ccgccg < 1:1185898/138‑1 (MQ=17) tcccACCCGCAGCCTGGTGCCGGGGGCTACGGTGT‑TTATTTCCCGCAAA‑AACCCGCGCTTTTA‑CCCAACCCC‑TTGTAAAAAAATGCTTC‑CGGGGGGGCTGTCGGCCAGTTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCagtagt < 2:1211666/148‑1 (MQ=255) tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTgg > 1:851415/1‑151 (MQ=255) ccggCCGGAGCCTTGTGCCGGCGGCTACCTTGA‑TTTTTTCCCCCAAA‑CACCCGCACATTTA‑CCCACCCCC‑TTTGGAAAAATTGCCTC‑CGCGGCGGCTGTCTGCAAGTTTTTTTCCCGCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTg < 2:1015942/146‑1 (MQ=255) cAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑g > 2:998059/1‑151 (MQ=255) aGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gg > 2:97630/1‑151 (MQ=255) ccGCAGCCCTGGGCCAGCGCCTACCGTAT‑TTCTTTCCCCCCAA‑CACCCGCGCTTTTA‑CCCACCGCC‑ATGTGAAAAAATCCCTC‑CGGGGCGGTTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑ggg < 2:642170/150‑1 (MQ=255) ccGCAGACCTGT‑TCAGCCGGAACCCTTATTTCATTCCCGCAAA‑CACCCGCG‑ATTTACCCCAACGCA‑AGGTAC‑AAAATGCCTTGCGCGGGGGCTGTCTGCAACCTATTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gggg < 2:1198551/151‑1 (MQ=255) cgcaacctggtgcCAGCGGCTCCCGTTT‑TTCATTCCCCCCAA‑CACCCGCGCTTTAA‑CCCAACGCC‑TTGTAAAAAAATCCCTC‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gggg < 1:1209561/141‑1 (MQ=255) cgaagcccTGTGTCAGGGGCTCCCGTGT‑TTCTTTCCCGCCAC‑CACCCGCGCTTTTA‑‑CCAACCCCAATGTGAAAAATTGCCTC‑CGGGGCGGCTTTCTGCTAGTTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGc < 2:1139967/145‑1 (MQ=255) cggccgtttgccggCGTCTACCGTAA‑TTCTTTCCCCCCAC‑CACCCGCGTTTTTA‑CCCACCCCC‑ATGAGAAAAAAGGCCTC‑CGGGGGGGTTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGc < 1:922202/136‑1 (MQ=255) gACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑gg > 2:187327/1‑144 (MQ=255) cggcgcCTACCGTAA‑TTCTTTCCCCCCAA‑CACCCCCCCTTTTA‑TCCAACCCA‑ATGTGCAAAAATGCCTT‑CGGGGCGGCTGTCGGCCCGTTATTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAg < 2:515583/147‑1 (MQ=255) gcctccgttGA‑TTCTTTCCCCCCAA‑CACCCCCGCATTTA‑CCCACCGCC‑ATGTGAAAAAATCCCTC‑CGCGGCGGTTTCCGGCCAGTTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATa < 1:115389/136‑1 (MQ=255) gcctcccGTGA‑TTCTTTCCCCCCAA‑CACCCGCGCTTTTA‑CCCACCGCC‑ATGTGAAAAAATGCCTC‑CGGGGGGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCt < 2:825477/142‑1 (MQ=255) tACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAGAttt > 1:584735/1‑146 (MQ=255) cccGTGT‑TTCTTTCCCCCCAA‑CACCCGCCCTTTTA‑‑CCAACCCCAAGGTGAAAAAATCCCTC‑CGGGGGGGCTGTCGGCAAGTTATTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGccc < 1:139969/150‑1 (MQ=255) cccGTGT‑TTC‑TTTCCCCCAACCACCCGCGCTTTTA‑CCCACCGCC‑TGGTAAAAAAATGCCTC‑CGGGGCGGTTGTCGGCAAGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGc < 1:1061724/148‑1 (MQ=255) aCCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAATTTaaaacc > 2:556077/1‑147 (MQ=25) ccGTGT‑TTTTTTCCCCCCAA‑CACCCGCGCATTTA‑CCCAACGCC‑ATGTGCAAAAATGCCTC‑CGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGccc < 2:159503/150‑1 (MQ=25) ccGTGA‑TTCTTTCCCCCCAA‑CACCCGCGCATTTA‑CCCACCGCC‑TTGTGAAAAAATGCCTT‑CGCGGCGGCTGTCGGCCAGTTTTTTTTCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGc < 2:788332/148‑1 (MQ=255) tGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGAGgaaaaagttt > 2:592861/1‑134 (MQ=25) gA‑TTCTTTCCCCCCAA‑CACCCGCCTATTTA‑CCCACCCCC‑AGGTCCAAAAAGCCCTC‑CGGGGGGCTTGTCCGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGt < 1:171552/148‑1 (MQ=255) ttCTTTCCCGCAAA‑CACCCGCGATTTTA‑CCAAACGCC‑ATGTGAAAAATTGCCTC‑CGGGGGGGTTGTCGGCCACTTATTTTCCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg < 2:452240/150‑1 (MQ=255) ccccccAA‑CACCCCCCCTTTTA‑CCCACCGCC‑ATGTGAAAAAATCCCTC‑CGCGGCGGCTGTCTGCAAGCTTTTTTCCCGCCAAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg < 1:556026/147‑1 (MQ=255) acccGCGCATTTA‑CCCACCGCA‑ATGTGAAAAAATCCCTT‑CGGGGCGGTTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc < 2:845587/146‑1 (MQ=255) tttA‑TCCACCGCC‑ATGAAAAAAAGTCCCTG‑CGCGGGGGTTTCCGCCCAGTTTTTTTCCCCCCACAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt < 1:786247/148‑1 (MQ=255) ccACCGCC‑ATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt < 1:370226/149‑1 (MQ=255) ccAACGCC‑ATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa < 2:954857/148‑1 (MQ=255) cGCC‑ATGTAAAAAAATGCCTT‑CGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGACAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt < 2:14748/147‑1 (MQ=255) cc‑ATGTGAAAAAAGCCCTC‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa < 2:449158/150‑1 (MQ=255) aaaaaTGCTTG‑CGTGGCGGATGACTGCAACAAATTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg < 2:462029/148‑1 (MQ=255) aaaaaGGCCTG‑CGCGGCGGCTTTCTCCCATCTTTTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt < 1:568480/149‑1 (MQ=255) aaaccctcccgcgCCGGTTGTCAGCCGTTTTTTTTCCCCCCCCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg < 1:278901/142‑1 (MQ=255) gggggggcTGTCGGCCGGTTTTTTTCCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc < 2:1066980/142‑1 (MQ=255) cggcTTTCCGCCAGCTATTTCCCCGCCACAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa < 1:812295/150‑1 (MQ=255) cTGCCAGTTTTTTTCCCCCCCCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg < 1:1143432/144‑1 (MQ=255) ggccggttATTTTCCCGCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 1:1218756/139‑1 (MQ=255) ccAGTTTTTTTCCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 1:145731/145‑1 (MQ=255) ccactttttttCCCACCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa < 1:1002127/145‑1 (MQ=255) tatttttcCCGCCGAAAAAAAC‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa < 1:1230355/127‑1 (MQ=255) ttttttcCCGCCAAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa < 1:1064137/145‑1 (MQ=255) ttttttcCCCCCGAAAAAAAT‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt < 1:899604/146‑1 (MQ=255) tttttccccgccgCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg < 1:604086/143‑1 (MQ=255) tttttcCCCCCGCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg < 2:551856/148‑1 (MQ=255) ttttcCCGCCGAAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa < 2:25172/149‑1 (MQ=255) ttttcCCGCCGAAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg < 2:694655/150‑1 (MQ=255) ttcccgccgCAAAAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt < 2:124324/147‑1 (MQ=255) ccgcAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAAATTTAACCCGCCCTGTTTTGTGGGCATTGCTGCTGGTAACAATATGGGAAGGGGAGAGGGGGGGAACAGACCCAGGACCCCGCCCAATGAGttt > 2:73755/1‑144 (MQ=255) caacaaAA‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAGGAAAAAAAATTTAAG‑CCGTCCTGTTGTGTGGGGAACGCTGCTGGGGACAAAATGGGAAAGGGAGAGGGGGGGATCAGCCCCCCAACCCCGCCAATTGAGTTTGCCCCCATg > 2:984988/1‑151 (MQ=255) | AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑CCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATG > NC_000913/357745‑358037 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |