Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,927,661 | 0 | A | T | 48.5% | 27.3 / 20.0 | 33 | S401T (TCA→ACA) | edd | phosphogluconate dehydratase |
Reads supporting (aligned to +/- strand): ref base A (12/5); new base T (0/16); total (12/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.26e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCAC > NZ_CP009273/1927582‑1927745 | aGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGcc < 2:258613/90‑1 (MQ=255) ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGcc > 1:148769/1‑84 (MQ=255) ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGcc < 2:148769/84‑1 (MQ=255) ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc > 2:365353/1‑90 (MQ=255) ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc > 1:131274/1‑90 (MQ=255) ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc > 1:649879/1‑90 (MQ=255) ccGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTCGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTcc > 1:770844/1‑90 (MQ=255) ttAACCCTTTTGCCCCCCAATGATGAGAAAAAGTTTGTTCAAAGGAAGCGATCACATTGCTGTCGAGTGTTTTTCCCGCCCCTCCCCGcc < 1:785391/90‑1 (MQ=255) tcTTTTGCCCCCCCAGGAGGAAAAAAAGGTTGTTCAAAGGAAGCGATCACTTTGCTGTGGAGTTATTTTCCCGCCCCTTCCCGccagtcc < 1:349223/89‑1 (MQ=255) gTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCa > 1:22949/1‑90 (MQ=255) accaTGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCatt > 2:671944/1‑90 (MQ=255) gagaAAGTTTGTTCGAAGAAGGCGATCACTTTGCTGTCGAGTTTTTTTTCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt < 1:365353/90‑1 (MQ=255) gagaAAGGTTGTTCAAAGGACGGGATCCTTTTCCTTTGGTGTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt < 1:398247/90‑1 (MQ=255) gaaaaaGGTTGTTCGAAGGAGGCAATCAATTTGTTGTCGTGTGTTTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCAt < 1:705643/87‑1 (MQ=255) gaAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATgg > 1:339866/1‑90 (MQ=255) aagggTTGTTCGAAGGAAGCGATCCCATTGCTGTCGAGTGTTTTTTCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGt < 2:495182/87‑1 (MQ=255) gtttgttCGAAGGAAGCGATCACTTTGCTGTCGAGTTTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCa < 1:58895/88‑1 (MQ=255) ggTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCa < 2:339866/90‑1 (MQ=255) tttgttCGAAGGAAGCGATCACTTTGCTGTGGAGTGTTTTTCCCGCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCaa < 2:649879/89‑1 (MQ=255) gtttttCGAAGGAAGGCATCACATTGCTGGCGGGTTTTTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCaa < 2:312972/90‑1 (MQ=255) tgttCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAgg > 2:558785/1‑90 (MQ=255) ttCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGt < 2:321488/90‑1 (MQ=255) gCGATCCCATTGCTGTGGAGTGTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac < 2:784288/90‑1 (MQ=255) gCGATCCAATTGCTGTCGAGTGTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac < 2:461212/90‑1 (MQ=255) gCGATCACTTTGCTGTGGAGTTTTTTTCCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac < 1:252355/90‑1 (MQ=255) gCGATCACATTGCTGTCGAGTGTTTTTTCCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac < 2:22949/90‑1 (MQ=255) gCGAACACTTTGCTGTGGAGTTTTTTTACCCCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac < 2:409951/90‑1 (MQ=255) ttGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCaaaa > 2:625413/1‑90 (MQ=255) gCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac < 2:44529/79‑1 (MQ=255) gCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGagac > 1:44529/1‑79 (MQ=255) gAGTGTTTTTTCCGCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCAc < 2:304653/90‑1 (MQ=255) gAGTGTTTTTCCCCCCCCTCCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCAc < 1:792007/90‑1 (MQ=255) gAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCAc > 1:754888/1‑90 (MQ=255) | AGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCACATTGCTGTCGAGTGATTTTTCCGCCCCTTCCCGCCAGTCCAGTTCACCATTATTCAGCCATGGTTCAAGGGTATAACGAGACAGACCAAAACCTGCCAC > NZ_CP009273/1927582‑1927745 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |