| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 993,913 | 0 | G | C | 52.9% | 9.6 / 19.3 | 17 | G44G (GGC→GGG) | ssuB | aliphatic sulfonate ABC transporter ATPase |
| Reads supporting (aligned to +/- strand): ref base G (4/4); new base C (0/9); total (4/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.94e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-03 | |||||||||||
CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCT > NC_000913/993763‑994054 | ccACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCgg < 2:114740/151‑1 (MQ=255) aCGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGcc < 1:424458/151‑1 (MQ=255) cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACaaa > 1:246672/1‑151 (MQ=255) tcatTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg > 2:123049/1‑151 (MQ=255) ctgcggttttttCCTGGATTTCCACCCAGGGGTGGGGGCCCCGCAAACACTCCGCCGGGGGTGGGCGTTCCC‑AGCCCTGCCGGCAGGCGCAGCGGGTTACTTTCCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGaa < 1:376428/141‑1 (MQ=255) cttgtCTCCCTAAATTT‑CACCAAGCGGTGGGGGGCCCGTTAACACACCG‑CTGGGGGTGGCCTTTCCAAGCCCCCCCAGCAGGCCAAGCAGTGTACTTTTGCCACCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATCTGt < 2:430680/149‑1 (MQ=255) ttCCAGATTTC‑CACCCAGGGGGGGGGCCCCCGCTAACACCCCGCGTGCGGTTGGCTTTCCC‑GGCCCTGCCAGCAGGCGCAGCAGGTTATTTTGCCCCCCCCCGCTGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc < 1:189072/151‑1 (MQ=255) gtggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt > 1:217209/1‑151 (MQ=255) ggggCCCGCTAACCCACCGCCGGGGGTTGGCGTTCCC‑AGCCCGGCCAGCAGGCGCAGCGGGGTTTTTTTCCCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc < 2:398517/147‑1 (MQ=255) cGGTGGGCTTTCCC‑AGCCCTCCCAGCGGGCGCAGCAGGGTTCTTTCCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt < 1:454822/150‑1 (MQ=255) tttccc‑AGCCCTCCCAGCAGGCGCAGCAGGGT‑CTTTTCCCACCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa < 1:400968/147‑1 (MQ=255) ccggccccgCCGGCAGCCGCAGCGGGTTCCTTTCCCCACCCCGGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg < 1:105026/143‑1 (MQ=255) tttctttcccccccaccGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt < 1:115611/139‑1 (MQ=255) ttcccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc < 2:426277/142‑1 (MQ=255) ttGCCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc < 2:111982/151‑1 (MQ=255) tcccccccccGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc < 2:302666/143‑1 (MQ=255) caccGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc < 1:290000/151‑1 (MQ=255) ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct > 2:429606/1‑150 (MQ=255) | CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCT > NC_000913/993763‑994054 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |