| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,536,370 | 0 | C | G | 30.0% | 21.3 / 15.0 | 20 | 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 (7/7); new base G (3/3); total (10/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCA > NC_000913/2536238‑2536507 | ccTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAtttttt > 6:154633/1‑149 (MQ=255) tCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGTTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGAGCCATTTTTTTATGCt > 4:106749/1‑149 (MQ=255) tGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAAt < 5:154633/149‑1 (MQ=255) gTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTAc < 6:100526/149‑1 (MQ=255) gTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCg > 1:95309/1‑149 (MQ=255) gTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCg > 6:87533/1‑149 (MQ=255) tAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCgg < 5:87533/149‑1 (MQ=255) ttCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCa > 2:59481/1‑113 (MQ=255) ttCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCa < 1:59481/113‑1 (MQ=255) gCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGc > 2:74648/1‑148 (MQ=255) gcataaaaaaaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTc > 5:41519/7‑149 (MQ=255) ataaaaaaaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTccc < 6:41519/145‑1 (MQ=255) taaaaaaaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCa < 3:27710/99‑1 (MQ=255) taaaaaaaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCa > 4:27710/4‑102 (MQ=255) gTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAATTCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCgg < 5:132992/135‑1 (MQ=255) gTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAATTCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCgg > 6:132992/1‑135 (MQ=255) aaaaaaaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAAc < 4:99267/137‑1 (MQ=255) aaaaaaaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAAc > 3:99267/3‑139 (MQ=255) aaaaTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCCGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCa < 8:126986/149‑1 (MQ=255) aaaaTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTgg < 2:112923/79‑1 (MQ=255) | CCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCA > NC_000913/2536238‑2536507 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |