Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 276,800 | 0 | A | C | 57.1% | 12.2 / 14.9 | 14 | Y299D (TAT→GAT) | yagA | IS481‑like element ISEc19 family transposase |
Reads supporting (aligned to +/- strand): ref base A (5/1); new base C (0/8); total (5/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.00e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGA > NZ_CP009273/276716‑276874 | tGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACcgcg < 1:280428/90‑1 (MQ=255) cTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGcc > 1:340872/1‑90 (MQ=255) gACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCat > 2:448020/1‑90 (MQ=255) tCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCa > 1:368574/1‑90 (MQ=255) cccccTTCACCGTTTCCCGGGGGGGTTGTGTTCCCCCTGTCCCGCCGCGCTGACGGCGGATCCCGGGAGCCCGGTACCGCCATATCCAgc < 2:123412/89‑1 (MQ=255) ggggCTTTGTGTTCCCGCTGTCCTCCCGCGCTGACGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:49463/90‑1 (MQ=255) ggggCGTTGTGTTGCCGCTTTCCCGCCGCGCTGACGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:254789/90‑1 (MQ=255) ggggCGTTGTGTTGCCGCTGTCCTGCCGCGCTGCCGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:88614/90‑1 (MQ=255) ggggCGTTGTGTTGCCGCTGTACTGCCGCGCTGCCGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:178490/90‑1 (MQ=255) ggggCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATCCCGCGAGCCGGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:127774/90‑1 (MQ=255) ggggCGTTGTGTTGCCGCTGTACTGCCGCGCTCACGGCTGATCCCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCaa < 1:147714/90‑1 (MQ=255) ttcgCCGCTGCCCTGCCGCGCGGACGGCTGATCCCGCGCGCCCGGTCCCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGAc < 2:302520/87‑1 (MQ=255) gCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCg > 2:84960/1‑90 (MQ=255) gcgcTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGa > 2:454312/1‑90 (MQ=255) | TGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGA > NZ_CP009273/276716‑276874 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |