| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,734,269 | 0 | T | G | 51.8% | 26.8 / 37.5 | 27 | V402G (GTG→GGG) | avtA | valine‑‑pyruvate transaminase |
| Reads supporting (aligned to +/- strand): ref base T (1/12); new base G (14/0); total (15/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.48e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATT > NZ_CP009273/3734187‑3734358 | gggCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGGAGAtt > 2:165669/1‑90 (MQ=255) cTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTg > 1:131679/1‑90 (MQ=255) aaCCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAag < 1:77516/90‑1 (MQ=255) tGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATa < 1:62836/90‑1 (MQ=255) cGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGaaa < 2:304513/90‑1 (MQ=255) gCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGGAGATTCTGGCGGAAGAGATAGAAag > 2:393409/1‑90 (MQ=255) aTACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGGAGAGATAGAAaaag > 2:221522/1‑90 (MQ=255) cGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCt < 1:353158/90‑1 (MQ=255) cGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCt < 1:355291/90‑1 (MQ=255) aTGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGa > 1:88871/1‑90 (MQ=255) tATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAg > 1:93911/1‑90 (MQ=255) tGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGAt < 1:203701/67‑1 (MQ=255) tGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGAt > 2:203701/1‑67 (MQ=255) aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTaa > 1:67802/1‑90 (MQ=255) aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGGCACTaa > 1:45487/1‑90 (MQ=255) aCGTACCAGAGCCGGAGAAAATTGAGGCGTGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGc > 2:374379/1‑90 (MQ=255) aCGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAAAGCCTGGGCTGGAAGTCACTAAAAATTGc > 2:131898/1‑90 (MQ=255) gTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGGAAAGAAAGAAAAAACCTGGGGTGGAAGGCACTAAAAATTGCCt > 2:4729/1‑90 (MQ=255) gTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAAAGGCTGGGCTGGAAGGCACTAAAAAATGCCt > 2:234746/1‑90 (MQ=255) aCCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGa < 1:120425/90‑1 (MQ=255) agCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgcgc < 2:58269/90‑1 (MQ=255) agCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgcgc < 1:358909/90‑1 (MQ=255) agCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgcgc < 1:162172/90‑1 (MQ=255) gCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGaaa < 1:388978/50‑1 (MQ=255) gCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGaaa > 2:388978/1‑50 (MQ=255) gCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAAATGCCTGATGCGCt > 2:22270/1‑90 (MQ=255) aaTTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATc < 1:221898/90‑1 (MQ=255) tGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGAtt < 1:36186/90‑1 (MQ=255) | GGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATT > NZ_CP009273/3734187‑3734358 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |