Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,067,745 | 0 | A | G | 47.2% | 33.8 / 10.5 | 36 | L511P (CTG→CCG) | yihQ | putative sulpholipid alpha‑glucosidase; alpha‑glucosyl fluoride glucosidase |
Reads supporting (aligned to +/- strand): ref base A (10/9); new base G (0/17); total (10/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.40e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTT > NC_000913/4067600‑4067890 | ggCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGggt < 2:25225/150‑1 (MQ=255) aaaTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGgtgt < 2:997610/151‑1 (MQ=255) aaaTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGgtgt < 2:202415/151‑1 (MQ=255) aaaTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGgtgt < 1:425301/151‑1 (MQ=255) aaaTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTC‑TTGCTGCGCTTCATCTCAAACAGGGTGgtgt < 1:100975/150‑1 (MQ=255) cGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGc > 2:776086/1‑150 (MQ=255) cTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCg > 2:978421/1‑149 (MQ=255) ccgaaaaccTGCTAGTTGTCGCCAGGCGGTTTCCCTTGGGGGGGCCGCATCATGGGCGGGAG‑GGCGCTAAAATCG‑CCCCA‑GCGCAGCAGCAGCTCTTTGCGCCCTTTCATCTAAACCAGGGTGGTGTACCCGCCAATGTCGCTGTGGTGc < 2:801576/142‑1 (MQ=255) cgcaaaccgccgtgTTGGCGCCGGGACGTTTCCCTTCGGGGGTGCCCCTCCGCGGCGAGAA‑GGCGCTAAAATCG‑CCCCG‑GCGCAGCCGCACCTTTTGGCTGGGTTTCATCAAAAACGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAg < 2:877497/138‑1 (MQ=255) aaaccTGCCTGTTGGCGCCAGGACGGTTCCCTTGGTGGGTGCGCATCATGGGCGTGAA‑GGCGCTGAATTCG‑CACCA‑GCGCACCAGCCGCTCTTTGCTCCGTTTCCCCCCAAAAGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGcc < 1:851147/147‑1 (MQ=255) aaaaCGGCCGGTTGCCGCCGGGACGGTTCCTTGGGGGGGGGCGCCTCAGCGGGGGAAG‑GGCGCTAAAATCC‑CACCA‑GCGCAGCAGCAGCTTTTTGCTGCGCTTCATCTAAAACGGGGTGGTGTACCCGCCAATGTCGCTGTGGTGCAGGcc < 2:396370/150‑1 (MQ=255) tgccgtttGCCGCCGGGCCGGTTCCCTCGGGGGGTGCGCATCATCGGGGTAAG‑GGCGCGAAAACCG‑CACCG‑GCGCAGCAGCAGCTTTTTGCTGCGCTTCATCCCAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc < 2:81527/145‑1 (MQ=255) cgcagttggccGCCGGGACGGTTACCTTCGTGGGTCCGCATCTTCGGCGGAAG‑GGCGCTAAAACCG‑CACCG‑GCGCAGCAGCAGTTTTTGGCGGCGTTTCATCTAAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc < 1:369837/141‑1 (MQ=255) gccgttggccGCAAGGACGGTTCCCTGGGGGGGTGCGCATCATCGGGGTGAG‑GGCGCTAAAGCCG‑CCCCA‑GCGGAGCAGCGGTTTTTTCCTGCGCTTACTCTCAAACGGGGGGGTATAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc < 1:210698/141‑1 (MQ=255) gCCGGTTGTCGCGGGGAGGGTTCCCTGGGTGGGTGCGCACCAGCGGCGTGAG‑GGCGCTAAAATCG‑CCCGG‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACGGGGGGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc < 1:408215/150‑1 (MQ=255) gCCATTTGTCGCAAGGGCGGTTACCTTGGGGGGGGGGCATCAGCGG‑GGGGA‑GGCCGTGAAAAGGACAACA‑GCGCAGCAGCACCTCTTGGCTCCGTTTCACCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCAGGTc < 2:913329/150‑1 (MQ=255) cccgttgggcGCCGGGACGTTTCCCTTCGGGGGTGCGCATCACCGGGGTAAG‑GGCGTTGAAACCG‑CACCA‑GGGCGGCAGCGGCTTTTTGCGGGGCTCCATCCAAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc < 1:648686/141‑1 (MQ=255) ttGTCGCCGGGGCGGTTCCCTGCGGGGGTGCGCATCCTCGGCGGGAA‑GGCGCTAAAATCG‑CCCCG‑GCGCAGCAGCAGCTTTTGGCTGCGCTACATCCAAAACGGGGTGGGTAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTc < 1:547182/150‑1 (MQ=255) cccAGGCCGTTTCCCTTGGGGGGTGCGCATCTTCGGCGTGAA‑GGCGCTGAAACCG‑CACCG‑GCGCAGCAGCAGCTTTTTGCTGCGTTTCATCAAAAACGGGGGGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGcc < 1:474115/150‑1 (MQ=255) ccAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCa > 1:1033525/1‑151 (MQ=255) cgtttccctcggTGGGTGCGCATCACGGGCGGGAA‑GGCGCGAAAACCG‑CACCG‑GCGCAGCAGCACTTTTTTGCTGCGCTTCATCTCAAAGGGGGGGGTGAACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGAc < 1:439280/139‑1 (MQ=255) gCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGggt < 2:118995/67‑1 (MQ=255) gCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGggt > 1:118977/1‑67 (MQ=255) ggTGAAGGGCGCTAAAAT‑G‑CACCAGGCCCAGCGGCCGCTTTTTGCTGCGTTTCATCTAAAACGGGGTGGTGAACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGc < 1:766105/149‑1 (MQ=255) cGTGAG‑GGCGCTGAAATCG‑CACGG‑GCGCGGCAGCGGCTTTTTGCTGCGCTTCATCCAAAACGGGGTGGTGTACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGcc < 1:195716/150‑1 (MQ=255) cgttaaaaCCG‑CCCCG‑GCGCACCAGCCGCTTTTTGCGGCGCTTCTCCCAAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCgag < 2:174940/146‑1 (MQ=255) aaaaccgcccaggcgcAGAGGCGGCTTTTTGCTGCTTTTCATCTAAAAGGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACt < 2:796662/138‑1 (MQ=255) cggagccgcacttttttGCTGCCCTTCCCCCAAAACGGGGGGGGTTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTc < 2:26559/136‑1 (MQ=255) cagcagcTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcg > 2:859950/1‑151 (MQ=255) gcTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcccg > 1:857304/1‑150 (MQ=255) cgcTTCATCACAAACGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCAtgg < 2:158992/151‑1 (MQ=255) ctttaaTCCCAAAGGGGGTGGTGTACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCAtggt < 2:256136/146‑1 (MQ=255) aTCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGt > 2:29895/1‑151 (MQ=255) aTCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGt > 1:778701/1‑151 (MQ=255) aTCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCCGGTCGCCCGCCCACATCATGGGGGAGt > 2:1052554/1‑151 (MQ=255) cAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGCCGAGACTCCAGTCGACGTTCTGGGCCCCCGCCCACATCATGGCGGAGTAttt > 2:510083/1‑151 (MQ=255) | GGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAA‑GGCGCTGAAATCG‑CACCA‑GCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTT > NC_000913/4067600‑4067890 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |