Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,678,419 | 0 | C | T | 39.5% | 51.5 / 12.3 | 38 | intergenic (‑69/+1) | yhjG/yhjH | putative inner membrane‑anchored periplasmic AsmA family protein/cyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated |
Reads supporting (aligned to +/- strand): ref base C (20/3); new base T (0/15); total (20/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.27e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CAATCAAAATGATCGCCACGACAATCAACAACAAGAAAGCCCCTGAAATCGCAGCGGTTATTTTGCCTGCCTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAAT > NC_000913/3678272‑3678565 | cAATCAAAATGATCGCCACGACAATCAACAACAAGAAAGCCCCTGAAATCGCAGCGGTTATTTTGCCTGCCTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTta < 2:540839/151‑1 (MQ=255) cAAAATGATCGCCACGACAATCAACAACAAGAAAGCCCCTGAAATCGCAGCGGTTATTTTGCCTGCCTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAgc > 2:266867/1‑151 (MQ=255) cAAAATGATCGCCACGACAATCAACAACAAGAAAGCCCCTGAAATCGCAGCGGTTATTTTGCCTGCCTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAg < 2:514863/150‑1 (MQ=255) tCGCCACGACAATCAACAACAAGAAAGCCCCTGAAATCGCAGCGGTTATTTTGCCTGCCTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AAc > 2:514979/1‑151 (MQ=255) cGCCACGACAATCAACAACAAGAAAGCCCCTGAAATCGCAGCGGTTATTTTGCCTGCCTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AAcc > 2:13215/1‑151 (MQ=255) gCCACGACAATCAACAACAAGAAAGCCCCTGAAATCGCAGCGGTTATTTTGCCTGCCTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AAccg > 1:444923/1‑151 (MQ=255) aaaCCCCCTAAAATCGCGGCGTTTTTTTTGCCTCCTTTGCTCATCTTTCGCTTCCTCATATTCTCCCTGTGCCAGTCCTAAAAATAGCCCAGCCA‑GGCGGAAAATGAGGCAGTTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATc < 2:420633/150‑1 (MQ=255) ccTGAAATCGCAGCGGTTATTTTGCCTGCCTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGcc > 1:22614/1‑150 (MQ=255) aaaCCGCGGCGTTTATTTTCCCTGCTTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGGGCCAGCCCAAAAAATATCCCACCCG‑GGGGAAAAATGGGGCGGTTTATAGCGCAAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCggg < 2:313233/149‑1 (MQ=255) aaTCGCAGGGGTTTTTTTGCCTCCCTTCCTCATCGTTCGCTCCCTCATTTTCTCCCTGTCCCAGCCCAAAAGATAGCCCGGCCA‑GGAGAAAAAGGAGGCAGTTTATGGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGt < 2:218273/151‑1 (MQ=255) gcggCGTTTTTTTGGCCTGCTTTGCTCTTCTTTCGCTCCCTCATATTCTCCCTGTGCCGGTCAAAAAATTAGTCCGGCCG‑GGGGAAAAAGGAGCCGGTTTATGGCGCAAA‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTgaaa < 1:364498/148‑1 (MQ=255) tttttttCCCTGCCTTGCTCATCTTTCGCTTCCTCATTTTCTCCCTGTCCCAGTCCTAAAGATAGCCCAGCCG‑GGGGAAAAATGAGGCAGTTTATAGC‑CCAGAAACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTgaaag < 2:512203/142‑1 (MQ=255) ttATTTTGCCTGCCTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAAcc > 1:605519/1‑150 (MQ=255) ttttGCCTGCCTTGCTCATCTTTCGCTCCCTCATATTCTCCCGGTGCCAGTCCAAAAGATAGTCCAGCCA‑GGGGAAAAAGGGGGCAGTTTATGGCGCAAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGc < 2:74279/150‑1 (MQ=255) ttttGCCTGCCTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGcc > 2:204528/1‑151 (MQ=255) ttttCCCTCCCTTCCTCATCGTTCTCTTCCTCTTTTTCTTCCGGTCCCATTCCAAAAGATAGCCCAGCCG‑GGGGAAAAATGGGGCAGTTTATAGCGCAAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAg < 2:491068/149‑1 (MQ=255) tttGCCTGCCTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCt > 2:171085/1‑151 (MQ=255) gcctTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTgcgg > 2:184775/1‑151 (MQ=255) cATCTTTCGCTTCCTCATATTCTCCCTGTCCCATTCCAAAAGATAGCCCACCCG‑GGCGAAAAATGAGGCGGCTTATAGCGCAAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGaa < 2:13040/147‑1 (MQ=255) ttCGTTTCCTCATATTCTTCCTGTGCCGCCCATAAGGATCGCCCAGCGG‑GGCGAAAATGGGGGCAGTTAATAGCGCCAA‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGc < 1:90799/149‑1 (MQ=255) atTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATc > 2:628567/1‑150 (MQ=255) atTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATc > 1:386811/1‑150 (MQ=255) ttttcccTGTCCCCGTCAAAAAGATAGACCAGCCA‑GGCGAAAAATGAGGCAGTTTATAGCGCAAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAAc < 1:580022/141‑1 (MQ=255) gtcccAGTCAAAAAGATAGCCCAGCCG‑GGGGAAAAAGGAGGCAGTTTATGGCGCAAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACg < 2:477311/141‑1 (MQ=255) gggCCAGTCCAAAAGATGTCCCAGCGG‑GGGGAAAAATGGGGCGGTTTAGGGCGCAAA‑ACCCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCAct < 2:576576/147‑1 (MQ=255) cTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGtt > 2:475610/1‑151 (MQ=255) aaaaaATAGCCCAGCCG‑GGCGAAAAAGGAGGCATTTAATAGCCCAAA‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGtt < 2:204558/149‑1 (MQ=255) aaGATACTCC‑GCCAGGGCGGAAAATGGGCCGGTTTAAAGCGACAA‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCt < 1:583921/151‑1 (MQ=255) aGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGcc > 2:625840/1‑151 (MQ=255) gTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGccc > 2:181452/1‑151 (MQ=255) ccccgcca‑gGGGAAAAAGGGGGCAGTTTATGGCCCAAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGcc < 2:117582/145‑1 (MQ=255) cccacccgggggaaaaaTGAGGCAGTTTAAGGCGCCAA‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGCCCt < 1:268961/138‑1 (MQ=255) agcca‑gGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCAGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGCCCTCg > 2:676851/1‑150 (MQ=255) cGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCAc > 1:181980/1‑151 (MQ=255) ggAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCAcc > 2:325427/1‑151 (MQ=255) aaaagggggagtcTTATAGCCCAAA‑ACCCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCAcccc < 1:614318/139‑1 (MQ=255) aaaaTGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCAcccc > 2:598591/1‑151 (MQ=255) aaTGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCcgc > 1:308405/1‑151 (MQ=255) cAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAAt > 2:511059/1‑150 (MQ=255) | CAATCAAAATGATCGCCACGACAATCAACAACAAGAAAGCCCCTGAAATCGCAGCGGTTATTTTGCCTGCCTTGCTCATCGTTCGCTTCCTCATATTCTTCCTGTGCCAGTCCTAAAGATAGTCCAGCCA‑GGCGGAAAATGAGGCAGCTTATAGCGCCAG‑AACCGCCGTATTCAGCGTTTCTATCGGTGCCGGGCGTGAAAGAAACCAGCCTTGTGCGGCGAATGCGGGCGAGTTCTGAACATCACGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAAT > NC_000913/3678272‑3678565 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |