Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,012,931 | 0 | A | C | 50.0% | 16.1 / 12.9 | 26 | T138P (ACA→CCA) | ssnA | putative aminohydrolase SsnA |
Reads supporting (aligned to +/- strand): ref base A (1/12); new base C (13/0); total (14/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.69e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTC‑GCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGC > NZ_CP009273/3012848‑3013018 | gTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGG‑GCGCTCCCCAtttc > 2:215848/1‑88 (MQ=255) gAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTC‑GCTCTCCACATTGCGC‑GACG‑Ca < 1:1208755/90‑1 (MQ=255) aaGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTC‑GCTCTCCACATTGCGC‑GACG‑CAt < 1:23991/90‑1 (MQ=255) aGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGGCCCCCCCCATT‑CGC‑GACG‑ccct > 2:572192/1‑87 (MQ=255) aTTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTC‑GCTCTCCACATTGCGC‑GACG‑CATTTTTa < 1:428925/90‑1 (MQ=255) aTTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTC‑GCTCTCCACATTGCGC‑GACG‑CATTTTTa < 1:926772/90‑1 (MQ=255) cGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGG‑GCCCTCCACCCCGCGCGGCCC‑CATTTTTTAAAATTg > 2:864041/1‑90 (MQ=255) aTCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGG‑GCTCTCCCCATTTCGG‑GGCG‑CGTTTTTTAAAATTGGGCTgcgcgc > 1:409381/1‑90 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGG‑GCTCTCCCCCATTCGC‑GGCG‑CCTTTTTTTAAATTTGGCTGcgcgcggg > 1:1151980/1‑88 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGC‑GGGCTCCCCCTTCCCC‑GGCG‑CCCTTTTTTAAATTGGCCTGcgcgcggg > 2:589075/1‑88 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGC‑GCCCTCCCCATTGCGC‑GACG‑CCTTTTTTAAAATTTGGCTGcgcgcggt > 1:610632/1‑88 (MQ=255) tATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGC‑GCTCTCCCCATTTCGC‑GGCG‑CCTTTTTTAAAATTTGTCCGCGCGCGGTGAc > 2:73168/1‑90 (MQ=255) aTCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTC‑GCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGAcc < 2:965286/90‑1 (MQ=255) aTCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTC‑GCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGAcc < 2:409381/90‑1 (MQ=255) tCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTC‑GCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCt < 2:610632/90‑1 (MQ=255) tCACCATGCCTCTCCGGCGTATATCGGCGGGTC‑GCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGttt < 2:32520/90‑1 (MQ=255) aTGCCTCTCCGGCGTATATCGGCGGGGG‑GCCCTCCCCCCTTCCC‑GACG‑CCTTTTTTAAAATTGGCCTGGGCGGGGGGGCCCGTTTTGgaa > 2:867695/1‑90 (MQ=255) gCCTCTCCGGCGTATATCGGCGGGTC‑GCTCTCCCCATTGCGC‑GACG‑CATTTTTTAAAATTGGGCTGCGCGCGGTGACCCGTTTTTAAact > 1:435388/1‑90 (MQ=255) ctcCGGCGTATATCGGCGGGTC‑GCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTg < 1:821755/90‑1 (MQ=255) ctcCGGCGTATATCGGCGGGGC‑GCTCTCCCCACTGCGC‑GCCGCCCTTTTTAAAAATGGGCC‑GCGCGCGGGGGCCCGTTTTTTAaatactt > 2:299748/1‑89 (MQ=255) ccGGCGTATATCGGCGGGGG‑GCTCTCCCCATTTCGC‑GCCG‑CCCTTTTTTAAATTGGGCCGCGCCCGGTGGCCCCTTTTTTAaatacttaa > 1:1324380/1‑87 (MQ=255) gggCGTATATCGGCGGGTC‑GCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGAcc < 1:1350952/89‑1 (MQ=255) gCGTATATCGGCGGGGC‑GCTCTCCCCCTTTCGC‑GGCG‑CCTTTTTTAAAATTGGCCTGCGGGCGGGGGCCCGTTTTTTTAATACTTACCGt > 2:1257738/1‑90 (MQ=255) gTATATCGGCGGGGC‑GCCCTCCCCCCTTCGC‑GCCC‑CCTTTTTTAAAAATTGCCTGCGCGGGGGGGCCTGGTTTTTAAATACTGACCCTaa > 2:1222251/1‑90 (MQ=255) gTC‑GCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCaa < 2:348820/90‑1 (MQ=255) ccACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGc < 1:1079658/90‑1 (MQ=255) | GTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTC‑GCTCTCCACATTGCGC‑GACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAACGGCATCAAAGAGTTGC > NZ_CP009273/3012848‑3013018 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |