Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,061,105 | 0 | A | G | 52.0% | 4.7 / 25.9 | 25 | L511P (CTG→CCG) | yihQ | sulfoquinovosidase |
Reads supporting (aligned to +/- strand): ref base A (5/7); new base G (0/13); total (5/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.49e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTT‑CGTGGGTGCGCATCATCGGCGTGAAG‑GCGCTGAAAT‑CGCACCAGCGCAGCAGCAGCTCTTTGCTG‑CGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACAT > NZ_CP009273/4060970‑4061235 | gCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTT‑CGTGGGTGCGCATCATCGGCGTGAAG‑GCGCTGAAAT‑CGCACCAGCGCAGCAGCAGCTCTTTGCTG‑CGCTTCATCTCAAACAggg < 2:922251/139‑1 (MQ=255) gTCAAACTGCCGTTTGTCGCCAGGACGGTTACCTT‑CGGGGGTGCGCATCATCGGCGTGAAG‑GCGCTGAAAT‑CGCACCAGCGCAGCAGCAGCTTTTTGCTG‑CGTTTCATCTAAAACAGGGTGGTTTAACCGCCAATGTCg < 1:131870/139‑1 (MQ=255) aaccTGCCTGTGGTCGCCGGGACGTTTCCCTG‑CGTGGGGGCGCATCATCGGCGGGAAG‑GCGCAGAAAT‑CGCCCCAGCGCAGCAGCAGCTTTTTGCTG‑CGCTTCATCACAAAGAGGGTGGTGTAACCGCCAATGTCGCtg < 2:890797/136‑1 (MQ=255) aaccTGCCGGTTGCCGCCGGGCCGTTTCCCTC‑CGGGGGTGCGCATCATCGGCGGGAAG‑GCGCAGAAAT‑CGCCCCGGCGCAGCAGCAGTTCTTTGCTG‑CGCTTAATCCAAAACGGGGTGGTGTAACCGCCAATGTCGCtg < 2:646932/136‑1 (MQ=255) aaccTGCCGGTTGCCGCCAGGCCGTTTCCCTT‑CGGGGGTGCGCATCATCGGCGTGAAG‑GCGCAAAAAT‑CGCACCAGCGCAGCAGCAGCTTTTTGCTG‑CGCTTCATCTAAAAGGGGGTGGTGTAACCGCCAATGTCGCtg < 1:769289/136‑1 (MQ=255) gCCAGTTGTCGCC‑GGGGGGTTTCCTTGCGTGGGGCCGCATCATCGGGGGGAAG‑GCGCTAAAAT‑CGCACCGGGGCAGCAGCAGTTTTTTGCGG‑GGCTTCACCCAAAACGGGGTGTTGTAACCGCCAATGTCGCTgtggtg < 2:585349/139‑1 (MQ=255) gCCAGTTGCCGCCGGGACGTTTCCTTT‑CGTGGGTGCGCATCATCGGCGTGAAG‑GCGCTGAAAT‑CGCCCCAGCGCAGCAGCAGCTTTTGGCTG‑CGCTTCATCTCAAACGGGGGGGTGAAACCGCCAATGTCGCTgtggtg < 2:384393/139‑1 (MQ=255) ccAGTTGTCTCCGGGACGGTTCCCTT‑CGGGGGGGCGAATCTTCGGCGGGAAG‑TCGATGAATA‑CGAGCCACCGCAGCAGCTGTTTTTTGCGC‑CGTCTCTTCCCAAGCAGGGGGTGGTAACCGCCAATGTCGCTgtggtg < 2:850989/138‑1 (MQ=255) ccAGTTGTCCCCAGAGCGGTTCCCTG‑CGTGGGTGCGCATCATCGGTGTGAGG‑GCGCTGAAAG‑CGCACCAGCGCAGCAGCAGCTCTTGGCTG‑CGCTTCATCTAAAACGGGGGGGTGGAACCGCCAATGTCGCTgtggtg < 1:773045/138‑1 (MQ=255) ttGTCGCCAGGACGGTTACCTC‑CGTGGGTGCCCATCCTCGGCGTGAAG‑CCGCTAAAAC‑CCCCCCGGCGCAGCAGCGCCTTTTTGCTG‑CTCTTCACCTCAACCGGGTTGGTGAAACCGCCAATGTCGCTGTGGTGCAGGc < 1:490843/139‑1 (MQ=255) gccGCCAGGCCGTTACCCTG‑GGTGGGGGCCCACTACCGGAGTGAAG‑GCCTTAAAAC‑CGCGCCGGCGCAGCGGCGGCT‑TTTGGTGCCGCTTCACCACAAACAGGTGGGTGTAACCGCCAATGTCGCTGTGGTGCAGGcc < 2:488391/136‑1 (MQ=255) ttACCTT‑CGTGGGTGCGCATCATCGGCGTGAAG‑GCGCTGAAAT‑CGCACCAGCGCAGCAGCAGCTCTTTGCTG‑CGCTTCATCTCAAACa > 2:694905/1‑88 (MQ=255) ttACCTT‑CGTGGGTGCGCATCATCGGCGTGAAG‑GCGCTGAAAT‑CGCACCAGCGCAGCAGCAGCTCTTTGCTG‑CGCTTCATCTCAAACa < 1:694905/88‑1 (MQ=255) gggTGCGCATCAGC‑GCCTGAAGCGCGC‑GAAATCCCCCCCCGCGCAGCAGCAGCTTTTTGCTG‑CGCTTCATCACAAACAGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGa < 2:919393/136‑1 (MQ=255) gcATCTCCGGCGGGAAG‑GCGCTAAAAC‑CCCCCCAGCGAAGCAGCAGCTCTTTGCTG‑CGTTTCATCCCAACCGGGGTGTTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGcgccgc < 2:97722/139‑1 (MQ=255) atcatcGGCGTGAAG‑GCGCTGAAAT‑CGCACCAGCGCAGCAGCAGCTCTTTGCTG‑CGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGccgccg > 2:368812/1‑139 (MQ=255) agcGGCGTGAAG‑GCGTTAAAAC‑CCCCCCAGCGCAGCAGCAGCTTTTTGCTG‑CGCTTCATCTCAAACAGGGTGGTATACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGa < 1:129875/137‑1 (MQ=255) acgggggTGAAG‑GCTATAAAAC‑CGCCCCGGCGCAGCAGCGGTTTTTTCCTG‑CGCTTCATCAAAAACGGGGGGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGa < 2:223362/133‑1 (MQ=255) tGAAG‑GCGCTGAAAT‑CGCACCAGCGCAGCAGCAGCTCTTTGCTG‑CGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGa > 2:946288/1‑139 (MQ=255) cTGAAAT‑CGCACCAGCGCAGCAGCAGCTCTTTGCTG‑CGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGccctgt > 1:223362/1‑135 (MQ=255) aaaT‑CGCACCAGCGCAGCAGCAGCTCTTTGCTG‑CGCTTCATCTCAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGcc < 1:190510/137‑1 (MQ=255) cgcgGCGGCAGCTTTTGGCTG‑CGCTTCATCTCAAACGGGGTGATATAGCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTcc < 1:336829/139‑1 (MQ=255) tttGCTG‑CGTTTCATCTCAAACGGGGGGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTg < 1:292319/138‑1 (MQ=255) cgcTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcccg > 2:149828/1‑139 (MQ=255) cTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcccgcc > 1:644304/1‑139 (MQ=255) cTTCATCCAAAGCGGGGGGGTTAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcccgcc < 2:617853/139‑1 (MQ=255) tcacAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCAcat < 1:494479/136‑1 (MQ=255) | GCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTT‑CGTGGGTGCGCATCATCGGCGTGAAG‑GCGCTGAAAT‑CGCACCAGCGCAGCAGCAGCTCTTTGCTG‑CGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACAT > NZ_CP009273/4060970‑4061235 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |