Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,392,087 | 0 | T | C | 53.8% | 0.6 / 12.1 | 13 | I265V (ATT→GTT) | uspE | universal stress protein UspE |
Reads supporting (aligned to +/- strand): ref base T (4/2); new base C (0/7); total (4/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.10e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAA‑CCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGA > NZ_CP009273/1391998‑1392173 | gCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAA‑CCACAAt > 2:806792/1‑90 (MQ=255) cTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAA‑CCACAATCCCCGCCTGTAAGTGCTCCGc > 2:716814/1‑90 (MQ=255) tgtTGCCGAGGAATGCTGCTAAAATCCCGGTGCGCCCTCCCG‑GCCCAAAACCCACACCCCCCGCCTGTAAGTGCTCCGCCAAATCAGGaa < 1:716814/90‑1 (MQ=255) gtTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAA‑CCACAATCCCCGCCTGTAAGTGCTCCGCCa > 1:1123226/1‑80 (MQ=255) gtTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAA‑CCACAACCCCCGCCTGTAAGTGCTCCGCCa < 2:1123226/80‑1 (MQ=255) aGGAATCCTGCTGAAATCCCGGTGCCCCCTCCCGCCCCCAAAA‑CCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACctct < 1:340644/90‑1 (MQ=255) aaaTCCGGGTGCCCCCCACCGTCCCCAAAC‑CCAAAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCtt < 1:466251/90‑1 (MQ=255) aaaTCCCGGTGCTCCCTCCCGTCCCCAGAC‑CCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCtt < 2:991718/90‑1 (MQ=255) ccGGTGCGTCCTACGGTCCCCAAAC‑CCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCt < 2:566360/90‑1 (MQ=255) gTGCTTCCTCCCGCGCCCACAA‑CCCCATCCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACg < 1:806792/90‑1 (MQ=255) tGCGTCCTACCGTGCCCAGAA‑CCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGt > 2:1035428/1‑90 (MQ=255) caACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGa < 2:745313/90‑1 (MQ=255) aTCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCa > 1:870113/1‑49 (MQ=255) aTCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCa < 2:870113/49‑1 (MQ=255) | GCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAA‑CCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGA > NZ_CP009273/1391998‑1392173 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |