Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 50,280 | T→C | F153S (TTT→TCT) | folA → | type 3 dihydrofolate reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 50,280 | 0 | T | C | 95.8% | 151.1 / ‑3.1 | 48 | F153S (TTT→TCT) | folA | type 3 dihydrofolate reductase |
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (18/28); total (19/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01 |
CCCGGATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGCCTTCCTA > NZ_CP009273/50197‑50366 | cccGGATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGAt > 1:521298/1‑90 (MQ=255) cGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcg < 1:174725/90‑1 (MQ=255) cGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcg > 1:325952/1‑90 (MQ=255) aGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcggc < 1:99408/90‑1 (MQ=255) aGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcggc < 2:383498/90‑1 (MQ=255) ccGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGcggcgg > 2:176032/1‑90 (MQ=255) cGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGt < 1:1117680/90‑1 (MQ=255) cGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGt > 1:759706/1‑90 (MQ=255) gACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt < 1:354375/90‑1 (MQ=255) gACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt < 1:553927/90‑1 (MQ=255) gACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt < 1:176032/90‑1 (MQ=255) cTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTg > 2:1005873/1‑90 (MQ=255) gAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGtata < 1:64081/90‑1 (MQ=255) aTTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAATTTTTATAGAAATTAc > 1:464621/1‑90 (MQ=255) tCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTTTATAGAATTTACgg > 1:1234018/1‑90 (MQ=255) gCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCCGAGATTCTGGAGCGGCGGTAATTTTTTATAGAATTTACGGCTa > 1:624507/1‑90 (MQ=255) aaTTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTTTATAGAATTTACGGCTAgcg > 1:257853/1‑90 (MQ=255) aaTTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTTCGGCTAacg > 2:977355/1‑90 (MQ=255) aaTTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGGAATTTTTTATAGAAATTTCGGCTaacc > 1:337320/1‑87 (MQ=255) aTTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAgcgc < 2:591110/90‑1 (MQ=255) ttCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTTTTTAGAATTTACGGCTAACGcc > 1:934170/1‑90 (MQ=255) aTGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGGTGCg > 2:173593/1‑90 (MQ=255) gATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGcc > 2:86389/1‑90 (MQ=255) cgcAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTc < 2:585628/90‑1 (MQ=255) cgcAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTc < 1:68836/90‑1 (MQ=255) cgcAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTc < 2:271275/90‑1 (MQ=255) aaCTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTc < 1:515838/90‑1 (MQ=255) ctcACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTTTTTATAATTTTCGGCTAGGGCCCGGAGCGGCGCCCGGGGCGGC‑TTAt > 1:1078590/1‑90 (MQ=255) ctcACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGGCGCGTC‑TTAt > 2:98306/1‑90 (MQ=255) aGCTATTGCTTTGAGATTCTGGAGCGGCGGTAAttt < 2:783707/36‑1 (MQ=255) aGCTATTGCTTTGAGATTCTGGAGCGGCGGTAAttt > 1:783707/1‑36 (MQ=255) aGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGc < 2:1234039/90‑1 (MQ=255) gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGcc < 2:381493/90‑1 (MQ=255) gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGcc < 1:1005873/90‑1 (MQ=255) gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGcc < 1:1110193/90‑1 (MQ=255) gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGcc < 1:1144128/90‑1 (MQ=255) gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGcc < 1:173593/90‑1 (MQ=255) gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGcc < 1:220129/90‑1 (MQ=255) gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGcc < 2:1078590/90‑1 (MQ=255) cTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTTTATAGAATTTACGGCTAACGCCGGGTGCGACGCCGGGCGCGGC‑TTTTTCCGCCt > 2:1110193/1‑90 (MQ=255) ttGCTCTGAGATTCTGGCGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGCCTTcc < 1:966831/90‑1 (MQ=255) ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTTTATAGAATTTACGGCTAGGGCCGGGTGGGACGCCGGTGGCGGCTTTTTCCGGCGTTc > 1:59389/1‑90 (MQ=255) ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGCCTTcc < 1:98306/90‑1 (MQ=255) ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGCCTTcc < 2:207817/90‑1 (MQ=255) ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGCCTTcc < 2:627737/90‑1 (MQ=255) ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGCCTTcc < 1:86389/90‑1 (MQ=255) ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGCCTTcc < 2:813002/90‑1 (MQ=255) gCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGCCTTCCta < 1:262499/90‑1 (MQ=255) | CCCGGATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTC‑TTATCCGGCCTTCCTA > NZ_CP009273/50197‑50366 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |