Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_vch_pgi | 860,157 | 0 | C | T | 40.0% | 1.1 / 13.1 | 15 | intergenic (‑129/+17) | ybiV/ybiW | sugar phosphatase; fructose‑1‑P/ribose‑5‑P/glucose‑6‑P phosphatase/putative pyruvate formate lyase |
Reads supporting (aligned to +/- strand): ref base C (7/2); new base T (2/4); total (9/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.36e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCAT > NC_000913_3_vch_pgi/860013‑860301 | tAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACgccgcc > 2:320001/1‑149 (MQ=255) tAACTTTTACGCTCATAGGAGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGAcgacgct < 1:372000/149‑2 (MQ=255) aCTTTTACGCTCATAGGTGTCAGCCTGGCGGGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCgg < 1:8006/149‑1 (MQ=255) aCTTTTACGCTCATAAGTGTTTTCCTAGCTCTAAAATGAAACCTTAGCAGTTTAATAAAAATTTTACTGAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATTTCAATGATTTAGCCCGGCGACGACGCCGCCgg < 1:127713/149‑1 (MQ=255) ttttACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGc > 1:126399/1‑149 (MQ=255) gCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaa > 1:161393/1‑143 (MQ=255) gCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaa < 2:161393/147‑5 (MQ=255) cTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaatca < 1:319321/149‑8 (MQ=255) cTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCCGAGCACAAATTATATTTCGTATGAAAGTAAGGATGAAATTCATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGCCCGAGGAGCTTACAGCAGATGtt > 2:204119/1‑149 (MQ=255) ggCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCa < 1:320001/149‑1 (MQ=255) taaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGAAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGtt > 2:263096/1‑148 (MQ=255) aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCgg < 1:368190/59‑1 (MQ=37) aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCgg > 2:368190/1‑59 (MQ=37) aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGAGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCcg > 2:331115/1‑149 (MQ=255) cGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGGGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGGCGGCGACACGCACTACCAGATCGCGATACTGAGCGGCGTGtttg > 1:68584/1‑148 (MQ=255) cGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGACCGGGATGtttt > 1:107959/1‑149 (MQ=255) gacgCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCAt > 1:351610/1‑149 (MQ=255) | TAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCAT > NC_000913_3_vch_pgi/860013‑860301 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |