| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,040,075 | 0 | C | G | 49.1% | 46.0 / 30.0 | 57 | L200V (CTG→GTG) | yedY | membrane‑anchored, periplasmic TMAO, DMSO reductase |
| Reads supporting (aligned to +/- strand): ref base C (7/22); new base G (28/0); total (35/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.55e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGG > NC_000913/2039930‑2040223 | gATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGaaa < 2:808530/151‑1 (MQ=255) ggtggtGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAata < 1:610294/150‑1 (MQ=255) gtggtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCCCCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATc < 2:13662/151‑1 (MQ=255) ggtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTTGGGGGAAatat > 2:18522/1‑148 (MQ=255) ggtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCGGGAATATCC‑t > 1:989852/1‑151 (MQ=255) gtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑tt < 1:239791/151‑1 (MQ=255) tgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGaa > 1:417812/1‑142 (MQ=255) gCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑TTa < 1:428616/150‑1 (MQ=255) gCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGAACGCTTTATCGGCGGTGGGCTGAAATATCC‑TTa < 1:864680/150‑1 (MQ=255) cGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCC‑tttttt > 1:308973/1‑147 (MQ=255) cGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGGTGAAATAACC‑TTTTGt > 2:137083/1‑151 (MQ=255) gTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATAC‑CTATGt > 2:1136961/1‑150 (MQ=255) tGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑TTATGTCGAAgg < 1:1025191/150‑1 (MQ=255) ttttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAAAATCC‑TTTTGTCGAAGGAtt > 2:1143487/1‑150 (MQ=255) cTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATAC‑TTTTTTTGGAGGAATGCg > 2:615978/1‑146 (MQ=255) caAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCC‑TTTTGTCGAAGGATTTCGTCTcgac > 2:966123/1‑148 (MQ=255) tGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATTC‑TTTTTTCGGAGGGTTGGGTCTCGACGAAACaa > 2:319415/1‑150 (MQ=255) ggCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCa < 1:345629/150‑1 (MQ=255) ggCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCAt < 1:292463/151‑1 (MQ=255) aaCCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATAC‑CTATTTCGGAGGGGTTCGGCTCCGCGAAGCAAAGCATcc > 1:52277/1‑142 (MQ=255) aaCCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACt < 2:160376/151‑1 (MQ=255) cTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGCAAAAACC‑TTTTTTGGGAGGAGTGCGTCTcggcgcaaacaagcagtc > 1:882316/1‑126 (MQ=255) tACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCC‑CTATGTCGAAGGGTTGTGTCTCCGCGAAGCAATGCATCCGCTCACACtgat > 2:413155/1‑150 (MQ=255) aCCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACtgatg < 2:648100/150‑1 (MQ=255) ccAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGCTGAAAAATAC‑TTATGTCGAAGGATTGCGTCCCGACGAAGCAAAGCATCCGCTCccaccg > 1:460749/1‑144 (MQ=255) aaTGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCC‑TTTTGGTGGAGGAGTGCGTCTCGACGAAGCAATGCATCCGCCCCCACCGGTGAcc > 1:894898/1‑147 (MQ=255) aTGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAgg < 1:946765/151‑1 (MQ=255) tGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCGGGAAAATCCTTTTTTTGGAAAGATT‑AGGCTCGACGACGCAATTCATCCGccccccctggtggcc > 2:37344/1‑131 (MQ=255) tGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGGTGGAAAATTC‑CTTTTTCGAAGGAATGTGTCTCGAAGGAGCAAAGCATTCCCCCCCACTGATGACCGTCgg > 1:163616/1‑150 (MQ=255) gCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGGTGAAATATTC‑CCCTGTCGAAAGAATTTGTCTCGGCGGAGGAAAACATACCCTCACACTCATGACCggatg > 1:1150967/1‑145 (MQ=255) aaGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGttt < 2:422472/151‑1 (MQ=255) aGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATAC‑TTTTTTTGGAGGATTTTGTCTCGACGGAGCAAAGCATCCGCTCACACCGATGACCGTAgggggtt > 2:163881/1‑145 (MQ=255) aTGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATg < 2:269200/150‑1 (MQ=255) aTGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATg < 1:972925/150‑1 (MQ=255) tCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCC‑TTATTTCGGAGGGTTTGGTCTCGGCGAAGCAATGCAACCCCTCCCACTGATGGACGGAGGGGttttt > 1:1050281/1‑144 (MQ=255) tCGAAACATTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCC‑TTATGTCGAAGGGTTGCGTCTCGACGAAGCAAAGCATCCGCTCACACTTATGACCGCAGGTGTTTTTTGGAAGGCggt > 2:528127/1‑149 (MQ=255) cGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCCGGAAAATTC‑TTATGTCGAAGGATTGCGTCTCGACGAAACAAAGCATCCGCCCACACTGATGACCGTAGGGGGTTTTTGCAAGGCggt > 2:933233/1‑148 (MQ=255) gAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCC‑TTTTGTCGAAGGGATTCGGCTCCGCGAAGCAAAGCATCCCCTCCCACTGATGACCGTAGGGGGTTATGGCAAGGCGGTa > 1:867869/1‑150 (MQ=255) aaCAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCC‑TTATGTCGGAGGATTGCGTCTCGACGAAGCAAAGCATCCGCCCACACCGATGACCGTAGGGGTTTTTTGGAAGGGGTTAccc > 2:853284/1‑150 (MQ=255) aCAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCC‑TTTTGTCGGAGGGTTTCGTCTCGACGAAGCAAAGCATCCGCCCCCACTGATGACCGTAGGGGGTTTTTGGAAGGCGGTACCGc > 2:1155736/1‑151 (MQ=255) aaTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGAAAAAACC‑TTATGTCGAAGGGTTGCGGCTCGGCGCAACAAAGCATCCGCTCACACTCATGAGCGTAGGGGGTTATTGCAAGGCGTTAcacccc > 1:701632/1‑146 (MQ=255) ttATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTAGCGGCGGTGGGCTAAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCCCACTGATGACCGTAGGTTTTTAGGGAAAGGCGTTACCGCcaca < 1:357772/150‑1 (MQ=255) aTGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATTC‑TTTTTTCGAAGGGGTTTGGCTCGCCGAAGGAAAACATACGCCCCCACTGATGACCGGAGGGGGTTTTGGGAAGGGGGTACCGCCACaaa > 1:950368/1‑150 (MQ=255) aTGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCC‑CTTTGTCGGAGGAGTTCGTCTCGCCGAAGCAAAGCATCCGCTCACACTGATGACCGTGGGGGGTTTTGGCAAGAGGGTACCGCCCCaaa > 2:861886/1‑150 (MQ=255) gCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGATGGGTGGAATATCC‑TTTTGTCGGAGGGTTGCGTCTCGACGAAGCAAAGGAACCCCTCACACTGATGAACGCAGGGGGTTATTGCCAGGGGGTAACCCCACaaaaa > 2:355241/1‑149 (MQ=255) aCAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCC‑TTATTTCGAAGGAGTGCGTCTCCACGGAGCAATGCATCCGCCCACCCTGATGACCGTAGGGGTTTTTGGCAAGGGGTTACCCCCCCCAAATGGGGCGcc > 1:784004/1‑150 (MQ=255) cAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCC‑TTTTGTCGAAGGGGTTTGTCTCGCCGCAGCAAAGCATCCGCTCACACTGATGACCGTAGGGGTTTTTGGGAAGGGGGTACCCCCCCAAAAAGGCGCGCCgg > 2:777851/1‑151 (MQ=255) aTGCCAGGCCAGCAGGACCGTTTTATCGGCGGTGGGCTAAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTGGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTg < 2:810092/150‑1 (MQ=255) aTGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCC‑TTATTTCGAAAGATTGCGTCTCGACGAAACAAAGCATACGCGCACACTGATGGCCGGAGGGGTTTATTGCAAGGGGTTTCCGCCACAAAAAAGCGCGCCgg > 2:196736/1‑148 (MQ=255) tGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCC‑CTTTGTCGAAGGGTTGGGTATCGACGCAACAAAGCATCCGCTCACACCGATGGGCGGAGGGGGTTTTTGCAAGGCGGTAACCCCCCAAAAAAGCGCGCCgggggg > 2:756818/1‑147 (MQ=255) ggcccGCTTTATGGGCGGGGGGCTAAAATTTCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGAAAGGCATCCGCCCACACTGATGACCGTAGGTTTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCg < 1:76483/148‑1 (MQ=255) gCTTTACCGGCGGTGGTCAAAAATATCC‑TTATGCCGAAGGATTGCGTCTCGACGAAGCAATGCTTCCGCTCCCCCTGTTGACCGTGGGTGTTTATGGAAGGGCGTTGCCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGa < 1:744337/150‑1 (MQ=255) cTTTATCGGCGGGGGGCTGAAATATCC‑TTTTTTCGAAGGGGTTTGTCTCGGCGAAGCAATGCATCCGCTCACACTGATGAACGGAGGGGTTTTTTGGAAGGGGGTAACGCCCCAAAAAAGCGCGCCGGGGCGGCTGGTTGTGGCggggga > 1:644582/1‑148 (MQ=255) aTCGGCGGTGGGCTAAAATATCC‑TTATGCCGAAGGATTGCGTCTGGACGAAGCAATGCATCCCCTCACACTGATGCCGGTGGGTGTTTATGGCAGGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAatat < 1:688488/150‑1 (MQ=255) ggcggTGGGCTAAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAatat < 2:312711/147‑1 (MQ=255) gggCTAAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACGGTGGGTTTTTATGGAAGGGGGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTaaa < 1:966092/150‑1 (MQ=255) ggCTGAAATATCC‑TTATGTCGAAGGCTTGCGTCTCGACGAAGCAATGCATCCGCTCAAACTGATGACCGTAGGGGTTTATGGCACGGCGTTACCGCCACAAAATGGCGCGCCGGCGCGACTGATTGTGCCGGGGAAATATGGCTTTAAAgg > 2:896391/1‑151 (MQ=255) | GATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCC‑TTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGG > NC_000913/2039930‑2040223 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |