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% | 10.8 / 27.2 | 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 = 5.02e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGT > NZ_CP009273/3436568‑3436714 | cccccccAACCCATGACGCGAGGGCCGCCGGGAACAAAATAATACGGGAAGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCa < 2:77036/88‑1 (MQ=255) cccccccAAACCATGACGCGATGGTCGCCGGGAAAAAAATAATACTGGAAGCAAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCa < 2:257603/88‑1 (MQ=255) gCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAg > 2:43592/1‑90 (MQ=255) cccccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGc > 1:683237/1‑90 (MQ=255) cccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCgg > 1:821240/1‑90 (MQ=255) taaCCATAGCGCGAGGGCCGCGGGGAAAAGAATAATACGGAAGGCAAAGATAGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTaa < 2:275943/89‑1 (MQ=255) aaaCCAGGACGCGATGGCCCCCGGAAAAAAAAAAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTaa < 2:216488/89‑1 (MQ=255) aTGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtg > 2:458517/1‑90 (MQ=255) aTGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtg > 2:30561/1‑90 (MQ=255) gcgGTGGCCGCCGGAAAAAAAATAATATTGGAAGAGAAGTTTGCGGGTATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc < 1:52432/90‑1 (MQ=255) gcgATGGTCGCGGGAAACAAAATAATACTGGAAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc < 2:18143/90‑1 (MQ=255) gcgATGGCCGCGGGAAAAAGAAAAATACGGGAAGCGAACATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc < 2:104669/90‑1 (MQ=255) gcgATGGCCGCGGGAAAAAAAATAATACGGGAAGCAAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc < 2:87962/90‑1 (MQ=255) gcgATGGCCGCCGGAAAAAAAATAATACGGGAAGCAAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc < 2:249566/90‑1 (MQ=255) gcgAGGGCCGCGGGGAAAAAAAAAATACTGGAAGCGAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc < 2:183167/90‑1 (MQ=255) gATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtg > 2:301924/1‑90 (MQ=255) gAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac > 2:24367/1‑90 (MQ=255) ataCTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCt > 1:181085/1‑90 (MQ=255) ggAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGAGCTCTGTGCAGCATAGACACGACGACCTTGCTg > 2:56737/1‑90 (MQ=255) ttGCGGGGATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGATACGACGACCTTGCTGACGTTTCGCGt < 2:574473/90‑1 (MQ=255) | CGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGT > NZ_CP009273/3436568‑3436714 |
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: ‑ |