Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,532,988 | 0 | G | A | 43.5% | 19.2 / 11.1 | 23 | G194S (GGC→AGC) | cysK | cysteine synthase A, O‑acetylserine sulfhydrolase A subunit |
Reads supporting (aligned to +/- strand): ref base G (6/7); new base A (10/0); total (16/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.49e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGC‑ACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCG‑CTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAA > NC_000913/2532848‑2533133 | cGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACCCCGGCTGGCGGCAGCCGCCACACTAAAGGC‑ACCcaagacaaga > 1:342325/1‑147 (MQ=255) aaaCCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGC‑ACCAAAGGCAAGACCGATCtt < 1:124107/151‑1 (MQ=255) ccTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCCGACTGGCGTCAGCCGCCACATTAAAAGG‑CCCAAAAGGAAGAACGGTCtt > 1:175483/1‑147 (MQ=255) ccTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGCTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTAAACCTGCTGGTGTCAGGCGCCACATTAAAGGC‑AACAAAAGCGAAAACGGTCtttt > 1:577284/1‑147 (MQ=255) gAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGGAAGCTGGCTGGGGTGCACCGCGACATTTAATGC‑GACAAAAGCAAAAACGAACTTAttt > 2:475387/1‑146 (MQ=255) ttCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGCACGCTGACTGGGGTCACCCGCCCCATTTAAGGA‑ACCCAAAGGAAGaaccggctt > 1:325795/1‑132 (MQ=255) cGAAAAGACCACCGGTCCGGAGATTTGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACACAGACTTGCGTATGCCCCTTCAATAAAAAC‑ACCAAAAAGCAGACCGATCTTTTATCTCTCGCCGtt > 2:127729/1‑151 (MQ=255) cGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCCGACCGGCGGGCGCCCCCACACTAAAAGCAACAAAAGGAAAAACCG‑TCTTATCTCTGTCCCCgct > 1:139179/1‑149 (MQ=255) aGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGAACCCCGGCTGGGGTTAGCCGCCACATTTAAGGG‑ACCAAAAGCAAAACCGATCTTATTTTTCTCGCCGTTgtgac > 2:275477/1‑147 (MQ=255) gTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGC‑ACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATtc < 2:321445/151‑1 (MQ=255) agaTGTGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGC‑ACCAAAGGCAAATCCGAGCTTATCTCTGTCGCCGTTGGGCCAACCGATTCTcc < 2:606183/148‑1 (MQ=255) agaTATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGCACGCTGATTTGCGTCAGCCGCTACCCTAAAAGC‑AACACAAGCAAGACCCAACTTTTTCCTTTCGCCGCTGAGCCAACCCATtctc > 1:482914/1‑147 (MQ=255) gATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGGACTCTGACTGGCGTTAACCGCTACATTAAAAGG‑AACAAAAGCAAAAACGATCTTATTTCTGTCGCCGCTGAGCCAACCGGGTCTCCCGTTATCGccc > 2:245701/1‑148 (MQ=255) tACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCCGACCGGGGGCAGCCGCCCCATTAAAAGC‑ACCAACAGCAAGAACCAGCTTTTCTCTTTCGCCCGTGAGCCAACCGCTTCTCCAGCGATCGccc > 2:441098/1‑146 (MQ=255) ccGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGC‑ACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCG‑Ct < 2:435003/151‑1 (MQ=255) tCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTAAGCTGAGAGGCGTCCGCCGCTACATTAAAAGC‑ACCAAAAGAAAAAACGATCTTTTTTCTGTCGCCGTTGAGGCAACCCATTCTCCACGTATCGCCCCGGCG‑CTGGCCggg > 2:182079/1‑150 (MQ=255) gTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGCAACCAGACTGGCGGGAGCCGCCACAATAAATGG‑ACCAAAGGCAAGACCCAAATTATTTTTGTTTCCCTTTAGACAAAAGATTTTCCCGTTATTGCCCCGGCGCCTGGCAg > 1:240403/1‑146 (MQ=255) gATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGC‑ACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCG‑CTGGCAGGTGAagag < 2:83166/150‑1 (MQ=255) tGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGC‑ACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCG‑CTGGCAGGTGAAGAGATTAAACCTgg < 1:76997/150‑1 (MQ=255) cGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGC‑ACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCG‑CTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGt > 1:35006/1‑151 (MQ=255) tCAGACGCTACATTAAAGGC‑ACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCG‑CTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTAt > 2:538322/1‑151 (MQ=255) aGCCGCTACATTAAAGGC‑ACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCG‑CTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATcc > 1:328562/1‑151 (MQ=255) gCTACATTAAAGGC‑ACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCG‑CTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGc > 1:331593/1‑151 (MQ=255) ccAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCG‑CTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCaa < 2:252717/151‑1 (MQ=255) | CGGCAAACCCTGAAATTCACGAAAAGACCACCGGTCCGGAGATATGGGAAGATACCGACGGTCAGGTTGATGTATTTATTGCTGGCGTTGGGACTGGCGGTACGCTGACTGGCGTCAGCCGCTACATTAAAGGC‑ACCAAAGGCAAGACCGATCTTATCTCTGTCGCCGTTGAGCCAACCGATTCTCCAGTTATCGCCCAGGCG‑CTGGCAGGTGAAGAGATTAAACCTGGCCCGCATAAAATTCAGGGTATTGGCGCTGGTTTTATCCCGGCTAACCTCGATCTCAA > NC_000913/2532848‑2533133 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |