Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,896,640 | 0 | G | C | 54.2% | 18.7 / 35.8 | 24 | A141G (GCT→GGT) | bglB | 6‑phospho‑beta‑glucosidase BglB |
Reads supporting (aligned to +/- strand): ref base G (10/1); new base C (0/13); total (10/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.61e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGATCCCCGCCTGCGCCAT > NZ_CP009273/3896558‑3896725 | cAGCCATAACGCCATTTTTTGTTGGTAGCGAGTAAAGCCCGTGGGGGCGTAATGCTCGAAGTGGTCGATGCCCGCTCGATTACCCCAAcc < 1:394338/90‑1 (MQ=255) ttttTGTTGGGGGCGAGAAAAGCCCGGGGGGGGCTAACGATCGAAGTAATCGATTCCCGCTAAATTACCCAACCCGCCGTAGTTTTTCAc < 2:318649/90‑1 (MQ=255) ttttTGTGGGTAGGGAGAAAAGCCCGGGGGGGCGTAATGCTCGAAGGGGCCGATGCCCGCTCGATTACCCCACCCGCCGTAGTTTTTCAc < 2:323065/90‑1 (MQ=255) tttATTTGGGTGGGGAGAAAAGCCCGGGGGGGCGTAATGCTGGAAGTGATCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCAc < 2:134900/90‑1 (MQ=255) gTTGGTAGGGAAAAAAGCCCGGGGGGGCGTAATGCTGGAAGTGGTCGATCCCCGCTCGATTACCCCACCCGCCGTAGTTTTTCACCAGcc < 2:285115/90‑1 (MQ=255) gTTGGTAGCGAGAAAAGCCCGGGGGGGGGTAATGCTGGAAGTGATCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGcc < 2:269866/90‑1 (MQ=255) gTTGGTAGCGAGAAAAGCCCGGGGGGGCGTAATGCTCGAGGGGGTCGATGCCCGCTCGATTACCCAAACCGCCGTAGTTTTTCACCAGcc < 2:437893/90‑1 (MQ=255) gTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCata > 1:358675/1‑90 (MQ=255) aaGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCATCAGCCCATATGGCATTTCg > 2:380162/1‑90 (MQ=255) gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAGTgg > 1:89204/1‑90 (MQ=255) gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg > 2:226916/1‑90 (MQ=255) gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg > 2:424826/1‑90 (MQ=255) gggggCGTAATGCTGGAGGGGGTCGATGCCCGCTCGTTTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGa < 2:149330/88‑1 (MQ=255) gCGGGCGTAATGCTCGAGGTGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGa < 1:142547/90‑1 (MQ=255) ggggCATAATGCTCGAGGGGGTCGATCCCCCCTCGATTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGAt < 2:444577/89‑1 (MQ=255) gTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGt > 1:289170/1‑90 (MQ=255) tCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGc > 1:180930/1‑90 (MQ=255) ggaggTGGTCGATGCCCGCTATTTTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCg < 2:259379/86‑1 (MQ=255) ggggtCGATGCCCGCTCTATTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCtt < 2:423998/86‑1 (MQ=255) gtCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGAt > 2:78853/2‑90 (MQ=255) gCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGATCCCCGCCTgc > 1:62384/1‑90 (MQ=255) tttACCCCACCCGCCGTAGTTTTTCCCCAGCCCATATGGCATTTCGTAATGGGATAACGTTAACAGCGGCTTGATCCCCGCCTGCGCCAt < 2:342886/89‑1 (MQ=255) aGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCtt < 1:105519/69‑1 (MQ=255) aGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCtt > 2:105519/1‑69 (MQ=255) | CAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGATCCCCGCCTGCGCCAT > NZ_CP009273/3896558‑3896725 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |