Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,679,388 | 0 | T | C | 56.5% | 10.1 / 28.0 | 23 | E1049G (GAG→GGG) | bcsC | cellulose synthase complex outer membrane protein BcsC |
Reads supporting (aligned to +/- strand): ref base T (7/3); new base C (0/13); total (7/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.89e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAC > NZ_CP009273/3679302‑3679466 | gCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTccc > 1:28384/1‑90 (MQ=255) gCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTccc > 2:52350/1‑90 (MQ=255) gCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTccc > 1:429173/1‑90 (MQ=255) aGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACg < 1:356354/90‑1 (MQ=255) ggtcggGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTAAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCa < 2:74361/90‑1 (MQ=255) aTTCATCAGCGGATAACGCGGCATGGTTTTGGTTCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTgcg < 2:153737/90‑1 (MQ=255) catcaGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTc > 1:307445/1‑90 (MQ=255) cGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCCCCCAGCCCCCACGCCCATTTTTCCGTGCGCTCCCGCca < 1:140430/90‑1 (MQ=255) cgcgGCTTGTTTTTGGTGGGTGAAGGCGCCCCCGAGCCAGCCGCCCCCAGCCCCCACGACCATTTTCCCGTGCGCCCCCGCCACATCAcc < 2:179611/90‑1 (MQ=255) gcttGGTTTTGGTGCGTGAAGGCCCCCCCGACCCAGCCGCCCCCACCCCCCCCGCCCATTTTTCCGTCCCCCCCCGCCACATCACCGGTa < 2:429173/87‑1 (MQ=255) ttttGGGGCGTAACTGCGCCCCCGGCCCAGCCCCCCCCACCCCCCACGACCAGTTTCCCGGGCGCTCCCGCCACATCACCGGTATGGCaa < 2:322812/90‑1 (MQ=255) tttgGGTGCGTAAATGCCCCCCCGCCCCAGCCCCCCCCCCCCCCCACCACCTTTTTCCCGTGCCCCCCCGCCACATCACCGGTATGGCaa < 2:424841/90‑1 (MQ=255) ggTGCGTAAAGGCCCCCCCGACCCAGCCCCCCCCCCCCCCCACGCCCATTTTCCCGGCCCCCCCCGCCACATCACCGGTATGGCAAACGa < 2:28384/90‑1 (MQ=255) tgggTGAATGCGCCCACGACCCAGCCGCCCCCAGCCCCCACGCCCTTTTTCCGGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACa < 1:356185/87‑1 (MQ=255) tGCGTGAATGCGCCCACGAGCCAGCCGCCCCCACCCCCCACGCCCATTTTTCCGTGCGCCCCCGCCACATCACCGGTATGGCAAACGACa < 2:202654/90‑1 (MQ=255) tGCGTGAATGCGCCCACGAGCCAGACGCCCCCACCCCCCACGCCCATTTTTCCGTGCCCTCCCGCCACATCACCGGTATGGCAAACGACa < 1:52350/90‑1 (MQ=255) gTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGt > 1:335930/1‑90 (MQ=255) gTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGt > 1:97379/1‑90 (MQ=255) ccACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCgg > 1:440127/1‑90 (MQ=255) ccACGACCCAGCCCCCCCCCGCCCCCCCGCCCTTTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCgg < 2:278796/90‑1 (MQ=255) ccGACCCAGCCCCCCCCCGCCCCCCCCACCATTTTCCCGGCCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGAc < 2:307445/89‑1 (MQ=255) ccGCCCCCAGCCCCCACCCCCTTTTTTCCGCGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAc < 2:243576/89‑1 (MQ=255) ccGCCCCCACCCCCCCCGACCATTTTTCCGTGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAc < 2:149561/89‑1 (MQ=255) | GCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAC > NZ_CP009273/3679302‑3679466 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |