Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 523,335 | 0 | T | G | 41.9% | 68.0 / 12.0 | 44 | G25G (GGT→GGG) | rhsD | Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor |
Reads supporting (aligned to +/- strand): ref base T (7/18); new base G (18/0); total (26/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.66e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT > NC_000913/523188‑523483 | ggTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGgtg < 2:509995/149‑1 (MQ=255) tgtgtgTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAGaa > 1:708829/1‑150 (MQ=255) cATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGgcgc < 2:736013/150‑1 (MQ=255) ttttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggc < 2:545807/150‑1 (MQ=255) tgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCCGGGGTCGGCCGGGGGAAGAATTGGCGCGCCCACCcgcg > 1:518060/1‑146 (MQ=255) tgttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtgg < 1:394742/149‑1 (MQ=255) gttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTTTCCCGGGTTCGGCAGGGGTAAGAAATTGCGGCCCCACCcgcgtgggg > 2:94206/1‑148 (MQ=255) gttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTCGGCAGGGGTAAGAAATTGCGCGCCCACCcgcgggggg > 2:422493/1‑148 (MQ=255) aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTc < 1:423537/149‑1 (MQ=255) atatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTTGGGAGGGGTAAGAAATTGCGCGCCCCCCCGCGTGGGGGGGTCGgggtggtccg > 2:679754/1‑146 (MQ=255) atGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcggg < 1:370788/150‑1 (MQ=255) aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCcgggcggg < 1:667101/151‑1 (MQ=255) tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTTGGGAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTTcggggggg > 1:641005/1‑146 (MQ=255) tgggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCAGGGGGAAGAATTTGCGCGCCCCCCCGCGTGGGGGGCTCGGGGTGTCCGGGCGGGGT‑GACtttcgg > 2:311291/2‑145 (MQ=255) cACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCCGGGTGTCCCGGCGGGGG‑GACTTCGGGCAAccc > 1:327074/1‑151 (MQ=255) cACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCACCGGCGTGGCGTGCTCGGTGGGGCCGGGGGGGGT‑GGCTTTGGGCAAacc > 2:473111/1‑151 (MQ=255) cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCAATGTCCAGGGGTCGGCAGGGGTAAGAAATGGGGCGCCCCCCGGCGTGGGGTGCTCGGTGTGTCCGGGGGGGGT‑GACGTCGGGCAAccc > 2:133169/1‑149 (MQ=255) ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCACTGTCCAGGGGTCGGGAGGGGGTAGAATTTGCGCGCCCACCGGCGGGGCGGGGTCGGGGGGTCCCGGCGGGGG‑GACCTCCGGCAAACCgg > 2:465600/1‑149 (MQ=255) tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCCGGGGTTGGCCGGGGTAAGAAATGGGGCGCCCACCCGCGTGGGGGGCTCGGGGTGGCCGGGGGGGGT‑GGATTCGGGCAAacc > 2:635907/1‑146 (MQ=255) tataTGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTTAGAATTGGCGCGCCCACCGGCGGGGCGTGCTCTGGGTGTCCCGGCGGGGT‑GGATTCGGGCAACCCggga > 2:197234/1‑148 (MQ=255) atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCCGGGGTAAGAATTTGGGCGCCCACCCGCGTGGCGGGCTCGGGGGGTCCCGGGGGGGT‑GGCTTCGGGGCACCCCGGaaa > 2:478967/1‑149 (MQ=255) atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAAt < 1:413637/150‑1 (MQ=255) atGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAAt < 2:42167/150‑1 (MQ=255) gAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCg‑ctgct < 1:559488/151‑1 (MQ=255) aaaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCCGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGTGTTCGGGCGGGGTGGCCTTGGGCCACCCCGGAAAATccccgc > 2:219937/1‑145 (MQ=255) aaaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTTGCGCGCCCACCGGCGCGGCGTGCTCGGGGTGGCCGGGGGGGGT‑GAACTCGGGGAAACCGGGAAAACCg‑c > 2:509160/1‑146 (MQ=255) aCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTgggg < 2:224673/151‑1 (MQ=255) cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAAGGGTCGGCGGGGGTAAGAATTGGCGCGCCCACCCGCGTGGCGGGGTCGGTGGGGCCGGGGGGGGT‑GGATTCTGGCAACCCGGGAAATTCG‑CTCCTggggg > 1:705472/1‑150 (MQ=255) cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGaa < 1:569539/151‑1 (MQ=255) cgcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGGGTGGCTGGGGTGTCCCGGCGGGGG‑GACTTCCGGGAACCCGGGGAATCCC‑CCGCTGGGGGCgccggg > 1:613763/1‑146 (MQ=255) cgTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCCCCCGCGTGGCGGGCCCGGTGTGTCCGGGGGGGGT‑GGCCTCGGGCAACCCCGGAAAACCG‑CCGCTGGGGGCGGAg > 2:61051/1‑147 (MQ=255) cgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCCCCGGCGTGGCGTGCTCGGTGTGTTCGGGGGGGGT‑GGCTTCGGGCAACCCCGGAAATCCC‑CTGCCGGGGGCGGAgggg > 2:303561/1‑148 (MQ=255) ggggaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGTGCCGACCGGCGTGGCGTGCTCGGTGTGTCCGGACGGGAT‑GCCTTCGGGCAGCCCGGTAAATCCG‑CTGCTGGGGAAGAAGGTTCTACCCg < 2:192076/148‑1 (MQ=255) ggAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTTCGGGCGGGGT‑GACTTTGGGCAACCCCGGAAAACCG‑CTGCTGGGGGGGAAGGGGCTGccc > 1:589909/1‑150 (MQ=255) tatGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGCGTAAGAATTGGCGCCCCCACCGGCGTGGCGTGCTCGGTGTGGCCGGGCGGGGT‑GACTTCGGGCAACCCCGCAAAACCG‑CTGGTGGGGGGGAAGGGGCCGCCCGGCg > 2:384122/1‑149 (MQ=255) atGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCgag < 2:312354/150‑1 (MQ=255) cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGa < 1:568731/149‑1 (MQ=255) cAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGGAAGAATTGGGGCGGCCACCCGCGTGGCGTGCGCGGTGTGTCCGGGGGGGGG‑TACTTTGGGGAACCCCGGAAATCCGCCTCCTGGGGG‑GGAGGTGCTGTCCCGCGGGGCgg > 2:285285/1‑148 (MQ=255) gTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCtt < 1:50633/151‑1 (MQ=255) ggCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcgc > 2:557184/1‑151 (MQ=255) cGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTg > 1:308134/1‑151 (MQ=255) tCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgccc > 1:156446/1‑151 (MQ=255) cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccg < 1:362501/151‑1 (MQ=255) tGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggccc < 2:768900/150‑1 (MQ=255) ccAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTg < 1:59193/151‑1 (MQ=255) ctgtAAGAATTCGCGCGCCCACCGGCGTGTCGTGCTCGGGGAGTCCGAGCGGGAT‑GACTTCCAGAAACCCCGCAAATCCG‑CTGCTGGGGGTGCAGGTGATGCCCGGCGAGACGGACCTTGCCCTGCCCGGCCCGCCGCCGATCCTtctct > 2:656688/2‑150 (MQ=255) | GGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGAT‑GACTTCGGGCAACCCGGTAAATCCG‑CTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCT > NC_000913/523188‑523483 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |