| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 523,335 | 0 | T | G | 54.3% | 35.2 / 51.0 | 36 | 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 (3/13); new base G (19/0); total (23/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.79e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-06 | |||||||||||
TGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTC > NC_000913/523193‑523484 | tgtgtgTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTTGGGAGGGGTAAGaa > 1:526825/1‑150 (MQ=255) gtgtCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGAGATGACTCAGTATGGCGGAACCCTTGTATCAGGTTTAGCAGGGGTAAgca > 2:209145/1‑145 (MQ=255) tttattttTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGGTCGGCAGGGGTAAGAATTGGGGCGccc > 2:129061/1‑151 (MQ=255) tttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtg < 2:443098/150‑1 (MQ=255) gCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGc < 2:794366/150‑1 (MQ=255) tAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCCGGGGGAAGAATTGGGGCGCCCCCCcgcggggcggg > 1:164856/1‑145 (MQ=255) acaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTGCCCTTTTTCAGGGTTCGGGGGGGGTAAGAATTGGCGCGGCCACCGGCGAGGGGGGGTCGgtgt > 2:724158/1‑151 (MQ=255) caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCCATGGCCCGGGTTCGGGAGGGGGGAGAAAACGCGCCCCCACCCGCCGGGCGGGCTc > 1:118746/1‑145 (MQ=255) tGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGGTCGGGGTGTCcggggg > 2:759012/1‑149 (MQ=255) aaCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTTGCGCGCCCACCCGCGGGGCGTGCTCGGGGTGTCcgggcggg > 2:641714/1‑151 (MQ=255) aTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGa < 2:757355/151‑1 (MQ=255) aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGTGGTCCGGGTGTTccgggggg > 2:627439/1‑145 (MQ=255) aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTc < 2:569890/149‑1 (MQ=255) aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTc < 1:209110/149‑1 (MQ=255) aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTC‑gg < 2:242731/151‑1 (MQ=255) aaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTC‑GGGc < 2:611522/151‑1 (MQ=255) ggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGGTGA‑CTTT‑GGGCaa > 2:345885/1‑150 (MQ=255) tttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGGGTGTCCGGGCGGGATGA‑CTTC‑GGGCAAccc > 2:171103/1‑148 (MQ=255) ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTTGGCCGGTGTAAGGAATGGCGCGCCCCCCCGCGTGGGGGGCTCGGGGTGTCCGGGGGGGGGGA‑CTTC‑CTGCAAcccc > 1:665737/1‑147 (MQ=255) tGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGTGGCGGGGTCCGGGTGTCCGGGGGGGGGGA‑CTTC‑GGGGAAACCGGGAAAAcc > 1:274930/1‑151 (MQ=255) ggAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTTGGAAGGGGTAAGAAATGGCGCGCCCACCGGCGTGGGGGGCTCTGTGTGTCCCGGCGGGGTGG‑CTTC‑TGGCAACCCCGCaaaatcgc > 2:553292/1‑143 (MQ=255) ggAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTACCCTTGTGCAGGGTTCGGCAGGGGTAAGAATTGGCGCCCCCACCGCGGGGGCGTGCCCGGGGTGGCCGGGGGGGGGGG‑CTTCGGGGCAACCCGGTaaccccc > 1:229590/1‑143 (MQ=255) cAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCCGGGGTCGGCGGGGGTAAGAATTGGCGCGCCCACCCGCGTGGCGGGCTCGGGGTGTCCCGGCGGGGGGG‑CATC‑GGGCAAACCGGGAAATCCCCCGCTggggg > 1:72569/1‑150 (MQ=255) aGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTTTCCAGGGGTCGGCAGGGGGAAAAAATTGGGCGCCCACCGGCGTGGGGGGCTCGGGGTGGCCGGGGGGGGTGG‑ATTC‑GGGGCACCCCGTAAAACCCCTGGCggggg > 2:434080/1‑149 (MQ=255) cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTTGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGGGTGGTCGGTGTGTTCCGGGGGGGTGACCTTG‑GGGGAACCCGGGAAA‑CCGCTGCTtgggggg > 1:678422/1‑147 (MQ=255) cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCCCGGGGGGTCCGGGCGGGGTGA‑CTTC‑GGGGAACCCCGGAAAACCCCTGCTGGGGGCg > 2:671237/1‑149 (MQ=255) ggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAg < 2:55274/151‑1 (MQ=255) cAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGGGGGGTCGGGGGGTTCGGGGGGGGTGA‑CTTT‑CGGGAACCCCGGAAAACCCCCGCTGGGGGGGGAgggg > 1:333058/1‑145 (MQ=255) cAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgctg < 1:230545/150‑1 (MQ=255) gaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCg < 1:691592/151‑1 (MQ=255) atGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCgaga < 2:674736/151‑1 (MQ=255) cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAc < 1:481294/150‑1 (MQ=255) gTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGGATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCtt < 1:45696/151‑1 (MQ=255) aTGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGGTGA‑CTTC‑GGGCAACCCGGGAAAACCGCTGCTGGGGGGGGAGGTGCTGCCCGGCGAGACGGACCtt > 2:712245/1‑149 (MQ=255) aTGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGGCTCGGGGTGTCCGGGGGGGGTGA‑CTTC‑GGGCAACCCGGGAAAACCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGGGACGGACCtt > 2:240886/1‑149 (MQ=255) tgacccgggTTCGGCAGGCGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCGCGGTGGGTCCGGGCGGGGTGA‑CTTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGGCCTTGCGCTgccccgccc > 1:701979/7‑150 (MQ=255) gtgtAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctctc < 1:373919/151‑1 (MQ=255) | TGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGA‑CTTC‑GGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTC > NC_000913/523193‑523484 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |