Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,436,636 | 0 | A | C | 55.0% | 7.9 / 27.3 | 20 | V274G (GTA→GGA) | secY | preprotein translocase subunit SecY |
Reads supporting (aligned to +/- strand): ref base A (9/0); new base C (0/11); total (9/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.95e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTT > NZ_CP009273/3436561‑3436719 | ccAGTCCCCCCCCCGACCCATGACGCGATGGCCGCCGGAAACAAAAAAATACTGGAAGCGAAGATTGCGGGGATT‑CCCCCCGCCATATTc < 1:226850/90‑1 (MQ=255) cccccccAAACCATGACGCGATGGTCGCCGGGAACAAAATAATACTGGAAGCGAAGATTGCCGGGATT‑CCCCCCGCCATATTCACTTTCa < 1:18715/88‑1 (MQ=255) cccccGACCCGTGCCGGGGGGGCCGCGGGAAAAAAAATAATAGGGGAGGGGAAGTTTGCGGGGATC‑CCCCCCGCCATATTCACTTTCAGc < 2:10475/90‑1 (MQ=255) cGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTa > 2:408242/1‑90 (MQ=255) ccgcgAGGGCCGCGGGAAACAAAAAAATACGGGAAGCGAAGTTTGCGGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGc < 1:259156/89‑1 (MQ=255) gggccGCCGGAAAAAAAAAAATACGGGAAGCAAAGATTGCGGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtgtg < 2:13945/86‑1 (MQ=255) ggccGCCGGAAACAAAAAAATACTGGAAGCGAAGATTGCCGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc < 1:221044/87‑1 (MQ=255) ataCTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCt > 1:772386/1‑90 (MQ=255) ggAAGCGAAGATTGCCGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc < 1:336448/61‑1 (MQ=255) ggAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc > 2:336448/1‑61 (MQ=255) ggAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTg > 2:80046/1‑90 (MQ=255) ggAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTg > 2:127508/1‑90 (MQ=255) ggAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTg > 1:742274/1‑90 (MQ=255) ggAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTg > 1:464135/1‑90 (MQ=255) gAAGCGAAG‑TTGGCGGGTTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa < 2:42275/90‑1 (MQ=255) aaGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGAc > 1:687717/1‑90 (MQ=255) gtttGCCGGAATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgc < 1:37628/88‑1 (MQ=255) ttGCCGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt < 1:39746/90‑1 (MQ=255) gggTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt < 1:293115/90‑1 (MQ=255) gggATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt > 2:410895/1‑90 (MQ=255) | CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTT > NZ_CP009273/3436561‑3436719 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |