Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,596,465 | 0 | T | C | 71.4% | 61.0 / 27.8 | 42 | T366A (ACG→GCG) | livK | L‑leucine/L‑phenylalanine ABC transporter periplasmic binding protein |
Reads supporting (aligned to +/- strand): ref base T (6/6); new base C (16/14); total (22/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTAAATCTTTCACCAGCGCCAGCGGCTCA > NC_000913/3596337‑3596601 | cgTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGt > 1:537944/1‑140 (MQ=255) tAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCg < 1:635559/140‑1 (MQ=255) ggTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGc < 2:34269/140‑1 (MQ=255) gTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCg < 2:537944/140‑1 (MQ=255) tACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGt < 2:424781/140‑1 (MQ=255) cGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGaa > 2:210979/1‑140 (MQ=255) tGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCa < 1:691048/132‑1 (MQ=255) tGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCa > 2:691048/1‑132 (MQ=255) cGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGacac > 1:497965/1‑140 (MQ=255) ttAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCaaaa > 2:329721/1‑140 (MQ=255) aCATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAaaatcaa < 1:281334/140‑1 (MQ=255) cATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAGCCGTCGGCGTGCCACTGGAAGACACCAaaatcaaa < 2:12730/140‑1 (MQ=255) cATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAaaatcaaa < 1:398281/140‑1 (MQ=255) tCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAaaatcaaatc < 1:329721/140‑1 (MQ=255) gctgcAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCt > 1:308279/1‑140 (MQ=255) gctgcAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCt > 1:720903/1‑140 (MQ=255) gctgcAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCt > 2:832398/1‑140 (MQ=255) aGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGAt > 1:388683/1‑140 (MQ=255) aaaaCTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCttttt < 2:388683/140‑1 (MQ=255) gCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATccc > 2:101994/1‑140 (MQ=255) aCCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAAt > 2:891959/1‑140 (MQ=255) tttCTCAACCCGCCCGCATTTTACGGCCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCAc < 2:851559/140‑1 (MQ=255) ttCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACg < 2:720903/140‑1 (MQ=255) ttCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACg > 1:604845/1‑140 (MQ=255) tCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCaa > 2:298074/1‑134 (MQ=255) tCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCaa < 1:298074/134‑1 (MQ=255) tCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACgg < 1:210979/140‑1 (MQ=255) cccgcATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACTGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCAcc > 2:350764/1‑140 (MQ=255) gcATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGtt < 2:497965/140‑1 (MQ=255) gcATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGtt < 1:625323/140‑1 (MQ=255) cGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTa < 1:454523/140‑1 (MQ=255) gTGGGATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTAAATCtt > 2:870870/1‑140 (MQ=255) gggATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTaaa > 1:525294/1‑134 (MQ=255) gggATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTaaa < 2:525294/134‑1 (MQ=255) gggATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTaaa < 2:842937/134‑1 (MQ=255) gggATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTaaa > 1:842937/1‑134 (MQ=255) gggATGATCACTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTAAATCTTTc > 2:368378/1‑140 (MQ=255) aCTTGGCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTAAATCTTTCAccagcgcc > 1:432144/1‑140 (MQ=255) ttGGCTGCCGCGGATGAACCGTCTGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTAAATCTTTCAccagcgccag > 2:273/1‑140 (MQ=255) ggCTGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTAAATCTTTCAccagcgccagcg > 1:233541/1‑140 (MQ=255) tGCCGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTAAATCTTTCACCAGCGCCAGCGGCt > 2:216506/1‑140 (MQ=255) ccGCGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTAAATCTTTCACCAGCGCCAGCGGCTCa < 1:870870/140‑1 (MQ=255) | CGTAGGTACTGCCCAGCGTGACGCCGTTAAACATCTGCTGCAAGAAATACAAAAACTGCTCAGACATAAGGTAACCTTTCTAAACCCGCCCGCATTTTACGGGCGGTGGGATGATCACTTGGCTGCCGTGGATGAACCGTCGGCGTGCCACTGGAAGACACCAAAATCAAATCCCTTAAGATCGCCTTTTTCATCCCAGTTCAGCGGCCCAATCACGGTGTTTGCACCGTTAGCTTTTAAATCTTTCACCAGCGCCAGCGGCTCA > NC_000913/3596337‑3596601 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |