Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,685,956 | 0 | G | C | 33.3% | 18.9 / 19.3 | 24 | P331P (CCC→CCG) | PP_4988 | signal transduction histidine kinase CheA/CheY‑like receiver |
Reads supporting (aligned to +/- strand): ref base G (10/6); new base C (5/3); total (15/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTGCGTCAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTT > NC_002947/5685808‑5686099 | aGTGCGTCAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAgg < 1:627095/150‑1 (MQ=255) atcaATCGGCGGCTCTTCGTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTg > 1:84434/1‑150 (MQ=255) atcaATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAgcccttg > 1:114649/1‑144 (MQ=255) aaTCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCg < 2:114649/150‑1 (MQ=255) tCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCggg < 1:293493/150‑1 (MQ=255) gcggcTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGcc > 2:758060/1‑150 (MQ=255) ttCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGa > 1:601962/1‑150 (MQ=255) tccAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTcc > 1:234674/1‑150 (MQ=255) cAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCCACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAg < 1:580959/150‑1 (MQ=255) ttcaaccaacGCTCCACATAGCATGCGCAGCCGACAAACAGCTCGGCCCCGGCCGGCAGCATCAACATGGCCAATACCCCGCGGATCTCCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCg > 2:775377/3‑150 (MQ=255) aaaCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCggg > 1:370519/1‑92 (MQ=37) aaaCAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAACGCCTGCAGCCCTTGCCggg < 2:370519/92‑1 (MQ=37) cAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGa > 2:478364/1‑150 (MQ=255) aGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCg > 1:154524/1‑86 (MQ=255) aGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCg < 2:154524/86‑1 (MQ=255) aGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGat < 2:84434/150‑1 (MQ=255) cgcgGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGc > 1:550766/1‑150 (MQ=255) tCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCa < 1:775377/150‑1 (MQ=255) tggCCCCGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTg > 2:145421/3‑150 (MQ=255) ccccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGAttcattca > 2:327329/1‑98 (MQ=255) ccccGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGAttcattca < 1:327329/98‑1 (MQ=255) cccGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCtt > 2:171656/1‑65 (MQ=25) cccGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCtt < 1:171656/65‑1 (MQ=21) ccGGCAACGCCTGCAGCCCTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGt > 2:237305/1‑150 (MQ=255) cGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGTACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGtt > 2:153999/1‑150 (MQ=255) | AGTGCGTCAGCTCATCAATCGGCGGCTCTTCCTCCAGGTGCTGCAACCAACGCTGCACATAGCTTGCGCAGCCGACAAACAGCTCGGCCACGGCCGGCAGCATCAACATGGCCAATGCCCCGCGGATCTGCTGCAACAGCCCCGGCAAGGGCTGCAGGCGTTGCCGGGGCCAGTCGGACTCCAGGTAATCGCCGATCAGATCCTTGGCCTGCTGCAACACGTTCAGCGATTCATTCAGCACCAGCTGGCGTATCTCGGCCAGGTCCGAGCCCGGCAGGCTGCCTTGGCCGTT > NC_002947/5685808‑5686099 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |