Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,612,067 | 0 | A | G | 54.5% | 18.1 / 41.5 | 33 | E161G (GAG→GGG) | yjjJ | type II toxin‑antitoxin system HipA family toxin YjjJ |
Reads supporting (aligned to +/- strand): ref base A (3/12); new base G (18/0); total (21/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.28e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC > NZ_CP009273/4611985‑4612147 | ataCGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGGTGGGCGGGGGGGGGGAt > 1:428185/1‑90 (MQ=255) ctctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAg > 2:315451/1‑90 (MQ=255) ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGGTGGGCGGGGGGGGGGATTATTAGCg > 2:571178/1‑90 (MQ=255) ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGATTATCAggg > 1:356665/1‑88 (MQ=255) ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTGGTCGGGGGGGGGGATTATTAggg > 1:393854/1‑88 (MQ=255) ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTGGTCGGGGGGGGGGATTATGAGCg > 1:60432/1‑90 (MQ=255) tGGCAGGAAGAAGATGTGCTCTCCGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGa < 1:256357/90‑1 (MQ=255) cAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATgg < 2:835568/90‑1 (MQ=255) cAGGAAGAAGATGTGCTCTACCCCCTGATCGTATTTACCGGCGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTACCAGCGATgg < 1:249477/90‑1 (MQ=255) aGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGGTGGGGGGGGGGATTATTAGCGATggg > 1:549439/1‑89 (MQ=255) gaagaaGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTGGGGTGGTGGGGGGGGGGAAGTATCAGAgattgatt > 2:734710/1‑90 (MQ=255) agaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc < 2:56791/90‑1 (MQ=255) gtgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCacaa < 1:593956/90‑1 (MQ=255) tgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc < 1:202228/82‑1 (MQ=255) tgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTAc > 2:202228/1‑82 (MQ=255) gCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATTAGCGATGGGATACTTCacaaca > 2:167589/1‑90 (MQ=255) aCGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAAAACCCTg > 1:754490/1‑90 (MQ=255) cccTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGAATATTAGCGGGGGATTTCTTCACCACACCCCGCgg > 1:143546/1‑90 (MQ=255) aCCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAAtt < 2:487753/90‑1 (MQ=255) aTTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGAATTATCAGCGATGGGTTACTTCACAACACCCCGCGGGAAATTctct > 2:215614/1‑90 (MQ=255) gTGAATATACTGGCGGTTGGTTTGGCGGGGGGGGGAATTATTAGCGATGGGTTTCTGCACAAAACCCCGCGGAAATTCCTCTGGATCaaa > 2:580677/1‑90 (MQ=255) gAATATACTGGCGGTTGGGTGGGCGGGGGGGGGGATTATTAGCGATGGATTAATGCACAACACCCTGCGGGAATTTCTCTTGGTCaaaaa > 1:458140/1‑90 (MQ=255) gAATATACTGGCGGGTGGGTGGGCGGGGGGGGGGATTTTCAGCGATGGATTTCTGGCCAACACCCCGCGGGAAATGCTCTGGGTCaaaaa > 1:352340/1‑90 (MQ=255) aCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAcc < 1:442719/90‑1 (MQ=255) aCTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAcc < 2:599449/90‑1 (MQ=255) tGGCGGTTGGTTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAAAACCCCGCGGGAATTCCTCTGGGTCAAAAACTCACCCa > 2:505102/1‑90 (MQ=255) tGGCGGTTGGTTGGTCGGGGAGGGGGATTATCAGCGATGGATTACTTCACAACACCCTGCGGAAATTCCTCTTGATCAAAAAATCAcccc > 2:722807/1‑89 (MQ=255) gCGGTTGGTTGGTCGGGGGGGGGAATTATCAGGGATGGATTAATGCACAAAACCCCGCGGGAATTCCTCTTGATCAAAAACTCACCCAtt > 1:599449/1‑90 (MQ=255) ggtttgttTGTCGGGGAGGGGAGTTATCAGCGAGGGATTAATGCAGAAAACCCGGCGGGAATTCCTCTGGACAAAAAACGCAGCCactga > 1:136932/1‑86 (MQ=255) ggtttgttGGTCGGGGGGGGGAATTTTCAGCGATGGATTTCTGCACAAAACCCCGCGGGAATTCCTCTGGATCAAAAAATCCACCATTAc > 1:751963/1‑90 (MQ=255) ggttggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAAc < 1:289909/90‑1 (MQ=255) gttggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACa < 2:458140/90‑1 (MQ=255) tggtCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGc < 2:589813/90‑1 (MQ=255) | ATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC > NZ_CP009273/4611985‑4612147 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |