Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 4,034,021 | 0 | G | A | 26.1% | 56.7 / 26.9 | 46 | intergenic (+25/+245) | rrfA/mobB | 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B |
Reads supporting (aligned to +/- strand): ref base G (21/13); new base A (5/7); total (26/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.14e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGG > CP009273/4033901‑4034127 | tCCCCCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGt < 1:370295/138‑1 (MQ=21) aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAt < 1:927012/139‑1 (MQ=34) ccTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATc < 1:532493/139‑1 (MQ=34) tGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCtg < 1:2114756/139‑1 (MQ=34) aCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGACAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCtgtt > 1:3606394/1‑139 (MQ=21) aCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCtgtt > 1:4284366/1‑139 (MQ=34) cccaATGCCGAACTCAGAAGTGAAACGCCGTAGAGCCGATGGTAGAGAGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCACAGTCGGAAGACTGGGCCTTTCGTTTTATCtgttg > 1:3858999/1‑139 (MQ=11) cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCtgttgt > 1:1290191/1‑139 (MQ=34) aTGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgtttg < 1:3461492/139‑1 (MQ=34) tGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgtttg > 1:2152660/1‑138 (MQ=34) gttgcgaagaaacggccgggaggggggggggcaggacgcccgccataaactgccaggaattaattccCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGa < 1:1171139‑M2/72‑1 (MQ=255) tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAAc < 1:1688526/139‑1 (MQ=34) ggagtcaacggcagggaggggggcaggacggtcatccgccataaactgccaggaattaattccCCAGGCATCAAATAAAACGTAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg < 1:1963605‑M2/75‑1 (MQ=255) aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGa > 1:3675388/1‑98 (MQ=37) aaaCGCCGAAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTg > 1:631102/1‑139 (MQ=21) cccggggggtggggggcaggccgcccgccataaactgccaggaattaattccCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg < 1:210197‑M2/75‑1 (MQ=255) tggcgggcaggacgcccgccataaactgccaggaattaattccCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTa > 1:1443842‑M2/44‑130 (MQ=255) tggcgggcaggacgcccgccataaactgccaggaattaattccCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTa > 1:1066486‑M2/44‑130 (MQ=255) ggcgggcaggacgcccgccataaactgccaggaattaattccCCAGGCATCAAATAAAACGAAAGGCTCAGTGGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATcc > 1:273478‑M2/43‑139 (MQ=255) ggcgggcaggacgcccgccataaactgccaggaattaattccCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg > 1:2383214‑M2/43‑127 (MQ=255) aGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGa < 1:2803488/139‑1 (MQ=34) tgGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAACCGAAAGGCTCAGTCGGAAGCCTGGGCCTTTCGTTTTATCTGTTGTTTGCCGGTGAACGCTCTCCTGAGTAGGCCAAATCCGCCGGGAGCgg > 1:222707/1‑139 (MQ=255) cgccataaactgccaggaattaattccCCAGGCAGCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTGTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAttt < 1:4018466‑M2/112‑1 (MQ=255) tctcCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTCGGACAAATCCGCCGGGAGCGGATTTg > 1:3588729/1‑139 (MQ=21) tctcCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTg > 1:1153439/1‑139 (MQ=34) cataaactgccaggaattaattccCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGaa < 1:4247437‑M2/115‑1 (MQ=255) taaactgccaggaattaattccCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGTGAGCGGATTTGAACg < 1:2446393‑M2/117‑1 (MQ=255) gCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg < 1:2439073/92‑1 (MQ=33) gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGa < 1:1857305/87‑1 (MQ=33) ccaggaattaattccCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGaa > 1:1113297‑M2/16‑139 (MQ=255) ggaGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAGCTGTTGTTTGTCGGTGAACGCTCTTCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg < 1:2395713/134‑1 (MQ=12) tAGGGAACTGCCAGGCATAAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAACTGTTGTATGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAAc > 1:8262/1‑139 (MQ=12) aattaattccCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCa < 1:3939809‑M2/127‑1 (MQ=255) gggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg < 1:1942975/83‑1 (MQ=33) ggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg > 1:266144/1‑92 (MQ=33) ggAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACgg > 1:1774963/1‑138 (MQ=32) aaCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCt > 1:4027474/1‑87 (MQ=33) ccAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCt < 1:3037631/45‑1 (MQ=38) ccAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAgg > 1:470386/1‑139 (MQ=32) aGGCATCAAATAAAACGAAATGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg > 1:4112906/1‑83 (MQ=18) gCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTTTCTGTTGTTTGTCGGCGAACGCTCCCCTGAGTa > 1:774729/1‑83 (MQ=11) aaaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTgttgtt > 1:3914588/1‑54 (MQ=35) aaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAt > 1:865582/1‑45 (MQ=38) gAAAAGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAt < 1:162210/37‑1 (MQ=255) aaGGCTCAGTCGGAAAACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg > 1:3952943/1‑55 (MQ=255) aGGCTCAGTCGGAAGACTGGGCCTTTCGATTTATCTGTTGTTTGTCGGTGAACGCTCTCCTg > 1:2469070/1‑62 (MQ=25) | TCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGG > CP009273/4033901‑4034127 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |