Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,678,039 | 0 | T | G | 44.7% | 71.2 / 13.2 | 48 | N108T (AAC→ACC) | yphE | putative sugar ABC transporter ATPase |
Reads supporting (aligned to +/- strand): ref base T (4/22); new base G (21/0); total (26/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.47e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTG > NC_000913/2677893‑2678188 | gctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGTGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTc < 1:200845/151‑1 (MQ=255) ttttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGCGGCCGA‑GGgttttctcc > 1:416251/1‑145 (MQ=255) tttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCACCGAGGGGGA‑GGgttttctcccc > 2:94522/1‑144 (MQ=255) ttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGGAGAGGGTTTTCccccc > 1:779492/1‑149 (MQ=255) gcgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTg < 2:618014/150‑1 (MQ=255) cgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTc < 1:182735/151‑1 (MQ=255) gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTTACCGAGGCGAA‑GGTTTTCCGCCAcctctagc > 1:534006/1‑141 (MQ=255) cTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCGGGGCAGA‑GGGTTTTCCCCCCTGTCAGcccc > 1:926556/1‑149 (MQ=255) cTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGcc > 1:201243/1‑148 (MQ=255) tAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTc < 1:235005/150‑1 (MQ=255) aGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTTCGCCCCTGTCAGCCCTTcc > 1:954472/1‑150 (MQ=255) gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑Ac < 2:754107/151‑1 (MQ=255) aaaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACt < 1:94949/151‑1 (MQ=255) aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCCAGGCAGA‑GGGTTTTCCCCCCTGTCAGCCCCTCCCCCAGCACTTa > 2:729955/1‑149 (MQ=255) cAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCGCTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGa < 2:735272/149‑1 (MQ=255) aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCCAGGCAGA‑GGGTTTTCGCCACTCTCCGCCCCTCCACCCG‑ACTTAACACCTg > 2:527422/1‑147 (MQ=255) ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATaaa < 2:78465/151‑1 (MQ=255) ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATaaa < 1:1016984/151‑1 (MQ=255) gACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGAACGAGGCAGA‑GGGTTTTCGCCACTTTCAGCCCCTCCCCCCG‑CCCTAAATccctgtaaa > 2:710599/1‑143 (MQ=255) aCTCACGTCAACGCCCAGCGACTGTAAGCACCGTTGGGTATCCTCCGCCTTTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCGCTGACCGAGGCAGA‑GGTTTACCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATaaa < 2:419635/149‑1 (MQ=255) tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGGAGA‑GGGTTTTCGCCACTGTCAGCCCCTCCCCCAG‑ACTTAACTCCTGATAAAcc > 1:100263/1‑149 (MQ=255) aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCGGA‑GGGTTTTCGCCACTGTCAGCCCCTCCACCCG‑ACTTAAATCCTGATAAAcccc > 2:27134/1‑148 (MQ=255) gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCCGGGCAGA‑GGGTTTTTGCCCCTGTCTGCCCTTTCCCCCG‑ACCTAACTCCCGATAAACCgcgcgc > 2:617866/1‑151 (MQ=255) gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCgcgcgc < 1:33637/151‑1 (MQ=255) gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTCCGCCCCTGTCAGCCCTTCCACCCG‑ACTTAACTCCTGATAAAAcccgc > 1:519114/1‑149 (MQ=255) aaCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCCTCCACCCG‑ACTTAACTCCTGATAAACCGCGCGcccc > 2:954936/1‑151 (MQ=255) gCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGA‑GGGTTTTCCCCACTGTCCGCCCCTCCACCAG‑ACTTTAATCCTGATAAACCCCGCGCcccccc > 2:1018840/1‑151 (MQ=255) gcgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCCGGGAGA‑GGGTTTTCGCCACTGTCAGCCCCTCCACCCG‑ACTTAACTCCTGATAAAACCCGCGCCCCCCCCGtt > 2:234477/1‑150 (MQ=255) gcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCAGA‑GGGTTTTTGCCACTGTCAGCCCCTCCCCCAG‑AATTAAATCCTGCTAAAACACGCCCAccccc > 2:695140/1‑144 (MQ=255) cTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTCCGCCACTGTCAGCCCCTCCACCCG‑ACTTAACTCCTTATAAACCGCGCGCACCCCCCCTTCAGCg > 1:335571/1‑151 (MQ=255) gTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCg < 2:87708/149‑1 (MQ=255) cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCCTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACCGCGAGTc > 1:1012882/1‑146 (MQ=255) aTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTa < 2:646547/151‑1 (MQ=255) ccTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGA‑GGGTTTTCGCCCCTGTCCGCCCTTCCACCCG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCCGCGGCACGGCGAGTCCGc > 1:988405/1‑146 (MQ=255) tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTTCGCCCCCTTCAGCCCCTCCACCAG‑ACCTAACTCCCGATAAAACGCGCGCACCCCCCGGTCAGCGGCACGGCGAGGCAGCGTAg > 2:448785/1‑148 (MQ=255) gcgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTc < 1:1011891/151‑1 (MQ=255) cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCg < 2:46209/151‑1 (MQ=255) ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGGGGCAGA‑GGGTTTTCCCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCGAGTCAGCCGAGCTTCGt > 1:314100/1‑150 (MQ=255) ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTc < 2:1039180/151‑1 (MQ=255) aTTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCa < 1:625403/150‑1 (MQ=255) tGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCAGGGCGA‑GGATTTTCCCCACAGGCAGACCTTCCACCCC‑ACTTTACTCCTTGTAAACCCCGCGCACCCCCACGTTAACCGCACGGCGGGGCAGCGTAGCTTCCACACCtt > 1:881194/1‑151 (MQ=255) tcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCCCGCCCACCCCCCGTTCAGCGGCACCGCGAGGCAGCGGAGCTTCGTCACCCTCCAGCCGCGtc > 1:690114/1‑149 (MQ=255) aTGCCGTTGCGGCGGGGCCACTGACCGAGGGAGA‑GGGTTTTCGCCCCTGTCAGCCCTTTCACCAG‑ACTTTAATCCTGATAAAACGCGCGCACCCCCCGTTTAGCGGCACGGCGAGTCAGCGGAGCTTCGTCACCCTCCCGTCGCGCCTCGc > 1:965092/1‑151 (MQ=255) gcggcgGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAg < 2:736126/151‑1 (MQ=255) ggggCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat < 2:954295/150‑1 (MQ=255) cTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCt < 2:686556/150‑1 (MQ=255) gACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc > 1:589851/1‑151 (MQ=255) aGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGc > 1:756177/1‑151 (MQ=255) ggggTTTTCCGCCACTGTCAACCCTTCCGCCAG‑ACTTAACTCCTGAGAAACCGCGCGCGCCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCTCCGCTATCCGGGCGTTc < 2:719311/148‑1 (MQ=255) ggTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt < 1:262467/151‑1 (MQ=255) ttttCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTg > 2:189305/1‑150 (MQ=255) | GCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTG > NC_000913/2677893‑2678188 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |