| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_bme_pgi | 1,292,474 | 0 | T | . | 30.8% | 20.7 / 10.6 | 15 | intergenic (+109/+35) | galU/hns | glucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS |
| Reads supporting (aligned to +/- strand): ref base T (4/5); new base . (3/1); total (8/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.59e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCC‑GCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGT > NC_000913_3_bme_pgi/1292370‑1292621 | cGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCC‑GCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGa > 3:41404/1‑150 (MQ=255) accACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAAATGTCTTAAACCGGACAATAAAAAATACCGCC‑GCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgtc < 1:143484/149‑1 (MQ=255) aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCC‑GCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAg > 3:64216/1‑150 (MQ=255) tGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCC‑GCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTg < 4:64216/150‑1 (MQ=255) tttttttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCC‑GCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttt > 4:6557/1‑150 (MQ=255) tattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCC‑GCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCt < 1:15248/84‑1 (MQ=255) tattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCC‑GCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCt > 2:15248/1‑84 (MQ=255) aTTGTCTTAAACCGGACAATAAAAAATCCCGCC‑GCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCtt < 3:6557/150‑1 (MQ=255) taaaaTCCCGCCAGC‑GGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGttt > 1:12006/2‑150 (MQ=255) taaaaTCCCGCCAGC‑GGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGttt > 3:89856/2‑150 (MQ=255) taaaaTCCCGCC‑AGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTc > 1:125958/2‑91 (MQ=39) taaaaTCCCGCC‑AGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTc < 2:125958/90‑1 (MQ=39) aaaTCCCGCCAGC‑GGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAgg < 1:2932/63‑1 (MQ=37) aaaTCCCGCCAGC‑GGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAgg > 2:2932/1‑63 (MQ=37) gcTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGt < 1:128277/150‑1 (MQ=255) | CGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCC‑GCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGT > NC_000913_3_bme_pgi/1292370‑1292621 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |