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 | 89.3% | 68.6 / 2.4 | 28 | F153S (TTT→TCT) | folA | type 3 dihydrofolate reductase |
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (5/20); total (6/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.30e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 |
CCCGGATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTCCTA > NZ_CP009273/50197‑50366 | cccGGATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGAt < 1:140533/90‑1 (MQ=255) aGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAgcggc < 2:59928/90‑1 (MQ=255) cGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGt < 2:544030/90‑1 (MQ=255) cGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGt < 1:63846/90‑1 (MQ=255) gATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTaa < 2:304225/90‑1 (MQ=255) gACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt < 2:29220/90‑1 (MQ=255) cGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt < 1:53465/80‑1 (MQ=255) cGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt < 1:560754/80‑1 (MQ=255) cGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAAttt > 2:53465/1‑80 (MQ=255) cGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGAAAttt > 2:560754/1‑80 (MQ=255) tATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTa < 1:355360/90‑1 (MQ=255) tCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCg < 1:836273/90‑1 (MQ=255) cTGATGCGCAGAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTTCGGCTAACCCCCGATGCGACg > 1:521429/1‑90 (MQ=255) gAACTCTCACAGCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGt < 2:394907/90‑1 (MQ=255) aCTCTCACAGCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTTCGGCTAGCGCCGGGTGCGGCGCCGGTCGCGTCt > 2:55649/1‑90 (MQ=255) gCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc < 2:418934/90‑1 (MQ=255) gCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc < 1:188844/90‑1 (MQ=255) gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc < 1:100285/90‑1 (MQ=255) gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc < 2:500025/90‑1 (MQ=255) gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc < 1:243547/90‑1 (MQ=255) gCTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGcc < 2:740628/90‑1 (MQ=255) cTATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCt < 2:367833/90‑1 (MQ=255) tATTGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGGTGGGACGCCGGGCGCGTCTTTTCCGGCCtt > 1:653008/1‑90 (MQ=255) ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATATAATTTACGGCTAGCGCCGGAGGCGACGCCGGGCGCGTCTTATCCGGCCCTcc > 1:549904/1‑90 (MQ=38) ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTcc < 1:740430/90‑1 (MQ=255) ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTcc < 1:412143/90‑1 (MQ=255) ttGCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTcc < 2:776429/90‑1 (MQ=255) gCTCTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTCCta < 2:785102/90‑1 (MQ=255) | CCCGGATTACGAGCCGGATGACTGGGAATCGGTATTCAGCGAATTCCACGATGCTGATGCGCAGAACTCTCACAGCTATTGCTTTGAGATTCTGGAGCGGCGGTAATTTTGTATAGAATTTACGGCTAGCGCCGGATGCGACGCCGGTCGCGTCTTATCCGGCCTTCCTA > NZ_CP009273/50197‑50366 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |