Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 1,021,663 | (A)5→4 | coding (1737/2055 nt) | helD → | DNA helicase IV |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,021,659 | 0 | A | . | 94.7% | 67.7 / ‑1.7 | 19 | coding (1733/2055 nt) | helD | DNA helicase IV |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base . (14/4); total (15/4) | |||||||||||
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 |
GGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATC > NZ_CP009273/1021576‑1021744 | ggATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGaaaagcg > 1:941973/1‑87 (MQ=255) aGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGaaaa > 2:904987/1‑82 (MQ=255) aGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGaaaa < 1:904987/82‑1 (MQ=255) tctGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCt > 1:505011/1‑90 (MQ=255) aaCCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCa > 2:883423/1‑90 (MQ=255) aaCCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCa > 2:302116/1‑90 (MQ=255) gAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCg < 2:783797/90‑1 (MQ=255) gcgcATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACttt > 1:213716/1‑90 (MQ=255) ttCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGa > 2:366308/1‑90 (MQ=255) cTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATg > 1:302517/1‑90 (MQ=255) cTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATg > 1:456072/1‑90 (MQ=255) tGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGc < 2:802090/90‑1 (MQ=255) ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATg > 2:347801/1‑78 (MQ=255) ggCGCGTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATg < 1:347801/78‑1 (MQ=255) gTTACCATCACATGAGGCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCaa > 2:633798/1‑90 (MQ=255) ggCCTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGa > 2:982299/1‑90 (MQ=255) cTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTa > 1:975323/1‑90 (MQ=255) cTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTa > 2:621202/1‑90 (MQ=255) cTGCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTa < 1:621202/90‑1 (MQ=255) gCCAGCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACg > 1:939881/1‑90 (MQ=255) gCCTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGtcat > 1:77483/1‑90 (MQ=255) ccTGG‑AAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGtcatc > 2:698394/1‑90 (MQ=255) | GGATAAGCTCTCTGGTTATGCCAAACCGGAAGAGCGCATTCTGATCCTGGCGCGTTACCATCACATGAGGCCTGCCAGCCTGGAAAAAGCGGCAACACGCTGGCCGAAGTTGCAAATCGACTTTATGACCATTCATGCCAGCAAAGGGCAACAGGCGGATTACGTCATC > NZ_CP009273/1021576‑1021744 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |