Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,997,250 | 0 | C | . | 29.2% | 109.7 / 51.1 | 48 | intergenic (‑163/+37) | PP_5252/PP_5253 | hypothetical protein/arylesterase |
Reads supporting (aligned to +/- strand): ref base C (15/19); new base . (7/7); total (22/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.58e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAT > NC_002947/5997109‑5997391 | ggCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGT‑TGTAGGAg < 1:835530/150‑1 (MQ=255) gggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCggg < 2:169673/150‑1 (MQ=255) gggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCggg < 1:72925/150‑1 (MQ=255) gggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTATTCGCGGGTAAACCCGCTCCTACAGACGT‑TGTAGGAGCggg < 1:812429/150‑1 (MQ=39) aCTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTAc < 1:62964/150‑1 (MQ=255) gCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGATTCTCGCTGCCTGGGCAGCACCCCTCTCGGGTTAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAAt < 1:667814/150‑1 (MQ=255) ccGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTc > 1:1000984/1‑150 (MQ=255) ggTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTccc < 2:753265/150‑1 (MQ=255) gCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTa < 2:1000984/150‑1 (MQ=255) aaGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGACTGGGCCTGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTAcc > 2:741174/1‑150 (MQ=255) tGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGCGCGGGGTTACCCGCGAAGTTGGTCCACAATTACCGCTGaa > 2:281378/1‑150 (MQ=38) ttGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGT‑TGTAGGAGCggg > 1:973450/1‑110 (MQ=38) ttGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGT‑TGTAGGAGCggg < 2:973450/110‑1 (MQ=38) aTGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCaa < 2:767362/150‑1 (MQ=255) tgCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTc > 2:401779/1‑150 (MQ=255) gCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTcc < 1:401779/150‑1 (MQ=255) aagaagGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCATCTTCGCGGGTAAACCCGCTCCTACA‑AAGTCTCGGGGAGCGGGTTTACCGGCGAATGTTTCCTAAAATTACTGCTGAAGGAACGCCAGCAAGTCATAATTAAggtg > 1:507260/1‑150 (MQ=37) aagGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAg > 1:15232/1‑133 (MQ=255) aagGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAg < 2:15232/133‑1 (MQ=255) cccTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAg > 2:817134/1‑150 (MQ=255) ccTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGc < 2:90496/106‑1 (MQ=255) ccTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGc > 1:90496/1‑106 (MQ=255) cTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCg < 2:686824/104‑1 (MQ=255) cTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCg > 1:686824/1‑104 (MQ=255) tCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGa < 1:462332/104‑1 (MQ=255) tCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGa > 2:462332/1‑104 (MQ=255) cAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGc < 1:752807/150‑1 (MQ=255) aaGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAgctgc < 1:216202/132‑1 (MQ=255) aaGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAgctgc > 2:216202/1‑132 (MQ=255) ccTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATGACCGCTGCAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATa < 1:244630/150‑1 (MQ=255) gcgGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTc < 1:281378/150‑1 (MQ=255) tAAACCCGCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc < 2:138527/124‑1 (MQ=255) tAAACCCGCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc > 1:138527/1‑124 (MQ=255) cccGCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATg < 2:775601/101‑1 (MQ=37) cccGCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATg > 1:775601/1‑101 (MQ=37) cccGCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACAcc > 2:194058/1‑138 (MQ=255) cccGCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACAcc < 1:194058/138‑1 (MQ=255) ccGCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcgc > 1:884387/1‑150 (MQ=255) ccGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACt > 1:1012604/1‑107 (MQ=255) ccGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACt < 2:1012604/107‑1 (MQ=255) ccGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgcggc > 2:876479/1‑122 (MQ=255) ccGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgcggc < 1:876479/122‑1 (MQ=255) ccGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTc < 2:726752/140‑1 (MQ=255) ccGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTc > 1:726752/1‑140 (MQ=255) gCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCa > 2:1049560/1‑150 (MQ=255) tCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCaa > 2:207605/1‑72 (MQ=21) tCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCaa < 1:207605/72‑1 (MQ=21) tACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCAGCATACACCTTCAGCTCGGCGCCACGAAt > 2:223943/1‑150 (MQ=255) | GGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAT > NC_002947/5997109‑5997391 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |