Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 585,386 | 0 | A | G | 53.8% | 13.8 / 32.0 | 26 | L342P (CTC→CCC) | nfrA | bacteriophage adsorption protein NfrA |
Reads supporting (aligned to +/- strand): ref base A (1/11); new base G (14/0); total (15/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.55e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTT > NZ_CP009273/585298‑585471 | cAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAg < 2:1478097/90‑1 (MQ=255) aaTCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAggggg > 2:1682407/1‑88 (MQ=255) cGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGtttt < 1:675796/90‑1 (MQ=255) cAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGAGGCGGttttt > 1:723977/1‑90 (MQ=255) agagCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGG‑GCAGGTTTTgc > 2:369156/1‑90 (MQ=255) agCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTgcgc < 2:1748608/90‑1 (MQ=255) agCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTgcgc < 2:1331913/90‑1 (MQ=255) cTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTg < 1:1636326/90‑1 (MQ=255) cAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGcgt < 1:1682407/90‑1 (MQ=255) cAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGcgt < 1:852828/90‑1 (MQ=255) tGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGG‑GGAGGTTTTTTGCTGCGTCGGACt > 1:1564264/1‑90 (MQ=255) gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGG‑GGAGGTTTTTTGCTGCGTCGTACTGACCCTCtttt > 1:145805/1‑90 (MQ=255) gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGG‑GCAGGTTTTTTGCTGCGTCGTACTGACCTTCtttt > 1:1367752/1‑90 (MQ=255) gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGG‑GGAGGTTTTTGGGTGGGTCGTAGTTACCTTCcttt > 2:1058317/1‑90 (MQ=255) cTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGc < 1:273786/90‑1 (MQ=255) aTAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGG‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTTAAAGCACCGGTAg > 1:1364234/1‑90 (MQ=255) tAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGG‑GCAGGTTTTTTGCTGCGGCGTACTTACCTTCTTTTTACAACACCCGtctc > 1:448936/1‑87 (MQ=255) cctcAAGCATTTCATTGGCGGGGGGGGGGGGGG‑GCCGGTTTTTTGCCGCGGCCTTCTTACCCTTTTTTTAAAGCACCGGGAGCGTCgggc > 2:491437/1‑90 (MQ=255) cctcAAGCATTTCATTGGCGGGGGGGGGGGGGA‑GGAGTTTTTTTGCTGCGTCGTACTGACCTTTTTTTTAAAACACCCGGAGCGgcgcgc > 1:236714/1‑90 (MQ=255) aaGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCcaac < 2:238645/90‑1 (MQ=255) aGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCcaaca < 2:924964/90‑1 (MQ=255) cATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCaacaac < 1:315104/90‑1 (MQ=255) aTTTCATTGGCGGGGGGGGGGGGGG‑GCAGGTTTTTTGCCCCGCCGGGCTGTCCCTTTTTTATATGCACCGCTGGCGGCGCGCCAacaaca > 2:843388/1‑90 (MQ=255) ttCATTGGCGGGGGGGGGGGGGG‑GGAGGTTTTTCGCCGCGGCGTACTGACCTTCTTTTTAAAGCACCCGGAGGGTCGCGCCCACAACATa > 2:109661/1‑90 (MQ=255) ttCATTGGCGGGGGGGGGGGGGA‑GCAGTTTTTTCGCTGCGTCGTACTTACCTTCTTTTTACAACCCCGGGAGCGGCGCGCCCacaacaaa > 1:238645/1‑88 (MQ=255) ggCGGGGGGGGGGGGGA‑GCAGTTTTTTCGCTGCGTCGTACTGACCTTCTTTTTACAGCACCGGGAGCGGCGCGCCCACAACAAACTGGCg > 1:31617/1‑90 (MQ=255) gggggggggggggg‑GCGGGTTTTTCGCTGCGTCGTACTGACCTTCTTTTTACAGCACCGGGAGGGGGGCGCCAACAACATACTGGCgggg > 2:610901/1‑88 (MQ=255) | CAGCAATCGTGCCAGACGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTT > NZ_CP009273/585298‑585471 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |