Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,510,395 | 0 | C | A | 44.4% | 33.4 / 25.9 | 36 | S131R (AGC→AGA) | ydcR | putative DNA‑binding transcriptional regulator and putative aminotransferase |
Reads supporting (aligned to +/- strand): ref base C (19/1); new base A (0/16); total (19/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.33e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01 |
CGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTT > NC_000913/1510249‑1510532 | cGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCaaaa < 1:340999/151‑1 (MQ=255) aTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGc > 1:163368/1‑151 (MQ=255) gATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGAt > 2:624537/1‑149 (MQ=255) atatTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGATCGTTCCGCTTGTCTCGGCCTTT‑CCCGACCAGCGACTTT‑TTGCCTTCTAACTACTAAA‑CCGCTCGCGGGCGCAGGACCGCAAAACCGCCACAGCGATGATCGTGAt > 2:408216/1‑145 (MQ=255) gATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGc > 1:587481/1‑151 (MQ=255) atatGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGcc > 1:404812/1‑151 (MQ=255) atatGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGcc > 2:372871/1‑151 (MQ=255) tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg > 1:490923/1‑151 (MQ=255) cAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACg > 1:16514/1‑151 (MQ=255) gcccAGCCGCGTCCCGTGGGTTGTCCCGTTTCCCCGGGCTTTC‑CCCGCCCCGCGTTTTT‑CCCCCCCCAAAAAACAAAC‑CCGCCCGGGGGCGAAGTTAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc < 1:91398/148‑1 (MQ=255) gcccAGCCGCGCCCCGTCGGTC‑TTCCTTTTCCCCCCGCCTTTCCCCCCCCCGCGTTTTT‑CCCCCCCCAAAAAACAAAC‑CCCCCCGCGGGCGCAGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc < 1:662267/148‑1 (MQ=255) ccggccgcgccccGTGGGTCGTCCCTTTTCCCT‑GGCCTTTCCCCGCCCGGCGCTTTT‑CCCCCCCCCAAAAACAAAC‑CCCCCCGCTGGCGCAGGTAACAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAg < 1:31341/139‑1 (MQ=255) ccagccccaccccgccgGTCGTCCCTTTGCCCCCGCCCTTC‑CCCCCCCCGCGTTTTT‑CCCCCCCCAAACAACAAAC‑CCGCTCGGGGCCGAAGATAAGAAAATCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACg < 2:566616/136‑1 (MQ=255) gtcggtcgTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCt > 2:415918/1‑151 (MQ=255) tcggtcgTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGa > 1:607887/1‑113 (MQ=255) cggtcgTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAt > 2:47464/1‑151 (MQ=255) ttttGCCCCGGCCTTTCCCCGCCCCGC‑ACTTTCCCCCCCCCCAAAAAAAAAC‑CCGCTCGGGGGCGAGGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAt < 2:628437/141‑1 (MQ=255) ttggCCGCGCCCTTC‑CCCCCCCCGC‑ATTTT‑TCCCCCCCCAAAAAAAAAACCCTCCCCTTGGCGGAGAAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCtcgt < 1:517623/144‑1 (MQ=255) ttccccccgCCTTTCCCCGCCCCGCTTTTTT‑CCCCCCCCCAAAAACAAAC‑CCGCTCGCGGGCGAGGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGc < 1:706297/129‑1 (MQ=255) ggCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAC‑CCGCTCGCTGGCGCAGGTAAGCAACACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCtt > 2:285066/1‑150 (MQ=255) gCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTAc > 2:233736/1‑151 (MQ=255) cTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTAc > 1:298582/1‑149 (MQ=255) ccccccgctttttccccccccccAACACCAAAC‑CCGCCCGCGGGCGCAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCAt < 1:469412/135‑1 (MQ=255) cgcactttcccccccTCCAACACAAACC‑CCGCGCGCGGCGGCAGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCAcc < 2:423392/139‑1 (MQ=255) cTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtctc > 1:244062/1‑151 (MQ=255) ttccccccccccAACAAAAAAC‑CCCCCCGCGGCCGCAGAAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtct < 2:352438/145‑1 (MQ=255) ttccccccccCAAACAACAAAC‑CCGCCCGCGGGCGCAGGAAACAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtc < 1:488018/144‑1 (MQ=255) ccccccCAAACAACAAAC‑CCGCCCGTGGGCGAAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgat < 1:520320/151‑1 (MQ=255) ccccccCAAACAACAAAC‑CCGCCCGGGGGCGCAGAAAACAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgat < 2:201185/151‑1 (MQ=255) ccccccCAAACAAAAAAC‑CCCCCCGCGGGCGAAGAAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAttt < 2:536805/143‑1 (MQ=255) ccTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgatga > 2:344070/1‑150 (MQ=255) caacaaCTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTg > 2:379102/1‑151 (MQ=255) aacTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCAt > 2:12061/1‑151 (MQ=255) ggtcaaagtaaaGAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTGCCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCg < 1:549383/134‑1 (MQ=255) ggCGCAGGAAAGAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCg < 1:75446/148‑1 (MQ=255) cgcAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGtt > 1:188617/1‑148 (MQ=255) | CGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGACCCGCGACTTT‑TCCCCCTCCAACAACTAAA‑CCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTT > NC_000913/1510249‑1510532 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |