Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB122 | 2,308,321 | 0 | C | G | 35.4% | 68.5 / 15.3 | 48 | A230P (GCC→CCC) | fucA | L‑fuculose phosphate aldolase |
Reads supporting (aligned to +/- strand): ref base C (5/26); new base G (17/0); total (22/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.20e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCA > ACB122/2308186‑2308449 | cggCACCTGTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGcc < 2:286429/136‑1 (MQ=255) aGGCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGcc < 5:374746/137‑1 (MQ=255) ccTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑ttt < 6:181027/137‑1 (MQ=255) ccTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑ttt < 6:1862/137‑1 (MQ=255) tttGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGc < 6:205680/137‑1 (MQ=255) tttGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGc < 5:346641/137‑1 (MQ=255) ccGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCGC‑TTTGCCCac > 6:58267/1‑136 (MQ=255) gttgttGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTg < 4:48737/137‑1 (MQ=255) gAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCgg < 6:244301/137‑1 (MQ=255) aCAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGG‑GTTGACCCGGTCTTCTGGGATGgggggg > 4:27198/1‑137 (MQ=255) aCAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGC‑GTTGTCCAGTGCTTGGGGGATTgggggg > 5:231939/1‑137 (MQ=255) aCAGCTTCGATCAGCCCAGGCAGTCGCCATCGCCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTTGCCACTGCTACGGGGATTGGGCgg > 8:155929/1‑137 (MQ=255) gCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGC‑GTTTCCCAGTGCTTGGGGGGTGggggggttt > 3:81453/1‑134 (MQ=255) ttCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTc < 7:305337/137‑1 (MQ=255) cGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTC‑TTTGGCCCGTGCGTCGGGGGTGGGGGGGGTGTCAc > 5:62138/1‑137 (MQ=255) aTCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGC‑GTTGCGCAGTGGGTCGGGGATGGGGGGGGTGTGAccc > 5:56321/1‑136 (MQ=255) tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTC‑TTTGGCCAGTGCTTGGGGGTTGGGGGGGGTGTCACCgg > 3:285837/1‑137 (MQ=255) tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTC‑TTTGCCCCGGGGTTGGGGGTTGGGGGGGGTGTGACCgg > 5:67251/1‑137 (MQ=255) tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTC‑TTTGCCCAGTGCGTCGGGGTTGGGGGGGGTGTCACCgg > 6:189413/1‑137 (MQ=255) tCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGC‑TTTGGGCAGTGGGTGGGGGTTGGGGCGGGTGGGACCgg > 7:250206/1‑137 (MQ=255) aGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTC‑TTTGGCCAGTGGTTAGGGGTTTGGGCGGCTGTGACCGGCGCCCaa > 7:162624/1‑136 (MQ=255) cAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGc < 1:112151/137‑1 (MQ=255) cGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGTGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCcagc < 2:173500/137‑1 (MQ=255) ccATGTGCACGCAAGGCACTGCGCTCGTAGTATCTGAAGGCAGCACCGTGACGTTTGGGTGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCcagcag < 5:210/137‑1 (MQ=255) ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGGGC‑TTTGGCCAGTGGGTGGGGGGTGGGGGGGATGTGACCGGCGCCCATGGCccgcag > 7:104330/1‑137 (MQ=255) ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGC‑TTTGCCCAGTGGGTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCcaccag > 2:229569/1‑137 (MQ=255) ccATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCcagcag < 2:159471/137‑1 (MQ=255) cATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGGGC‑TGTGGCCAGTGGGTTGGGGGTTGGGGGGCTGTGGACGGCGCCCATGGCcacccgt > 3:331259/1‑132 (MQ=255) cATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGC‑GTTGGCCCGTGGGTCGGGGCTGGGGGGGGTGTGACCGGCGCCCATGGGCCGCAGt > 5:219077/1‑137 (MQ=255) cATGGGCACCCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGC‑CATTACCACGGCGACCGGGAGGCGGCGGGTATCACCCGCGCCCCTGGCccgcaga > 6:13048/1‑136 (MQ=255) aTGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCTG‑TTTGGCCCGTTCTTAGGGGATAGGGGGGTTGTAACCCGCGCCCATGGCcagaaatt > 6:210/1‑132 (MQ=255) cacggCACCGCGCGCGTACTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGTGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTTCATCTTtg < 8:157377/134‑1 (MQ=255) cACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTgcg < 2:265425/137‑1 (MQ=255) gcgcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGGCGCGTTTGCCC‑GTGGGGGGGGGGTGGGGGGGCTGTGACCGGCGCCCCAGGCCCGCCGGTGGCTCTTTTGTgcgcgcggg > 1:107090/1‑134 (MQ=255) gcgTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGaa < 3:27198/137‑1 (MQ=255) agCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGGCGTGTGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCaga < 3:114000/137‑1 (MQ=255) aGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAgagcgagc < 5:189413/137‑1 (MQ=255) cagcaCCGTGACGTTTGGGGGTGGAAATCCAGTCGGGGTCCTC‑TTGGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACg < 8:104771/137‑1 (MQ=255) cGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTTACCAGTGGTTCGGGGATGGGGCGGATGTCACCGGCGCCCAAGGCCAGCAGTTGCATCTTTTCTTCGCGCTCGAACAGCAGAGCGCGCACCCAGGCGTACtc > 6:166854/1‑137 (MQ=255) cGTTTGGGGGTGGAAATCCAGTCGTGGGCCGC‑TTTTGCCAGTGCGTCGGGGGTTGGGGGGATTTGACCCGCGCCCAAGGCCAGCAGTTGCATATTTTCTTCGCGCGCGAACACCAGACCGAGCCCGCCCGCTTccac > 4:191691/1‑135 (MQ=255) ttGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCtcgat < 6:275619/137‑1 (MQ=255) ggTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAAc < 4:199816/137‑1 (MQ=255) ggTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAAc < 5:162778/137‑1 (MQ=255) ggTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAAc < 1:159477/137‑1 (MQ=255) ggTGGAAATCCAGTCGTGGGCCGC‑TTTTGCCAGTGGTTAGGGGATTGGGGGGATGTCACCCGCGCCCCAGGCCAGCAGTTGGAACTTTTATGCGCGCGCGCAAACCAGAGCGAGCACCCAGGCTTACTCGATAGaaa > 4:170965/1‑136 (MQ=255) taaaTCCAGGCGTGGGCCTC‑TTTGGCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTGCTCGATCGAACGccc < 7:81073/136‑1 (MQ=255) aaaTGCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCa < 1:106156/137‑1 (MQ=255) cAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGAt < 3:191691/137‑1 (MQ=255) aGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCATCAGTTGGATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCAGAGCGATc < 5:352530/137‑1 (MQ=255) gTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATca < 5:58267/137‑1 (MQ=255) | AGGCACCTTTGCCGTTGTTGAGGAAAGAGTACAGCTTCGATCAGCCCAGGCAGTCGCCATGGGCACGCAAGGCACGGCGCGCGTAGTAGCTGAAGGCAGCACCGTGACGTTTGGGGGTGGAAATCCAGTCGTGGGCCTC‑TTTGCCCAGTGCTTCGGGGATGGGGCGGATGTCACCGGCGCCCATGGCCAGCAGTTGCATCTTTGCTGCGCGCTCGAACAGCAGAGCGAGCACGCAGGCTTCCTCGATCGAACGCCCAGCGATCA > ACB122/2308186‑2308449 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |