Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 585,386 | 0 | A | G | 45.2% | 29.3 / 12.6 | 31 | L342P (CTC→CCC) | nfrA | bacteriophage adsorption protein NfrA |
Reads supporting (aligned to +/- strand): ref base A (0/17); new base G (14/0); total (14/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.77e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTT > NZ_CP009273/585309‑585471 | ccAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCa < 1:1274053/90‑1 (MQ=255) agCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTgcgc < 1:416981/90‑1 (MQ=255) cAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGcgt < 1:83097/90‑1 (MQ=255) cAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGATGcgt < 1:252950/90‑1 (MQ=255) gTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGGGGGGGGGGGGCAGGTTTTTTGATGCGGCGTACTg > 2:252950/1‑90 (MQ=255) ttACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGa < 1:993839/90‑1 (MQ=255) aCGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGAcc < 1:1392582/90‑1 (MQ=255) aCGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGAcc < 1:1510000/90‑1 (MQ=255) cGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGAAGTTTTTGCGCTGCGTCGTACTGACct < 2:729009/90‑1 (MQ=255) cGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGGGGGGGGGGGGCAGGTTTTTCGCCGCGTCGTACTGAACTTCttt > 2:319277/1‑90 (MQ=255) gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGGGGTTTTTGGGC‑GCGTCGTACAGGCCTTCtttt > 2:1032246/1‑90 (MQ=255) gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGGGGGGGGGGGGGGGGTTTTTTGCTGCGGCGTGCTTTCCTTCtttt > 1:260723/1‑90 (MQ=255) gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGGGGGGGGGGAGCAGTTTTTTTGCTGCGTCGTTCTGACCTTCtttt > 2:1525496/1‑90 (MQ=255) cACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTaa < 1:853775/90‑1 (MQ=255) cACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTaa < 1:873476/90‑1 (MQ=255) tGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCa < 1:683206/90‑1 (MQ=255) cGCTCCTCAAGCATTTCATTGGC‑GGGGGGGGGGGGGGGGAGGTTTTTCGCTGCGTCGTACTTACCTTCTTTTTACAACACCCGGGGCGTc > 1:166699/1‑90 (MQ=255) tcctcAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTcgcg < 1:216922/90‑1 (MQ=255) cctcAAGCATTTCATTGGC‑GGGGGGGGGGGGGGGGAGGTTTTTTGCTGCGTCGGTCTTACCCTCTTTTTACAGCACCCGGAGCGgcgcgc > 1:420318/1‑90 (MQ=255) cctcAAGCATTTCATTGGC‑GGGGGGGGGGGGGGGGAGGTTTTTCGCTGCGGCGTACTGACCTTCTTTTTACAGCACCCGGAGCGTCgcgc > 2:765548/1‑90 (MQ=255) cctcAAGCATTTCATTGGC‑GGGGGGGGGGGGGGGCAGTTTTTTTGCTGCGTCGTTCTGACCTTCTTTTTAAAACACCCGGAgcggcggcc > 2:974029/1‑87 (MQ=255) cctcAAGCATTTCATTGGC‑GGGGGGGGGGGGGGGCAGGTTTTTCGCTGCGTCGTACTGACCTTCTTTTTAAAGCACCCGGAGGGgcgcgc > 2:1094935/1‑90 (MQ=255) ctcAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGcc < 1:1094935/90‑1 (MQ=255) ctcAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGcc < 2:1533915/90‑1 (MQ=255) aaGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCcaac < 2:1472541/90‑1 (MQ=255) aaGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCcaac < 1:423707/90‑1 (MQ=255) aTTTCATTGGC‑GGGGGGGGGGGGGGGCCGGTTTTTCGCTGCGTCGTACTGACCTTCTTTTTACAGCACCCGGAGGGGCGCGCCAacaaca > 1:1533915/1‑90 (MQ=255) tttCATTGGC‑GGGGGGGGGGGGGGGGAGGTTTTTCGCTGCGGCGTGCTGACCTTTTTTTTAAAGCACCGGGAGGGTCGCGCCAacaacaa > 1:997881/1‑89 (MQ=255) ttCATTGGC‑GGGGGGGGGGGGGGGCAGGTTTTTTGCTGCGTCGTACTGACCTTCTTTTTACAGCACCGGGAGCGTCGCGCCCACAACATa > 2:423707/1‑90 (MQ=255) ttCATTGGC‑GGGGGGGGGGGGGGGCAGGTTTTTCGCTGCGTCGTACTGACCTTCTTTTTAAAACACCGGGAGCGGCGCGCCCACAAAATa > 2:793120/1‑90 (MQ=255) ttCATTGGC‑GGGGAGGGTGGCGAGCAGTTTCTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATa < 1:765548/90‑1 (MQ=255) ggggggggggggGAGCGGTTTTTTTGCTGCGTCGTACTGACCCTCTTTTTACAGCACCCGGAGCGGCGCGCCCACAACATACTGGCgggg > 1:23637/1‑88 (MQ=255) | CCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGC‑GGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTT > NZ_CP009273/585309‑585471 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |