| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | KT2440 | 3,995,418 | 0 | T | G | 52.4% | 9.9 / 21.0 | 21 | intergenic (+90/‑44) | dmoA_3/allE | Dimethyl‑sulfide monooxygenase/(S)‑ureidoglycine aminohydrolase |
| Reads supporting (aligned to +/- strand): ref base T (0/10); new base G (3/8); total (3/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.14e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-04 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCC > KT2440/3995288‑3995554 | tGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGtt < 6:200928/137‑1 (MQ=255) taaCGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGGGCGGTTCAAgg < 5:490404/136‑1 (MQ=255) cGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCAt < 7:76271/137‑1 (MQ=255) gCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGGGCGGTTCAAGGCATgg < 8:181644/137‑1 (MQ=255) ggTTGGACGTCAGGTCAGCATTCAGTGATCAACTACCGGGTATTCCAACATTTAATTTGTGGAGGCCACGAACTCGACCCTATTCTGTGGCCACTTCCCCTGCGCCGCACCGGCGGGGGCGGTTCAAGGCATGGGga < 5:375014/136‑1 (MQ=255) aCGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGaac < 8:357609/137‑1 (MQ=255) aTCAAATACCGGTTATTCAAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATg < 5:161/137‑1 (MQ=255) aTCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGGGGGGTTTAAGGCATGGGGATAGAACAACAAATAGAGATAGTGCAtt > 6:113662/1‑136 (MQ=255) tCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGt < 4:355575/137‑1 (MQ=255) aaaaTACCTGGTATTCAAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTc < 8:133492/136‑1 (MQ=255) aaaTACCGGTTATTCCAACATTTAATTTGTGGATGCCACGCACTCTACCCTATTCTGTGGCCACATCCCCTGCGCCGCACGGGCGGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCg < 4:408464/137‑1 (MQ=255) aaaTACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCg < 7:374877/137‑1 (MQ=255) aaaTACCGGGTATCCCAACATTTAATTTGTGGAGGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCg < 8:302169/137‑1 (MQ=255) cGGGAATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCAt < 2:147141/137‑1 (MQ=255) accAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACt < 6:167724/136‑1 (MQ=255) ggATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTggg < 6:354756/137‑1 (MQ=255) gcttGCCACGCACGCCACCCTATTCTGTGGCCACATCCGCTGCGCCCCACGGGCTGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTggg < 1:272982/134‑1 (MQ=255) ggggCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGctgac < 7:392656/135‑1 (MQ=255) ccACATCCCCTGCGCCGCCCCGGGCGGGGGGGGTGAACGGATTGGGGGAGAACAAAAAGTCGCGCTCGTGCCTGTGGCACCCTTTATAAATATCCCCCCCCCGTGGGCCCCCGGGTCCAACCGCGCTGCtggctgac > 8:372033/1‑137 (MQ=255) cacaTCCCCTGCGCCGCCCCGGCCGGGGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGAcc < 1:485948/137‑1 (MQ=255) cggccggGGGGGGTGAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTGCATGTCGCACCATTAATAATTAGCCCCCCCCCGTGGGCACCCCGCGCAAACCGACCTGCTGACTGACCGTGGCCTGTTAACCGCAGc > 2:83446/1‑137 (MQ=255) tGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGcc < 4:82434/137‑1 (MQ=255) | TGCCCGAACGGCACGTTGGACGTCAGGTCAGCATTCAGTGATCAAATACCGGGTATTCCAACATTTAATTTGTGGATGCCACGCACTCGACCCTATTCTGTGGCCACATCCCCTGCGCCGCACCGGCCGGTGCGGTTCAAGGCATGGGGAGAGAACAACAAGTCGAGATCGTCCATGTCGAACCATTCATACTTCGCCCCCCACGGTGGGCACCCGGCTCAAACCGAGCTGCTGACTGACCGTGCCATGTTCACCGAAGCCTATGCC > KT2440/3995288‑3995554 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |