| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,040,081 | 0 | T | A | 43.6% | 47.8 / 32.0 | 55 | Y202N (TAT→AAT) | yedY | membrane‑anchored, periplasmic TMAO, DMSO reductase |
| Reads supporting (aligned to +/- strand): ref base T (3/28); new base A (24/0); total (27/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.18e-12 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01 | |||||||||||
GGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGC‑AAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGG > NC_000913/2039936‑2040223 | ggtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATcc < 2:815100/150‑1 (MQ=255) gCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTa < 1:172490/150‑1 (MQ=255) gCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTAt < 2:423096/151‑1 (MQ=255) ccGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATg < 1:426128/151‑1 (MQ=255) gTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTc < 1:1025530/151‑1 (MQ=255) tGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTCGAAgg < 1:1137848/150‑1 (MQ=255) ggTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTCGAAGGAt < 1:794842/151‑1 (MQ=255) ggTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTCGAAGGAt < 1:419775/151‑1 (MQ=255) ttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTCGAAGGATTGc < 2:654270/150‑1 (MQ=255) ttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTCGAAGGATTGc < 2:929039/150‑1 (MQ=255) ccGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTCGAAGGATTGCGt < 1:45438/150‑1 (MQ=255) aTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GGAAAATTCTTTTTTCGCAGGGTTTCGTTTCGGCGaa > 1:1037118/1‑148 (MQ=255) gcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCg < 2:125097/151‑1 (MQ=255) gcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑CGGGGT‑GGAAAAAACCTATGTTGGAGGATTTCGTCTCGCCGGAACAATGCCTcc > 2:263408/1‑150 (MQ=255) gcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑CGGGCT‑GGAGTATTCTTCTGTCCAAGGGCTGCGGTTCGAAGGAGCaaatcctc > 2:59754/1‑143 (MQ=255) cGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GAAAAATTCTTTTTTTGAAGGGTTTTGTCTCGCCGGAGCAATGCATCCGCCcac > 1:194395/1‑150 (MQ=255) gAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑AGGGGT‑GGAAAATCCTTATTTCGGAGGGTTGCGTCTTGACCAAACAAAGCATCCGCTCacac > 1:219185/1‑151 (MQ=255) ccTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GAAAAATGCTTATTTCGGAGGATTGGACTTCTACGAGGAAATGCAtacccccac > 1:1211826/1‑138 (MQ=255) cTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GAAATATCCTTATTTCGGAGGATTGCGTCTCGACGGAGCAATGCATCCGCTCCCACtg > 1:813523/1‑149 (MQ=255) tACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GGAAAATTCTTTTTTCGAAGGATTGCGGCTCGGCGAAACAAAGCATCCCCTcac > 1:1178712/1‑144 (MQ=255) ccAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GGAAAATCCCTATGTCGAAGGAGTGAGTGTCGGCGAAGCAAAGCATCCGCCCACACATATGAc > 2:350610/1‑151 (MQ=255) gCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGTGGAAAAATCTTTTTG‑CGAAGGATTGCGTCTCCGAGGAGCAAAGCATCCGccccca > 2:216629/1‑135 (MQ=255) aaTGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GGAAAATCCTTTTTTTGAAGGAATTCGTCTCGGCGAAGCAAAGCAACCCCCCACACTGATGAcc > 1:110012/1‑147 (MQ=255) cGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GAAATATTCTTTTGTCGGAAGGTTGCGTCTCGACGAAACAAAGCATCCGCTCACACCGATGAACGTAgggg > 1:907480/1‑148 (MQ=255) aaGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGttt < 2:344762/151‑1 (MQ=255) tATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑GGGGCT‑GAAAAATCCTTATTTCGCAGGAATGTGTGTCGACGAAGCAAAACATCCGCTCCCACCTATGAACGCAGGAGttttt > 1:347819/1‑148 (MQ=255) gCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCTGGGC‑‑GAAAAAGCCTTATTTCGGAAGGTTATGGCTCGTCGAAGGAAAGCATCCGCGCACACCGATGTCCGGTGGTGttttttg > 2:1159095/1‑142 (MQ=255) gCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GAAAAATTCCTTTGTTGAAGGGGTGCGGCTCGGCGAAACAAAGTATTCGCCCACACTGATGACCGTGGGGGTTTTTgggcagg > 2:500063/1‑146 (MQ=255) tCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGC‑AAGGCGt < 1:500088/150‑1 (MQ=255) cGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAAAAACCTTATTTCGAAGGGGTGCGTCTCTACGAAGCAATGCATCCCCTCCCACTGATGACCGTAGGGGGTTATGGC‑CAGGGGt > 1:263157/1‑149 (MQ=255) gAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GGAAAATCCTTATGTTGAAGGATTTCGTCTCCCCGAAGCAAAACATCCCCCCACACCGATGACCGTAGGGGGTTTTGGC‑AAAGCg > 1:467372/1‑147 (MQ=255) gAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GAAAAATCCTTTTGTCGAAGGGATTCGGCTCGACGAAGCAAAGCATCCCCCCACACCGATGACCGTAGGGGGTTTTTGG‑AAGGCggt > 1:897810/1‑147 (MQ=255) gAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GAAAAAACCTTATTGCGAAGGAATGCGCGCCTTCGGAGAAATGGATACCTTCCCACCGATGACCGTAGGTGtttcttgcgaagggtta > 2:147432/1‑136 (MQ=255) aaTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑GGGGGT‑GAAAAAAACCTTTGGGGTAAAAATGAGGCTCGACAAAAGAAATCATCCCCTCACACCCATGACCGCAGGGGTTTTTGGG‑AAGGCggt > 1:692484/1‑142 (MQ=255) aCCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAAAATCCTTTTGTCGGAAGATTGCGTCTCGGCGAAACAATGCAACCGCTCACACTGATGAGCGGAGGGGTTTTTGGC‑AAGGCGGTACCCCCCCaaaa > 1:113158/1‑147 (MQ=255) cgggAACAGATGCCAGCCCGGCAGGCCTGCTTGATCGGCGG‑TGGTTT‑AAAATCTCCTTATGTCGAAGTATTGCGTCCGGACAAACCAATCCATCCGCTCCCACTGTTGCCTTTGGTTGTTTATGGA‑AAGGCGTTACCGCCACAAAATGgc < 2:639924/148‑1 (MQ=255) ccGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TTGGGT‑GTAAAAATCTTTTGGCGCAGGAATTCGGTTCCACGAAGAAATACATCCACTTCCCCCGATGGACGTCGGGGTTTTTTTCAAAGGCTTTACCCCCAaaaa > 2:898193/1‑146 (MQ=255) gaaaCAGATGCCAGGCCAGCAGGCCCGTTTTATCGGCGG‑GGGGTT‑GAAATACCCTTATGTCGAAGGATGGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTGGGTTTTTATGGC‑AAGGCTTTACCGCCACAAAATGgcgcg < 2:1080416/149‑1 (MQ=255) gAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GAAAAATCCTTATTTCGGAGGAGTGTGTCTCGACGAAGCAAAGCATCCGCCCCCACCGATGACCGGAGGGGGTTTTGGG‑AAGGGGGTAACGCCACaaaaaagcg > 2:145056/1‑143 (MQ=255) aCAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGGT‑GGAAAATTCTTTTTTCGGAGGGCAGTGCCTCTACGGAGCAAAAGAATCTCTCCCACTGATGCCCCTGGGGGGTTCTTGG‑AAGGGGGTACCGCCCCAAAAAAgcgcggc > 2:835314/1‑148 (MQ=255) agagcccaggccaggAGGACCGTTTTATCGGGGT‑TGGGTT‑GAAATATCCTTATGCGGAGGGATTGCTTCCGGACGAAGCAATGCACCCGCTCACATTGATGCCGGTAGGTTTTTATGGC‑AGGGCTTTACCGCCACAAAATGGCGCGCCGGt < 1:957658/146‑1 (MQ=255) aGATGCCAGGCCAGCAGGCCCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTGGGTGTTTATGGC‑AAGGCGTTACCGCCACAAAATGGCGCGCCGGt < 1:64455/151‑1 (MQ=255) caggccagCAGGACCGCTTTATCGGCGG‑TGGGCG‑GAGAAATCCCTAATTTGAAAGGTTGCGTGTCGACGGAACCAAGCATCCGCCCACCCAGCTTTCCGGGGGGGTTTTTTGG‑CAGGCGATACCGCCCCAAAACATCGCGCCCGGGCg > 2:686505/1‑148 (MQ=255) ggccagCGGGCCCGCTTTAGCGGCGG‑TGGGCT‑GAAATACCCTTATGTCGAAGTATTGCGTCTCGACGAAGCAATGCTCCCGCTCACACTGATGACCGTAGGTGTTTAGGGG‑AGGGTGTTACCGCCACAAAATGGCGCGCCGGTGCGACTg < 2:811677/150‑1 (MQ=255) gcaaagcagcGCCTCTTTTTGGGGGG‑TGGGTT‑AAAATCTCCTTAGGTCGAAGGATTGCGCCCGCTCGAACAAATGATCCCCCTCTCACTGTTGACCGTGGGTGTTCTTGCA‑AAGTGTTTTCCGCCACAAAATGGCGCGCCGGTGCg < 2:525595/142‑1 (MQ=255) ccagCAGGACCGCTTTATCGGCGG‑TGGGGG‑GAAAAATCTCTTTTTTGCAGGATTTCGTCTCGACGCAGCAAATCATCCCCCCACACTAATGAAACTAGGTGGTTTTTTC‑AAGGGGTTACCCCCACAAAAAAGGGCGCCGGGGCggc > 2:956425/1‑144 (MQ=255) ggACCGCTTTATCGGGGG‑GGGGTT‑AAAATATCTTTATGTGGAAGGTTTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAGGACGGTGGGTGTTTAGGGC‑AGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGcc < 2:1037049/150‑1 (MQ=255) tttACCGGCGG‑TGGGCT‑AAAATATCTTTATGTCGAAGGTTTGCGTCTCGACGAAGCAAGGCATCCGCTCACCCTGATGCCCGTAGGTTTTTATGGC‑AAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGaaa < 1:932738/150‑1 (MQ=255) ttttCGGGGG‑TGGGTT‑AAAATATCTTTATGTCGAAGGTTTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGAGGACCGTGGGTTTTTTTGGA‑AAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGACTGTGCCGTGGAAat < 2:46914/147‑1 (MQ=255) cggcgg‑TGGGCT‑GAAAAATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGC‑AAGGCGTTACCGCCACAAAATGGCGCGCCGGGGCGACTGATTGTGCCGTGGAAATATGGc > 2:254117/1‑151 (MQ=255) gggggggggCG‑GAAATTTCTTTATGTGGAGGGATGGCGGCCCGGCGACGCAATGCACCCGCTCAGCGTGGTGGCCGGGGGGGTTGAGGGA‑AGGGCGTTGCCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAatagg < 2:344818/142‑3 (MQ=255) g‑TGGGCT‑AAAATATCTTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCCCCCTGATGACCGTAGGTGTTTATGGC‑AAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTaa < 1:218592/151‑1 (MQ=255) gggCT‑GAAATATCCTTATTCCGAAGGTTTGCGTCTCGACGAAGAAATGCAGCCGCCCACACTGATGACCGTAGGTGTTAAGGGC‑AGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAg < 1:394278/151‑1 (MQ=255) gggCT‑AAAATATCCTTATGTCGAAGGTTTGCGCCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTAGGTGTTTATGGC‑AAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTaaa < 1:390455/150‑1 (MQ=255) ggCT‑GAAATATCCTTATGTCGAAGGTTTGCTTCTCGACGAAGAAATGCCTCCGCCCACACTGATGCCCGTGGGTGTTTATGGC‑AGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAgg < 2:1178636/151‑1 (MQ=255) | GGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGG‑TGGGCT‑GAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGC‑AAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGG > NC_000913/2039936‑2040223 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |