Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | ACB143 | 4,841,918 | C→G | 19.8% | intergenic (+7/‑99) | CJ019_04539 → / → hemN | hypothetical protein/Oxygen‑independent coproporphyrinogen IIIoxidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 4,841,918 | 0 | C | G | 19.8% | 49.9 / 3.8 | 48 | intergenic (+7/‑99) | CJ019_04539/hemN | hypothetical protein/Oxygen‑independent coproporphyrinogen IIIoxidase |
Reads supporting (aligned to +/- strand): ref base C (4/30); new base G (6/4); total (11/37) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.22e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 |
TGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCC‑ACACCAGCACTGGCTGATGGGC‑CA‑TTGACGGGGGC‑TGC‑AA‑A‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACC > ACB143/4841788‑4842043 | tgCCTAACGGATGACCGCCCTGTTGCGGCGACGCAGGGTGCGAGTGGCTGGGGGTGTGCTGGTTATGATATGCGTTATCTGGACCCTGCCGGGCCA‑ACACCACCACTGGCTGATGGGC‑AA‑TTGACGGGGGC‑TGC‑aa < 3:353342/136‑1 (MQ=255) gCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCC‑ACACCAGCACTGGCTGATGGGC‑CA‑TTGACGGGGGG‑Tggca > 2:104005/1‑132 (MQ=255) gCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCC‑ACACCAGCACTGGCTGATGGGC‑CA‑TTGACGGGGGG‑TGC‑aa > 2:298623/1‑135 (MQ=255) gCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCC‑ACACCAGCACTGGCTGATGGGC‑CA‑TTGACGGGGGC‑GGC‑ga‑a > 5:13520/1‑136 (MQ=255) gATGAACGCCGTGGTGCGGCGAGGCAGGGGGGGTGTGGCTGGGGGTGTGGTGGGGATGCTTTTGGGTAGGTGGCCCCTGCCGGGCCC‑ACACCAGCAGGGGCTGAGGGGC‑CA‑TTGAGGGGGGG‑TCC‑AA‑C‑GCAGcccc < 6:267767/136‑1 (MQ=255) aaCGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCC‑ACACCAGCACTGGCTGATGGGC‑CA‑TTGACGGGGGC‑GGC‑AA‑C‑CCAGcccccaa > 4:338528/1‑132 (MQ=255) gACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCC‑ACACCAGCACTGGCTGATGGGC‑CA‑TTGACGGGGGG‑GGG‑CA‑A‑GCACGCCCAAACCCCGTGATACaaaac > 5:344964/1‑133 (MQ=255) cGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCC‑ACACCAGCACTGGCTGATGGGC‑CA‑TTGACGGGGGG‑GGC‑GA‑A‑GCAGCCCCACTGCCCGGGGTAACAATCAACACAGt > 7:229977/1‑136 (MQ=255) gCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCC‑ACACCAGCACTGGCTGATGGGC‑CA‑TTGACGGGGGG‑GGG‑CG‑A‑GCAGCCCCAAAGCCCATGATAAAAAAAAAAACACttt > 3:34116/1‑135 (MQ=255) gAGTGGCCGGGGGTGTGGTGGTTTTGCTATTGGTTTTCTGGACCGTGCGGGGCCA‑ACACCCGAAGGGGGT‑TTGGGC‑CATTTGCGGGGGGC‑TCA‑AA‑A‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCt < 8:221087/136‑1 (MQ=255) ggggCTGGGGGTGGGGTGTTGTTGCTATGCGTTACCTGGCCCCGGCGGGCCAC‑ACACCAGCTGTGGGT‑TTGGGCACT‑TTGCCGGGGTA‑TCA‑AA‑A‑CCCCCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTAc < 6:229515/134‑1 (MQ=255) gTGGCTGGGGGGGGGGTGGTTATGCTTTTGGTTATCTGGCCCCGGCGGGCCCA‑ACACCCGCAGGGGGTGTGGGGC‑CT‑TTGCGGGGGGC‑GCC‑AA‑C‑GCACCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTAc < 6:269410/136‑1 (MQ=255) tGGCCGGGGGTGTGGTGTTTATGTTATTGGTTTTCGGCACCCTGCGGGCCCC‑ACACCAGCTGTGGTTTTTGGGC‑CA‑TTTAGGGGGGC‑TGC‑AA‑A‑GCCGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTAc < 8:238244/135‑1 (MQ=255) ggggTGTGCGGGTGATGATATTCGGTATGGGGCCCCTGCGGGGCAC‑ACACCAGCAGTGGCTTAGGGGC‑CT‑TTGAGGGGGGC‑TCA‑AA‑C‑CCTGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGc < 1:298623/134‑1 (MQ=255) ggggtggCTGGTGAGGCTATTCTTTATCT‑GCCCCCGCCGGGCCCAACACACCCAGTGCAGTATGGGC‑TA‑TTGCGGGGGGA‑TGC‑AA‑C‑GCCCCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGcc < 6:196765/130‑1 (MQ=255) ttgtgCGGTTTATGCTATGGGGTATCTGCACCCTCCCGGGCCC‑ACACCCGAAGTGGCTGAGGGCC‑CA‑TTGAGGGGGGA‑GGA‑AA‑A‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCAt < 5:284056/135‑1 (MQ=255) gtgtGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCC‑ACACCAGCACTGGCTGATGGGC‑CA‑TTGACGGGGGG‑GGG‑CA‑A‑GCACCCCCCAAGCCCGGGGTACAAAAAAACACACTTTGCGACAGACCGccc > 7:105757/1‑134 (MQ=255) gTGATCCTATTCGTTATCTGGACCCGGCCGGGCCC‑ACACCAGCAGTGGCTTATGGCC‑CA‑TTTCCGGGGGC‑TGA‑AA‑C‑GCCGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCg < 7:381898/136‑1 (MQ=255) aTTCGGTATCTGGACCCTGCCGGGCCC‑ACACCAGCACTGGCTGATGGGC‑CA‑TTGACGGGGGC‑GGG‑CA‑A‑GCAGCCCCCAAGCCCATGCTAAAAAAAAAAACAGTTTACTACAGACGGCCCAACACACCGGTGcccgg > 2:73434/1‑132 (MQ=255) tcgggccccGGCCGGCCCC‑ACACCCGCACGGGCGGATGGGC‑CT‑TTGCCGGGGGC‑TGC‑AA‑C‑CCCGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCtt < 1:293933/130‑1 (MQ=255) tCTGGACCCTGCCGGGCCC‑ACACCTGCAGTGGCTGTTGGGC‑CT‑TTGAGGGGGGC‑TG‑‑AA‑A‑GCAGCCCCAATGCCCTTGATACAAATCAACACTGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCtt < 2:363588/135‑1 (MQ=255) gctgcccccGGCCGGCCCC‑ACACCAGCACTGGCTGAGGGCC‑CA‑TTGCCGGGGGC‑TCC‑AA‑C‑CCCGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCtt < 1:289251/130‑1 (MQ=255) cgggccccGGCCGGGCCC‑ACACCAGCAGTGGC‑GTTGGGG‑CATTTGCGGGGGGG‑TGA‑AA‑C‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCtt < 2:152457/130‑1 (MQ=255) cTGGACCCTGCCGGGCCC‑ACACCAGCACTGGCTGATGGGC‑CA‑TTGACGGGGGC‑TGA‑AA‑A‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCttt > 8:285149/1‑136 (MQ=255) gcccccTGCCGGCCCC‑ACACCAGCACTGGCGGATGGGC‑CA‑TTGACGGGGGC‑GGA‑AA‑C‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTcc < 6:278735/133‑1 (MQ=255) gcaCCCTGCCGGGCCC‑ACACCAGCAGTGGCTTAGGGGC‑AA‑TTGCCGGGGGC‑TGC‑AA‑A‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTcc < 2:54758/134‑1 (MQ=255) aCCCGGCCGGGCCC‑ACACCAGCACGGGCTGAGGGGC‑CA‑TTGACGGGGGC‑TGC‑TA‑A‑GCACCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCgg < 3:399663/136‑1 (MQ=255) ccgtGCGGGGCCA‑ACCCCAGCATTGGGTTTTGGGC‑AA‑TTGCCGGGGGC‑TGA‑AA‑A‑GCACCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCggg < 8:116155/133‑1 (MQ=255) ccTCCCGGCCCC‑ACACACGCTGTGGTTTTTGGG‑‑CA‑ATGCCGGGGGCATCC‑AA‑C‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCgggg < 2:46055/136‑1 (MQ=255) cggCCGGGCCC‑ACACCAGCAGGGGCTTTGGGGC‑CA‑TTGCGGGGGGC‑TCA‑AA‑A‑GCACCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGa < 4:122772/134‑1 (MQ=255) ggCCGGGCCC‑ACCCCCCGCCTGGGTGTTGGGC‑CA‑TTGCCGGGGGC‑TTC‑AA‑ACCCCCCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGa < 8:218156/135‑1 (MQ=255) gcgggcccc‑ACACCAGCTGTGGGTGTTGGCC‑AT‑TTGCGGGGGGC‑TTA‑AA‑C‑CCCCCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACt < 3:399157/130‑1 (MQ=255) ggCCC‑ACACCCGCACTGGCTGAGGGGC‑CA‑TTGCCGGGGGC‑TGC‑AA‑A‑GCCGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTAccc < 8:310086/136‑1 (MQ=255) gcccc‑ACACCCCCCGTGGGTGAGGGCC‑CT‑TTGCGGGGGGC‑TCC‑AA‑C‑GCCGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTAccc < 6:281214/134‑1 (MQ=255) gCATTGGTTGATGGGC‑CA‑TTGCGGGGGGC‑TGA‑AA‑A‑CCAGCCCCAATG‑CCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTcgacg < 1:7679/136‑1 (MQ=255) gCAGTGGCTGATGGGC‑CA‑TTGACGGGGGC‑TGC‑AA‑C‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTcgac < 7:56385/136‑1 (MQ=255) ccaggggCTGTGGGGC‑AA‑TTGCGGGGGGC‑TCC‑AA‑C‑CCACCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTcgac < 3:333799/131‑1 (MQ=255) cagtGGCTGTTGGGC‑CT‑TTGAGGGGGGC‑TGC‑AA‑A‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTcgac < 8:291468/132‑1 (MQ=255) caggggCTGTGGGGC‑AT‑TTGACGGGGGC‑TGC‑AA‑A‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTcgacg < 8:47155/132‑1 (MQ=255) aaCTGGGTGATGGGC‑CT‑TTGAGGGGGGC‑TGC‑AA‑C‑CCACCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTcgac < 4:34116/134‑1 (MQ=255) ggggCTTAGGGGC‑CA‑TGCAGGGGGGC‑TGA‑AA‑A‑CCACCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCgacgac < 4:162824/134‑1 (MQ=255) tGGCTGATGGGC‑CA‑TTGACGGGGGT‑TGC‑AA‑A‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGAcc > 7:116281/1‑136 (MQ=255) gggCTGGTGGGC‑CA‑TTGGCGGGGCC‑CCC‑AACC‑CCCCCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCgacgac < 1:404053/135‑1 (MQ=255) ggtttggggC‑CT‑TTGCGGGGGGC‑GCA‑AA‑A‑CCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCt < 7:14570/129‑1 (MQ=255) gggggtGGGC‑CC‑TTTGCGGGGGG‑CGCAAA‑A‑CCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCt < 2:176272/131‑1 (MQ=255) aggggC‑CA‑TTGACGGGGGG‑TGA‑AA‑A‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCg < 8:220736/132‑1 (MQ=255) gcc‑cA‑TGGAGGGGGGC‑TGC‑AC‑C‑GCGGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTggg < 2:120450/134‑1 (MQ=255) tgcggggggC‑TGA‑AA‑A‑CCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATAcc < 4:126627/131‑1 (MQ=255) | TGCCGAACGGATGAACGCCCTGTTGCGGCGACGCAGCGTGCGAGTGGCAGGCGGTGTGCTGGTGATGCTATTCGGTATCTGGACCCTGCCGGGCCC‑ACACCAGCACTGGCTGATGGGC‑CA‑TTGACGGGGGC‑TGC‑AA‑A‑GCAGCCCCAATGCCCTTGATACAAATCAACACAGTTTCCTACAGACGGCCATAGACTCCGGTGCACTGCTGCTCTTTCCGGGGACTACCCTCATGCTCGACGACCTCCGTTGGGATACC > ACB143/4841788‑4842043 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |