Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,820,642 | Δ1 bp | 35.3% | intergenic (+400/+112) | PP_2473 → / ← PP_2474 | hypothetical protein/glutathione S‑transferase family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,820,642 | 0 | A | . | 35.3% | 31.0 / 18.4 | 17 | intergenic (+400/+112) | PP_2473/PP_2474 | hypothetical protein/glutathione S‑transferase family protein |
Reads supporting (aligned to +/- strand): ref base A (4/7); new base . (3/3); total (7/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.44e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01 |
TTATTGAAGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAG > NC_002947/2820493‑2820732 | ttATTGAAGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATACCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGa > 2:379310/1‑150 (MQ=255) gggATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGt > 1:1778/1‑150 (MQ=255) cAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGtt < 2:1778/150‑1 (MQ=255) aaGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGa < 2:359437/150‑1 (MQ=255) tcGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCg < 1:426430/122‑1 (MQ=255) tcGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCg > 2:426430/1‑122 (MQ=255) gTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCgg < 1:239474/130‑1 (MQ=255) gTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCgg > 2:239474/1‑130 (MQ=255) tCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGtt < 2:48049/114‑1 (MQ=255) tCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGtt > 1:48049/1‑114 (MQ=255) tCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGctc < 1:379310/150‑1 (MQ=255) gggCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGCGCAAAGCTGCTCTCAGACCGGCGTTGGTTTCTTGGa < 2:314223/150‑1 (MQ=255) aaaaCCCGGGTTCG‑ACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAg > 2:131959/1‑104 (MQ=255) aaaaCCCGGGTTCG‑ACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAg < 1:131959/104‑1 (MQ=255) aaCCCGGGTTCG‑ACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAg < 1:123158/102‑1 (MQ=255) aaCCCGGGTTCG‑ACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAg > 2:123158/1‑102 (MQ=255) aCCCGGGTTCG‑ACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTa < 1:429136/56‑1 (MQ=25) aCCCGGGTTCG‑ACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTa > 2:429136/1‑56 (MQ=25) | TTATTGAAGGGATGCCCAGGGGTTAGTTCGAAAGTTTCTCCTCGTGGGTCAGATATATCCCATGATCGGTGGTCCGTGAGCTAGAAGTTGGGCTATTTGCTGGCTGATATAGAGAAAGTGCGTCTGTAAAAAAGCAAAAACCCGGG‑TCGAACCCGGGTTTTCGGTTGAGTCACTCGGTTGAGTCCTTTTGGGGTATGAATTGAGCAATGCTGCTCTCAGACCGGCGTTGGTTTCTTGGAG > NC_002947/2820493‑2820732 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |