| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 730,748 | 0 | T | C | 44.7% | 53.4 / 13.4 | 47 | I389T (ATC→ACC) | rhsC | Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor |
| Reads supporting (aligned to +/- strand): ref base T (17/9); new base C (0/21); total (17/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.34e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGC > NC_000913/730599‑730895 | tgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatc > 1:436585/1‑151 (MQ=38) cGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAcc > 2:136624/1‑151 (MQ=40) gATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccgc < 1:417769/151‑1 (MQ=40) aTGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccgcc < 2:828374/151‑1 (MQ=40) gcACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTc > 1:731442/1‑151 (MQ=40) cGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGc > 1:313566/1‑151 (MQ=255) tCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGcac > 2:217468/1‑151 (MQ=255) acGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGc > 2:89998/1‑151 (MQ=255) ccGACCGGAAATCTGTTCCCGTTACGACAGCGGGGGGCGGGTGAAAGAACAGCAAACCCC‑GGCAGGCTAAATTTACACGTACCAGTATGAGAAGGCCCGCACCCCCATCACCGCCAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGt < 2:544190/151‑1 (MQ=255) ccGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGt > 2:918680/1‑151 (MQ=255) gACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTg > 1:904194/1‑150 (MQ=255) aCCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGa > 2:193416/1‑150 (MQ=255) cggaaaaTCTTTTCCCGTTACGACAGCGAGGGGGGGGTGACAGAACAGCAAACCCC‑GGCAGGCTTAATTTACACGTATCAGTATGAGAAGGACCGCATCCCCATCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGa < 1:28075/145‑1 (MQ=14) ccGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGa > 1:305917/1‑149 (MQ=255) ccGAAAATCTGTTACCGTTACGACAGCGACGGGGGGGTGACAGAACAGCTAACCCC‑GGCGGGCTTAAGTTACACGTATCAGTATGAGAAGGCCCGCATCACCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAg < 1:749503/151‑1 (MQ=14) aattCTTTTCCCGTTACGACAGCGGGGGGGGGGTGACAAAACAGCAAACCCC‑GGCGGGCTTAAGTTACACGTATCAGTAAGAGAGGGCCCGCACCCCCACCCCCGACAGCTTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCgg < 2:644902/148‑1 (MQ=255) aaaTCTGTTCCCGTTACGACAGCGACGGGCGGGTGACAGAACAGCAAACCCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGCCCGCATCCCCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCg < 1:170236/150‑1 (MQ=255) cttttACCGTTCCTAAAGCGACGGGCGGGTGAAAGAACAGATAACCCC‑GGCGGGCTAAAGTTACACGTATCAGTAGGAGAGGGCCCTCACCCCCACCCCGGACGGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCg < 2:497220/143‑1 (MQ=255) ttcccGTTACGACGGCGAGGGGCGGGTGAAAAAACAGCAAACCCC‑GGCGGGTTTAAGTTACACTTACCAGTATGAGAAGGACCGCATCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGaa < 1:441159/147‑1 (MQ=255) ttcccGTTACGACAGCGACGGGCGGGTGACAGAACAGCAAACCCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGaa < 1:312559/147‑1 (MQ=14) ccGTTACGACAGCGACGGGGGGGTGACAAAACAGCAAACCCC‑GGCAGGTTTAATTTCCACGTATCAGTATGAGAAGGACCCCATCACCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCg < 1:556781/150‑1 (MQ=255) ttACGACGGCGCCGGGGGGGTGACAAAACAGCAAACCCC‑GGCAGGTTTAAGTTACACGTATCAGTATGAGAAGGCCCGCATCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCggg < 1:961504/149‑1 (MQ=255) aCAGCGAGGGGGGGGTGAAAGAACAGCTAACCCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGggtg < 1:1043085/146‑1 (MQ=255) aCAGCGACGGGGGGGTGACAGAACAGCTAAACCC‑GGCAGGTTTAAGTTACACGTATCAGTATGAGAAGGACCGCACCACCATCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGa < 1:514189/150‑1 (MQ=255) aCAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGa > 2:543873/1‑150 (MQ=255) cAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGaa > 2:752796/1‑150 (MQ=255) gaggggggggTGACAGAACAGCAAACCCC‑GGCAGGTTTAATTTACACGTATCAGTATGAGAAGGCCCGCATCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGggtggt < 1:876011/140‑1 (MQ=14) gggggggTGAAAAAACAGCTAACCCC‑GGCGGGCTTAATTAACCGTTATCGGTAGGAGAGGCCCCGCACCCCCACCCCCGACAGCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAgg < 2:465020/148‑1 (MQ=9) ggggGAAAGAACAGCAAACCCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGCCCGCACCCCCATCCCCGACGGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACAcg < 2:188657/150‑1 (MQ=14) gggTGAAAGAACAGCTAACCCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGCCCGCATCCCCACCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACAcg < 1:627356/150‑1 (MQ=255) gggTGAAAGAACAGCAAACCCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGCCCGCACCCCCATCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACAcg < 1:688029/150‑1 (MQ=14) ggTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACAcgcg > 1:436417/1‑151 (MQ=255) ggTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACAcgcg > 1:359152/1‑151 (MQ=255) aGAACAGCAAACCCC‑GGCGGATTTAGGTAACACGTATCAGTAGGAGAAGGCCCCCATCCCCACCCCCGACGGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACg < 1:31910/149‑1 (MQ=14) aacccc‑GGCAGGTTTAAGTTACACGTATCAGTATGAGAAGGCCCCCTCCCCCACCCCCGACAGCCTAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCa < 1:1061916/146‑1 (MQ=14) acccc‑GGCGGGCTTAATTAACACGTATAAGTATGAGAAGGACCGCACCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACgg < 2:994449/138‑1 (MQ=14) aaCCCGGGCAGGTTAAAGTTACACGTATAAGTAGGAGAAGGCCCGCACCCCCACCCCCAACAGCTTAAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCAc < 1:399064/150‑1 (MQ=14) c‑GGCTGGCTTATTTTAAACGTAAAAGAATGAGAGGGACCGCACCCCCACCACCGACAGCCTGATCCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCag < 1:41139/149‑1 (MQ=14) cAGGTTTAAGTTACACGTACAAGTAGGAGAAGGCCCGCATCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTc < 2:313396/150‑1 (MQ=14) cAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTc > 1:739741/1‑150 (MQ=255) aGGTTTAAGTTACACGTATCAGTAGGAGAAGGACCGCATCCCCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTc < 1:588547/149‑1 (MQ=255) ggCTTAAGTTACACGTATCAGTATGAGAAGGCCCGCATCACCATCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAg < 1:122599/150‑1 (MQ=255) aGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGAcg > 2:455395/1‑151 (MQ=255) aGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGAcg > 1:170346/1‑151 (MQ=255) acGTATCAGTATGAGAAGGACCGCATCACCACCCCCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGt < 1:56383/149‑1 (MQ=255) ggagaAGGCCCGCACCCCCCCCACCGACGGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCt < 2:323192/146‑1 (MQ=14) gagaAGGCCCGCACCCCCACCCCCGACAGCCGGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAg < 1:854708/149‑1 (MQ=14) gaAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGc > 1:294693/1‑150 (MQ=255) cccacccccatcaccGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACg < 2:829655/143‑1 (MQ=255) catccccaccacCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGAt < 1:414156/151‑1 (MQ=255) atcaccatcaccGACAGCCTGAACCGCCGTGAAGTCCTGCTCACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATg > 1:693517/1‑151 (MQ=255) ccatcaccGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCgc > 1:1114932/1‑151 (MQ=40) catcaccGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCgc > 1:420344/1‑150 (MQ=40) | TGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCC‑GGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCTGCACACGCAGGGTGAAGGCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCGGTGGGCAGGCTCAGGGCACAGACGGATGCCGC > NC_000913/730599‑730895 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |