Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,598,904 | 0 | T | C | 22.9% | 57.6 / 11.2 | 35 | S436G (AGT→GGT) | lsrK | autoinducer‑2 (AI‑2) kinase |
Reads supporting (aligned to +/- strand): ref base T (15/12); new base C (4/4); total (19/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCAGGCGTTCTCCGGTTTCTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTA > NC_000913/1598766‑1599046 | ccAGGCGTTCTCCGGTTTCTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTCATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCACGAAAGATCAGCGAGAAATTGCCTCCAt > 1:114097/1‑143 (MQ=255) cAGGCGTTCTCCGGTTTCTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACCCCATa > 1:618293/1‑143 (MQ=255) cAGGCGTTCTCCGGTTTCTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAAGGGGACATCCTAATGCAGTCGCTTCTTTGACCGCCGGACTGTTGCCGGGTAATCCCGAGACATCAGCGAGACTTCGACTCCATa > 1:4175/1‑143 (MQ=255) tgttgtgcGGTCTCTGCCAGTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAAttt < 2:340333/136‑1 (MQ=255) gTTTCTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGa > 1:500530/1‑143 (MQ=255) cTGCCATTGATGAAAAACTTCCGGCACCGACGCCAGCGGCAATGGCACATCCTGATGCAGTGGCTTCTTGGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCt < 2:172210/143‑1 (MQ=255) gatgaAAAAATTCCGGCACCGACGCCATCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACCCCATAATTTCCCTTTTGAACCTCCGCCTGc < 1:217590/143‑1 (MQ=255) gaAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACCCCATAATTTCCCTTTTGAACCTCCGCCTGCaaa > 1:163661/1‑143 (MQ=255) aaaaaaTTCCGGCACCGACGCCAGCTGCAATGGCAGATCCTAATGCAGTGGCTGCTTTGAGCACCGGAAAATTGACGGGTAATCCCGAGACATAAGAGAGAATTTGCCTCCATAATTTCCCTTTTGAAACTCCGCCCGAAAAg > 2:43248/1‑143 (MQ=255) cACGGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGaa < 1:430927/143‑1 (MQ=255) gACGCCAGCTGCAATCGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGAt < 1:152396/143‑1 (MQ=255) aCGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATg > 1:663731/1‑143 (MQ=255) ccccAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACCCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGa < 1:657693/141‑1 (MQ=255) ccAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGTGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCGCCGCCTGCAAAGACTAACGAGGAAGGATGAat < 1:669126/143‑1 (MQ=255) ccAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAat < 2:473669/143‑1 (MQ=255) gCACTGGCACATCCTAATGCAGTGGCTTCTTTGACCACTGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCga < 2:711940/143‑1 (MQ=255) gCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCga < 2:569965/143‑1 (MQ=255) gCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCga < 1:316611/143‑1 (MQ=255) cAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGGCTAACGATGAAGGATGAATATTCgag < 1:234525/143‑1 (MQ=255) aaTGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATCATTTCCCTTTTGAACCTCCGCCTGCACAGTCGAACGCTGACGGCTGAATGTTCgcga > 1:695112/1‑143 (MQ=255) aaTGGCACATCCTAATGCAGTGGCGTCTCTGCCCAATGCACAAGTGACGTTTAATCCCGAGACAGCAGCGAGAATTTGACCCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCgaga < 2:760907/143‑1 (MQ=255) cacaTCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGACTTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAACGCCTAAGGATGCACGATGAATCTCCGACAAATCa > 2:388909/1‑143 (MQ=255) cAGTGGCTTCTTTGACCACCGGCATATTGACGGGTAATCCCGCGACATCAGCGAGCGTTTGACCCCATAATTTACCTTTTGAACCTACGCCTGCAGAGACTAACGATGAAGGATGAATCTTCGCGAACGCAGAAATTggctgc > 2:508111/1‑143 (MQ=255) gTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAACATTCGAGAAATCAGCAATTTGCTGCaa < 1:478125/143‑1 (MQ=255) cttcttTGACCACCGGAATATTGACGGGTAATCCCGAGACAGCAGCGAGAATTTGACTCCATAATTTCCCTCTTGCACTTGCTCTTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTa < 2:725774/143‑1 (MQ=255) accGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACCCCATAATTTCCCTTTTGAAGCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCCAGAAATCAGCACTTTGCTGCAACCTACACGCTCATAc > 1:660672/1‑143 (MQ=255) gAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTACCCTAGTGAAGCTGAGACTGGTGACACCAACGATGAAGGATGAATATTCGCGAAATCAGCACTTTGCTGTACGTTACACGCTGATAGAATc > 2:277010/1‑143 (MQ=255) gggTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGATACACGCCGATACAATCACCGCAtttt > 2:488152/1‑143 (MQ=255) aTCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTcc > 1:425346/1‑143 (MQ=255) ccGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATGCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAgc > 1:213321/1‑143 (MQ=255) cATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACgg > 1:761433/1‑143 (MQ=255) cATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACgg > 1:628009/1‑143 (MQ=255) aGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAAACTAACGATGAAGGATCACTCTTCGAGAAATCAGCAATTCGCTGCAAGTTACACGCTCACACAATCGCAGCCTTCTCTTCCACCGAAccggacg > 2:731219/1‑137 (MQ=255) gaATTTGACCCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTc < 1:566032/143‑1 (MQ=255) aaTTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCg > 2:578840/1‑143 (MQ=255) tCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCACTTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTa > 1:194196/1‑143 (MQ=255) tCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTa > 1:609413/1‑143 (MQ=255) | CCAGGCGTTCTCCGGTTTCTGCCATTGATGAAAAAATTCCGGCACCGACGCCAGCTGCAATGGCACATCCTAATGCAGTGGCTTCTTTGACCACCGGAATATTGACGGGTAATCCCGAGACATCAGCGAGAATTTGACTCCATAATTTCCCTTTTGAACCTCCGCCTGCAAAGACTAACGATGAAGGATGAATATTCGAGAAATCAGCAATTTGCTGCAAGTTACACGCTGATACAATCGCCGCATTTTCTTCCAGCGCACGGAACAATGTCGCTTTGTTA > NC_000913/1598766‑1599046 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |