Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,896,628 | 0 | A | C | 54.5% | 19.7 / 28.9 | 22 | V145G (GTC→GGC) | bglB | 6‑phospho‑beta‑glucosidase BglB |
Reads supporting (aligned to +/- strand): ref base A (8/2); new base C (0/12); total (8/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.41e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGAT > NZ_CP009273/3896543‑3896710 | gATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGc < 2:189009/90‑1 (MQ=255) ttAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTa < 2:58938/90‑1 (MQ=255) 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) | GATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTTGAT > NZ_CP009273/3896543‑3896710 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |