Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 108,789 | 0 | T | C | 57.1% | 8.0 / 23.5 | 21 | S100G (AGC→GGC) | zapD | cell division protein ZapD |
Reads supporting (aligned to +/- strand): ref base T (7/2); new base C (0/12); total (7/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.10e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCCA > NZ_CP009273/108706‑108858 | gATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCccgc < 2:713788/90‑1 (MQ=255) gCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCcgccgc > 1:705582/1‑90 (MQ=255) gccGCTGACGCACCAGAGCAATCAAACGCTCTCCCCGCAAAAATTCCCCGATACGCGGCGGGAAAATTAATCCCCCCCCCGCCGCTTTTa < 1:468241/88‑1 (MQ=255) agagCAATCAAACGTTCTCCCCGCAGAAATTGCCCGATACGGGGGGGGAAAATTAATCCGCCCCCCGCCGCTTTTAACTGCTGAATTAAt < 1:502853/90‑1 (MQ=255) agagCAATCAAACGATCTCCCCGAAAAAATTCCCCGATCCGCGGCGCGGAAATTAATCCCCCCCCCGCCGCTTTTAACTGCTGAATTAAt < 1:572031/90‑1 (MQ=255) aaTCAAACGATCTTCACGCAAAAATTCCCCGATCCGCGGCGCGGAAATTAATCCGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTc < 2:533560/90‑1 (MQ=255) cAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAAt > 2:668104/1‑90 (MQ=255) aaCGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATAc < 2:625687/90‑1 (MQ=255) aTCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGct > 1:765309/1‑90 (MQ=255) tCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGctc < 2:123168/90‑1 (MQ=255) cTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGctct > 1:264527/1‑90 (MQ=255) tCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTgg > 1:738538/1‑90 (MQ=255) aCGCAGAAATTGCCGGATCCGCGGCGCGAAAATTAATACCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTc < 2:476178/90‑1 (MQ=255) aaaaTTCCCCGATACGCGGCGCGAAAATAAATCCCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGc < 1:708471/89‑1 (MQ=255) tccccGATCCGCGGCGCGAAAATTAATCCCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGc < 2:725622/88‑1 (MQ=255) tccccGATCCGCGGCGCGAAAATTAATCCCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGc < 1:241594/88‑1 (MQ=255) gATACGGGGCGAGAAAATTAACCCCCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGcc < 1:447805/90‑1 (MQ=255) ggggcgggaAAATTAATCCGCCCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATcc < 1:331528/88‑1 (MQ=255) cggcgcggAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCCa > 2:68696/1‑90 (MQ=255) ggcgcggAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATcc < 2:143394/88‑1 (MQ=255) ggcgcggAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATcc > 1:143394/1‑88 (MQ=255) | GATGCTCAGTCGCTGACGCACCAGAGCAATCAAACGATCTTCACGCAGAAATTGCCCGATACGCGGCGCGGAAATTAATACGCTCCCCGCCGCTTTTAACTGCTGAATTAATGCTTCAATACGGCTCTGGTCCACGCCAGGCACGCCAATCCA > NZ_CP009273/108706‑108858 |
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: ‑ |