Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,799,638 | T→C | 12.3% | V65A (GTA→GCA) | PP_5664 → | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,799,638 | 0 | T | C | 12.3% | 77.8 / 3.4 | 33 | V65A (GTA→GCA) | PP_5664 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (14/14); new base C (2/2); total (17/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01 |
GATTTTCAGGCAGGAGGTGCCTCCTTCTAGCATAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACAA > NC_002947/4799489‑4799777 | gATTTTCAGGCAGGAGGTGCCTCCTTCTAGCATAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGt > 2:1282419/1‑150 (MQ=255) ggggcgtcctTCTAGCATAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAAGACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCtcattcat < 1:81968/144‑1 (MQ=255) tCTAGCATAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACt < 2:86970/150‑1 (MQ=255) gCATAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTgcgc > 2:1029001/1‑150 (MQ=255) tAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCg < 1:119390/150‑1 (MQ=255) gcTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGc < 2:401717/150‑1 (MQ=255) aaTTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATc < 1:422692/150‑1 (MQ=255) ccGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAg > 2:686420/1‑112 (MQ=255) ccGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAg < 1:686420/112‑1 (MQ=255) gaagaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAgg < 2:408102/114‑1 (MQ=255) gaagaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAgg > 1:408102/1‑114 (MQ=255) aagaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTa < 1:1239262/93‑1 (MQ=255) aagaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTa > 2:1239262/1‑93 (MQ=255) aagaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATc > 1:681227/1‑150 (MQ=255) gaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTAt > 1:370878/1‑139 (MQ=255) gaaACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTAt < 2:370878/139‑1 (MQ=255) gAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGcc > 2:424676/1‑150 (MQ=255) gACCGAGCACGGATGGTGCGGACATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGCAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGa > 2:227730/1‑150 (MQ=255) gACCGAGCACGGATGGTGCGGACATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGCAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGa > 2:228030/1‑150 (MQ=255) aCCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGAt > 2:175217/1‑150 (MQ=255) aCCGAGCACGGATGGTGCGGATATGGTGATGGAATATCCGGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGAt < 1:1231755/150‑1 (MQ=255) ggATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTg > 2:459875/1‑150 (MQ=255) tGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTAc < 1:175217/150‑1 (MQ=255) gTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACgg > 1:984792/1‑150 (MQ=255) atatGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGacaa < 1:1282419/150‑1 (MQ=255) tGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTgg > 2:345780/1‑107 (MQ=255) tGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTgg < 1:345780/107‑1 (MQ=255) gTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTa > 2:487690/1‑150 (MQ=255) tCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAg > 1:1041461/1‑150 (MQ=255) cAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGc > 2:992768/1‑150 (MQ=255) cAAGACTACGAGGGCGGAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACt > 2:690564/1‑150 (MQ=255) gACTACGAGGGCGCAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGCAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTca < 1:228030/150‑1 (MQ=255) gACTACGAGGGCGCAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTca < 1:227730/150‑1 (MQ=255) tACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACaa < 2:1041461/150‑1 (MQ=255) | GATTTTCAGGCAGGAGGTGCCTCCTTCTAGCATAGAGCGCTACAGCGGTAAATTGCCGAAGAAACTGTTGGAATACTGGACCGAGCACGGATGGTGCGGATATGGTGATGGAATATTCTGGACGGTCAATCCGCAAGACTACGAGGGCGTAGCCGCCTCATTCATAGAAGGCACTGCGCTCGAAGGGCGAGATACCTATCACGTCATCGCTCGGGGCGCCTTTGGTGATCTCTATCTGTACGGAGAAAATTCCGGTTTTTCATTAGACATCGTAGCCCACACTTCACAA > NC_002947/4799489‑4799777 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |