Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | AM260480 | 1,192,562 | +C | coding (163/594 nt) | cccA → | cytochrome c, mono‑or diheme variant |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AM260480 | 1,192,560 | 1 | . | C | 95.5% | 79.1 / ‑3.3 | 22 | coding (161/594 nt) | cccA | cytochrome c, mono‑or diheme variant |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (10/11); total (10/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GT‑CCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATG > AM260480/1192435‑1192688 | cGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGAcc < 1:711282/134‑1 (MQ=255) gctggTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTc > 1:545541/1‑134 (MQ=255) tATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGc > 2:51068/1‑134 (MQ=255) tATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGc > 2:340522/1‑134 (MQ=255) tATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGc > 1:35736/1‑134 (MQ=255) tGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATg > 1:859476/1‑134 (MQ=255) aaaCAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGAc < 2:314721/134‑1 (MQ=255) aacaacCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGa < 1:1007275/123‑1 (MQ=255) aacaacCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGa > 2:1007275/1‑123 (MQ=255) ggAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAgg > 1:101180/1‑113 (MQ=255) ggAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAgg < 2:101180/113‑1 (MQ=255) tCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCccgg‑cc < 2:1015227/69‑1 (MQ=255) tCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCccgg‑cc > 1:1015227/1‑69 (MQ=255) aTGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCccgg‑cc < 1:687577/66‑1 (MQ=255) aTGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCccgg‑cc > 2:687577/1‑66 (MQ=255) aaaCGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGcgac < 2:35736/134‑1 (MQ=255) cgccTGGTTGCCTTGCTCGCCCCTTGCGCGGC‑CCCTGCCGCCCCCCCGGCCCGGCT‑CGGGCTGGCCCATGGCGGCCTCCCGCCCCCGGCCGTCGGAACCCAAGCGTTACCACAACTGGGCGACGAATGGCGccc < 1:51068/134‑1 (MQ=255) cgccTGCTTGCCTTGCTCGCGCTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGccc < 2:859476/134‑1 (MQ=255) cgcTTTGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGc < 2:545541/134‑1 (MQ=255) tGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGt < 1:694452/132‑1 (MQ=255) tGCGC‑GTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGt > 2:694452/1‑132 (MQ=255) cgc‑gTCCCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATg < 1:340522/134‑1 (MQ=255) | CGGGCCTGGTCCGGCTGGCTGGTATGGCTCTTGCGGAAACAGGGCAAACAACCAAGACCCGTCCACTGTAGGAGACATCCATGAAAACGCAGCGCCGCCTGCTTGCCTTGCTCGCGCTTTGCGC‑GT‑CCTGACCGCCCTCCCGG‑CCGGCTCCGTGCTGGCCCATGGCGACGTCACGCCCCAGGCCGTCGATACCAAAGCGTTACCACAACTGGGCGACGAATGGCGCCCGGACAATCCCTACCGTAGCGGTGATG > AM260480/1192435‑1192688 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |