Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB122 | 3,063,521 | 0 | T | C | 32.3% | 38.1 / 12.6 | 40 | intergenic (+32/‑322) | adh_2/pqqD_2 | Aldehyde dehydrogenase/PqqA binding protein |
Reads supporting (aligned to +/- strand): ref base T (8/13); new base C (2/8); total (18/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.28e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGTTCCAGAACCGGTG‑AAGAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGT > ACB122/3063390‑3063654 | cGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCCCAccccc > 3:372939/1‑135 (MQ=255) ggTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCACCcacc > 1:6156/1‑136 (MQ=255) ccGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCAcc < 8:96263/137‑1 (MQ=255) ccGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCAcc < 6:213577/137‑1 (MQ=255) ccGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCAcc < 2:333775/137‑1 (MQ=255) gaCCCACAAGATGATGCTTGACCATTATCAGCAGACCAACAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGt < 5:330606/137‑1 (MQ=255) aCCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTgggt > 2:28318/1‑135 (MQ=255) acaAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGTTCCAg < 7:96124/137‑1 (MQ=255) caAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGTTCCAGa < 4:33650/137‑1 (MQ=255) caAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGTTCCAGa < 6:19000/137‑1 (MQ=255) aaGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTACCCGCTGGGTGCAgcc > 6:339027/1‑135 (MQ=18) gatgCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGTTCCAGAACCGGt < 4:111141/137‑1 (MQ=255) cTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCCGGC‑TTTCCCCCACCCCCA‑A‑GTTTACCCCGAGTTCCAGAACCCGGg‑aaa > 5:340508/1‑136 (MQ=18) cATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑C‑GTTAACCGCTGGTTCCAGAACCGGTG‑AAGAACCtg > 8:93030/1‑137 (MQ=255) tCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑GTGCCCCCCCCCCCC‑C‑GGTTACCCCGGGGTGCAGGACCCGTG‑GAGAAACCGTgggg > 3:50809/1‑137 (MQ=14) cagcagACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCGCCCACT‑G‑GTTCACCGCTGGGTCCAGAACCGGTG‑AAGAACCTGTggggg > 2:100887/1‑137 (MQ=255) aaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCAcc < 8:54069/94‑1 (MQ=255) aaGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCCCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCAcc > 7:54069/1‑94 (MQ=255) tgGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTTCCCCCGCCCACACG‑GTTAACC‑CGGGGTGGCCAACCCGGG‑TAGAACCGGTTGGGGCGGGCGCG‑CCCCCcaaaag > 1:5025/1‑135 (MQ=12) ggTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTTCCCCCACCCACA‑C‑GTTGAACGCAGGGTGCAGAAACCGGG‑AAGAACCAGTGGGGGCGGGCGCG‑CCCGCAAAAcgg > 5:149958/1‑137 (MQ=21) gTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCTTTGCCCCAACCCCCA‑‑‑GTTAAACCCTGGTTGCAGAAACCGTG‑AAAAAACCCTTGGGGCCGGCGCC‑CCCCCAAAACGGc > 7:30004/1‑137 (MQ=9) gTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGCATTGCCC‑CCCCCCCA‑GCGTTACCCCTTGGTTGATCAACCGG‑G‑GAGAACCGGTTGGGGGGGGCGCGCCCCCCAAAAggg > 5:365113/1‑137 (MQ=9) tGAGCTACCACATCAACCCGCTGGGCTTCTTCTAACCCGCGTCGCGGCCTTCGCGGGC‑TGGCCCGCTCCCACC‑G‑GTGCACCGCGGGTTCCAGAACGGGTG‑AAGAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGcc < 7:375922/137‑1 (MQ=255) tagatACCACATCAACCCGCGGGGCTCATTCAAAACCGCGGCGGGGCCTTCGCGGGC‑TTCCCCGCCCCCACC‑G‑TTCCCCCGCTGGTCCCACAACCGGTC‑AACAACCTGGGGGGGGGGGCGAG‑CCCGCAAAAGGGCCg < 4:222884/133‑1 (MQ=255) gAGCAACGCCAAAAAGCCGCTGGGCTTCTTATAACCCGCGGCGCGGCCTCCGCGGGC‑TGGCCCGCCCCCACC‑G‑GTTCACCGCGGGTTCCACAACCGGTC‑AATAACCTGGGGGGGCGGGCGAG‑CCCGCAAAAGGGCCg < 2:476/137‑1 (MQ=255) gCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTTCCCCCGCCCACA‑C‑GGTGACCGCGGGGTTCAGCACCCGGG‑AAAAAACAGTGGGGGCGGGGGCG‑CCCGCCAAAAGGCCCGc > 6:2131/1‑137 (MQ=16) atgACATCAATCCGCGGGGCTGCTACAAACCCGCGGCGCGGCCTTCGCGGGC‑TGGCCCGCCCCCACA‑G‑GTCCCCCGCTGGTCCCCCACCCGGTT‑AACAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGcaca < 1:231784/135‑1 (MQ=255) aaTCAATCCGCTGGGCTTCTTATAACCCGCGGCGCGGCCTCCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTCCACCGCTGGTACCACACCCGGGC‑AAGAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCAcatgca < 7:174196/136‑1 (MQ=25) atcacaccGCGGGGCTGCTTATAACCCGCGCCGCGGCCTCCGCGGGC‑TGGCCCGCCCCCCCC‑G‑GTGCACCGCGGGTCCCACACCCGGTC‑AACAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCAcatgcat < 4:211646/131‑1 (MQ=255) agaaccccGCGGGGCTTCTACAAACCCGGGCCGGGGCTTGGGCGGGC‑TGCCCCGCCCCCACC‑G‑GGCCCCGGCGGGTCCCACACCCGGTA‑AACAACGTGGGGGGGGGGGCGAG‑CCCGCAAAAGGGCCGGCAcatgcat < 4:50809/131‑1 (MQ=255) cccGGTGGGGTTATTATAACCCGCGGCGCGGCCTTCGCGGGC‑TGCCCCGCGCCCCCC‑G‑GTCCACCGCTGGTACAACACCCGGTT‑AATAACGTGGGGGGGGGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCt < 4:35321/136‑1 (MQ=255) ccGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTTCCCCCTCCCACA‑C‑GTTGACCGCTGGTTCCAGAACCGGGG‑AAGCACCAGTGGGGGCGGGCGCG‑CCCGCCAAATGGCCCGCCCATGCCTGGccga > 4:371003/1‑135 (MQ=255) gCTTCTTCTAAACGGCGGGGCGGCGTTCGCGGGC‑TTGCCCGGCCCCCCA‑G‑GTTCACCGCTGTTTCCACAACCGGTT‑AACAACGTGGGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTAc < 2:86166/137‑1 (MQ=255) cgcgGCCTTCGCGGGC‑TTGCCCCCTCCCACA‑G‑GGTGACCGCTGGGTGCAGAACCGGGG‑AAGAACCGGTGGGGGCGGGCGCG‑CCCGCAAAACGGCCCGCCCAAGGCTGGCCGACCCGGACCACCGCCCCCCCCCAGcg > 2:51891/1‑136 (MQ=255) cccTGCGGGGGC‑GTGCCCGCCCCCCCT‑G‑GTTCACCGGTGGTACCAGAACGGGTT‑AAGAACCGGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGtt < 5:293031/136‑1 (MQ=255) cgggggC‑TGGCCCGCTCCCACA‑G‑GTGCCCGGCTGTTTCCAGACCCGGTG‑AACAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCaaa < 7:238786/134‑1 (MQ=255) ccgGGC‑TGGCCCCCCCCCACC‑G‑GTCCCCCGCTGGTCCCA‑AACCGGGGAAACAACGGGGGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAg < 3:153244/136‑1 (MQ=255) ggttGCCGGCTCCCACA‑G‑GTTCACGGGGGGTTCCAGACCCGGTG‑AAGAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGc < 6:86256/135‑1 (MQ=255) gCCCCCTCCCACA‑C‑GTTGAACCCTGGGTGCAGAACCCGGG‑AAGAACCTGTGGGGGCGGGCGCG‑CCCGCAAAAAGGCCCGCCCAATCCTGGCCGACCCGGACCAACCCCCCCCCCCAGCGCGTTGCAAAGTAGCaata > 6:293031/1‑134 (MQ=255) cGCTCCCCCA‑G‑GTTGACCCCGGGGTCCAGAACCGGTG‑AAGAACCTGTGGGGGCGGGCGCG‑CCCGCAAAAGGGCCGGCCCAAGCATGGCCTACCCGGACCAACGCCCCCCGCAACCGCGTTCCAAAGTAGCATTCGGt > 8:206007/1‑137 (MQ=255) | CGGTGTTGGCCGTGAGACCCACAAGATGATGCTTGACCATTATCAGCAGACCAAGAACCTGCTGGTGAGCTACGACATCAATCCGCTGGGCTTCTTCTAAACCGCGTCGCGGCCTTCGCGGGC‑TTGCCCGCTCCCACA‑G‑GTTCACCGCTGGTTCCAGAACCGGTG‑AAGAACCTGTGGGGGCGGGCGAG‑CCCGCAAAAGGGCCGGCACATGCATGGCCTACCCGTACCAACGCACCCCGCAAGCTCGTTCCAAAGTAGCATTCGGT > ACB122/3063390‑3063654 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |