| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_vch_pgi | 2,536,370 | 0 | C | G | 24.0% | 46.3 / 12.9 | 25 | intergenic (+27/+16) | crr/pdxK | glucose‑specific enzyme IIA component of PTS/pyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase |
| Reads supporting (aligned to +/- strand): ref base C (9/10); new base G (3/3); total (12/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACT > NC_000913_3_vch_pgi/2536228‑2536517 | ccAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGgc > 1:284182/1‑149 (MQ=255) aaGTCCACCCTGACTCCGTTTGTTATCTACAACATGGACGAAATCAAAGAACTAACGAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGgcgc < 1:317394/149‑1 (MQ=255) tCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAt < 2:220118/149‑1 (MQ=255) tCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGCCAtttttcactgcg < 1:248497/149‑8 (MQ=39) tcCAACATAGACGAAATCAAAGAACTGATAAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGATATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAACAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGACAtcggcg > 2:168059/1‑149 (MQ=255) aTGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCaaaa > 1:258010/1‑149 (MQ=255) aCGCAATCAAAGATCTGATCAAACTGTCCGGTAGTGTAACCGTGGGTGAACCCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCCGCGGCGGCAAAATCaa < 1:227644/149‑1 (MQ=255) tCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGCCAtt > 1:66713/1‑110 (MQ=255) tCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGCCAtt < 2:66713/110‑1 (MQ=255) aaGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCACTTCATCGct > 2:246745/1‑149 (MQ=255) aCTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTgg < 1:134704/149‑1 (MQ=255) gTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGctc > 1:137073/1‑129 (MQ=255) gTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGctc < 2:137073/129‑1 (MQ=255) gTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAgcg > 1:213170/1‑149 (MQ=255) aGCGTAACCGTGGGTGAATCCTCGGTTATCCGCAACAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCa < 1:83514/149‑1 (MQ=255) gggTAAGCGTGGGGGAATCCCCAGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCAt < 2:16332/147‑1 (MQ=255) tAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCCGGGTGTCGCCCATCAc > 2:21508/1‑149 (MQ=255) ccGTGGGTGAAACCCCGGTTGTCCGCATCAAGAAGTACTTCTTGCCGCAGTGAAAGATGGCGCCGATCGGCGCCATTTTTTTATGCTTACGCCGGGGGCGGCAAAGTCAATTCGTCGCTCTCCTGCTGCTGGGTGTTGGGCATCACtcg > 2:32500/1‑147 (MQ=255) ccGTCGGTGAAACCCCGGTTATCAGCATCAAGAAGTAATTCTTGCTGCAGTGAAAAATGGCGTCCATCGGCGCCATTTCTTTATGCTACCGACAACGGAGGCAATATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTc < 1:273605/149‑1 (MQ=255) ccccGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCCGTCCGCGCac > 2:276006/1‑148 (MQ=255) aTTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGtt > 1:217941/1‑149 (MQ=255) tGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCtt < 1:160650/149‑1 (MQ=255) tGAAAAAGGGCGCCCATCGGCGCCATTTTTTTATCCTTCAGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCCTCCCtt < 2:24645/149‑1 (MQ=255) aTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGgt > 2:132412/1‑78 (MQ=39) aTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGgt < 1:132412/78‑1 (MQ=39) ccATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACt > 1:196608/1‑149 (MQ=255) | CCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACT > NC_000913_3_vch_pgi/2536228‑2536517 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |