| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_vch_pgi | 543,091 | 0 | G | C | 23.8% | 30.6 / 11.5 | 23 | intergenic (+58/+170) | glxK/allE | glycerate kinase II/S‑ureidoglycine aminohydrolase |
| Reads supporting (aligned to +/- strand): ref base G (6/10); new base C (3/2); total (9/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.11e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGC‑CATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATA > NC_000913_3_vch_pgi/542950‑543231 | cTCGCCAGGGGTGAAACCAATCTCTTCAACAGCGCGCGAAATACTGCCTGCGC‑CATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGc < 2:52434/149‑1 (MQ=255) ccAGCGGTGAACGCAATCTCTTCAACAGCGCGGGAAATATTGCCTGCGC‑CATTAAAATAGGTCAGGGAATTAAAAACTATCCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGa < 1:318025/149‑1 (MQ=255) gCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGC‑CATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa < 1:402578/146‑1 (MQ=255) gCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGC‑CATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGa > 2:402578/1‑146 (MQ=255) ggTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGC‑CATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCTGATGCGAGTTATATCGCATAAGAGTGCAGTACTCATCCCGGATGCGGCATGAGTAcc > 2:101201/1‑149 (MQ=255) gTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGC‑CATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAATGCAGTACTCATGCCGGATGCGACATGAGTACCa > 2:363573/1‑149 (MQ=255) aCCAATCTCTTCAACAGCGCGCGAAATATT‑CCAACTCTTTTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCAGAtcc > 2:397617/1‑149 (MQ=255) aTCTCTTCAACAGCGCGCGAAATATTGCCTGCGC‑CATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCGGTACTCATGCCGCATCCGGCATGAGTActgcact < 1:108755/144‑7 (MQ=39) aTCTCTTCAACAGCGCGCGAAATATTGCCTGCGC‑CATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCGGTACTCATGCCGCATCCGGCATGAGTActgcact > 2:108755/1‑138 (MQ=39) cAGCGCGCGAAATATTGCCTGCGC‑CATTAAAATAGGTCAGGGGATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATCAGTACCATATCCTTCCTGAAAATcgc < 1:455282/149‑1 (MQ=255) aGCGCGCGAAATATTGCCTGCGC‑CATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCCTCCTGAAAATcgcg > 1:297902/1‑149 (MQ=255) cTGCGC‑CATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGc < 1:195454/149‑1 (MQ=255) cgtcATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGATTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCaga < 2:241781/146‑1 (MQ=255) aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttc > 1:108843/1‑113 (MQ=255) aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtccttc < 2:108843/113‑1 (MQ=255) aaaTAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTGCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGt > 1:431671/1‑149 (MQ=255) ggTCAGGGAAGTAAAAACTAACGCTTAGCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAACTCGCGCAACTTCTATATATTGCAGAGATCATGTAGGcc < 2:540868/149‑1 (MQ=255) aaaaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCCGAAAATCGCAAAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcgaag < 1:131529/149‑1 (MQ=255) tggaCCTTTAAAGCGGAGGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGAGCATGTAGGCCTGATAAGCGAAGCGCATCAGGCa < 2:543499/146‑1 (MQ=255) ttACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGACGATGTATGCCTGATAAGCGAAGCGCATCAGGCaa < 2:536185/149‑1 (MQ=255) ttAAAGCGGATGGCATTTATATCGCATAAGAGTGCAGTATTCATGCCGTATTCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTAc < 1:46577/149‑1 (MQ=38) aaGAGTGGAGTGGTAATGGCGTATGCGGGATGAGTAACATATCCCTCCTGAAAATCGCGCAAATTCTCTTTAGTGTAGAGATTATGTAGGCCGGATAAGCGTAGCGCATCAGGCAATGTTATAAAAAAAGCCACGGTATAAACCGTGGc < 2:15195/149‑1 (MQ=25) tCATGCCGCATCCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCCACATa > 1:536278/1‑149 (MQ=255) | CTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGC‑CATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTATAAACCGTGGCAAAATCCAACATA > NC_000913_3_vch_pgi/542950‑543231 |
| 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: ‑ |