Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 2,088,271 | 0 | A | C | 60.0% | 2.3 / 14.7 | 15 | intergenic (‑265/+52) | htpX/alaA | Protease HtpX/Glutamate‑pyruvate aminotransferase AlaA |
Reads supporting (aligned to +/- strand): ref base A (0/6); new base C (9/0); total (9/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.00e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTG‑CGC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCG‑AATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAGGTCACCACGCGG > ACB143/2088141‑2088407 | gggATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTG‑CGC‑Ccccccaccc > 5:271356/1‑136 (MQ=255) gggATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTG‑CGC‑Cccacccccc > 3:89416/1‑136 (MQ=255) ggATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTG‑CGc‑ccccccccccc > 1:121252/1‑135 (MQ=255) tAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTGGATCTG‑CGC‑CCCCCCCCCAcccc > 8:67297/1‑137 (MQ=255) gTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTG‑CGC‑CCCCCCCCCCCCCCcgccgg > 5:66491/1‑134 (MQ=255) gTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTG‑CGC‑CCCCCCCCCAccccccccgc > 3:52675/1‑132 (MQ=255) tGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTG‑CGC‑CCCCCCCCCCCCCCcccggggg > 5:69155/1‑133 (MQ=38) cAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTG‑CGC‑CCCCCCCCCCCCCCcccggccgac > 5:113431/1‑131 (MQ=38) aGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTG‑CGC‑CCCCCCCCCCCCCCCGCGTGGGAcc > 4:215309/1‑137 (MQ=38) aGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCCGCCGCCCCGCGGGCGGCGCTCGTTCTGCCGCGCCCCCCCCCCCCCCcgcccgcgc > 1:127744/1‑131 (MQ=255) ggcacggCCGCTGGGTGTTTGCCGCACTGGGTGTGTGGCCTGTGACATCCCGCGCCGCGCGGGCGGCGCTCTATCTG‑CGC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCt < 4:342831/132‑1 (MQ=255) cgtgggggGTTTGCCGCATGTGGGGTGGTGCGGGTGAGAGCTCGCGCCCGGCGGGCGGCGCTCTATCTG‑CGC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCG‑aa < 7:343738/133‑1 (MQ=255) ggTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTG‑CGC‑CCCCCCCCCCCCCCCGCGGGGAACCCCGCGGAAGCCCCCGCGCCCCCGGCGCACCCCACAGCGCAATAAc > 8:343738/1‑137 (MQ=21) gcgcCGCGCGGGCGGCGCTCGAGCTG‑CGC‑CCCCCCACCACCCCCCCGGGGGACCTCTCGGCGCCCCCCGCGCGCCCGGCGACCCCGGATGCG‑CATAGCACGTCCGGATAATGCCGAGACGGCCAAATTCCGcgtacc > 4:339481/1‑132 (MQ=255) aTCTG‑CGC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCG‑AATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAgg < 5:166578/137‑1 (MQ=255) tCTG‑CGC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCG‑AATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAGGt < 4:52675/137‑1 (MQ=255) cgc‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCG‑AATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAGGTcacc < 7:164617/137‑1 (MQ=255) cgc‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGAG‑AATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAGGTcacc < 7:334864/137‑1 (MQ=255) caccaccacCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCG‑AATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAGGTCACCAcgcg < 5:96194/137‑1 (MQ=255) accaccaccCCCGGGGCGAACCTCTCGAAGCCCCCGGCGCCCCCGTCGAACCCTACTGCG‑AATAACCTTTCAGGAAATTGCCGATCAGGCCAATTGCCGCTTCCAGGTCATCCACCAGCGGCAAGGGCACCACGCgg < 2:195703/137‑1 (MQ=255) | GGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTG‑CGC‑CCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCG‑AATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCACCCGCGGCAAGGTCACCACGCGG > ACB143/2088141‑2088407 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |