Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 956,588 | A→C | 14.7% | intergenic (+18/‑96) | alaC → / → PP_0818 | aminotransferase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 956,588 | 0 | A | C | 14.7% | 84.1 / 11.8 | 41 | intergenic (+18/‑96) | alaC/PP_0818 | aminotransferase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (20/15); new base C (2/4); total (22/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.90e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01 |
GGGGATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTC > NC_002947/956441‑956731 | ggggATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAagc < 1:64297/150‑1 (MQ=255) ggggATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAagc < 1:642370/150‑1 (MQ=255) gATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACGGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAagcagc > 1:147400/1‑150 (MQ=255) cTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCatt < 1:637080/150‑3 (MQ=255) cTTTGGTGACTAAGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGACCATCCAATGCGGCTGCTTCGCAGCCCCatt < 1:337015/150‑3 (MQ=255) gTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGcat > 1:476754/1‑148 (MQ=255) tATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTgg < 1:929/150‑1 (MQ=255) tATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTgg > 1:203075/1‑150 (MQ=255) gggACGACAGTGTGCGCGTTGCGTTGATCGAGAAGCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTCCCGGGCTGATGGTTGGGGTCGCAAGCAATCCTAATGGGGCTGCGAAGCAGCCCCAACGGCCTCTTGTGTGa < 2:1073976/150‑1 (MQ=255) tgCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCgcctg > 1:1091515/1‑150 (MQ=255) tGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTg < 2:1091515/150‑1 (MQ=255) gTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCatt < 2:277349/114‑3 (MQ=255) gTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCatt > 1:277349/1‑112 (MQ=255) gTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGtt > 1:704766/1‑150 (MQ=255) aTCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCattgg > 2:114935/1‑108 (MQ=255) aTCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCTTCGCAGCCCCattgg < 1:114935/112‑5 (MQ=255) aTCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTAccc < 2:203075/150‑1 (MQ=255) tCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGa < 1:276040/123‑1 (MQ=255) tCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGa > 2:276040/1‑123 (MQ=255) ccAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTg < 2:1002229/135‑1 (MQ=255) ccAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTg > 1:1002229/1‑135 (MQ=255) cAGGAACGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTa > 1:224911/1‑150 (MQ=255) gACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCg > 1:1064470/1‑150 (MQ=255) tATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGc > 1:987877/1‑150 (MQ=255) aTTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGaa < 2:831255/81‑1 (MQ=255) aTTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGaa > 1:831255/1‑81 (MQ=255) gCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCa > 2:557333/1‑150 (MQ=255) cAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCGTCTGGTTACCCGTACGCCTACCGTCAGCACTCAtt > 1:264157/1‑150 (MQ=255) gcgGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCaaa < 2:147400/150‑1 (MQ=255) ggATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCTCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTAc < 2:793552/123‑1 (MQ=255) ggATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTAc > 1:793552/1‑123 (MQ=255) tCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATc > 1:1082562/1‑150 (MQ=255) gCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATccccc > 1:824409/1‑150 (MQ=255) gTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAgt < 2:224911/150‑1 (MQ=255) gTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAgt < 2:704766/150‑1 (MQ=255) ccGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTc < 2:1064470/150‑1 (MQ=255) ccGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTc > 1:725086/1‑150 (MQ=255) gggCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGc > 1:926317/1‑101 (MQ=255) gggCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGc < 2:926317/101‑1 (MQ=255) gCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCa > 2:121109/1‑150 (MQ=255) gTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTa > 1:64909/1‑150 (MQ=255) tCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGc < 1:557333/150‑1 (MQ=255) gCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCt > 1:51204/1‑150 (MQ=255) cTGCGAAGTAGCCCAAACTGCCGCTTGGGTGCGCTGTCAGCCGGAGTTGGGGGGAGCCGTAGGCCGACCGTGAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTc < 2:264157/150‑1 (MQ=255) gAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAgt > 1:785849/1‑104 (MQ=255) gAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAgt < 2:785849/104‑1 (MQ=255) | GGGGATCGGCTTTGGTGACTATGGGGACGACCATGTGCGCTTTGCGTTGATCGAGAACCAGGACCGTATTCGCCAGGCCATCCGCGGGATCCGGCAAATGTTCCGGGCTGATGGTTTGGCTCGCAAGTAAGCCCAATGGGGCTGCGAAGCAGCCCCAACTGCCTCTTGGGTGAGCTTTCCGCCTGCCTTCTGGTTACCCGTACGCCTACCGTCAGCACTCATTTCTCACCAAAGAGATCCCCCATGAGTGTCTTCACCGCCTACTTCTGCGGCACTGGCTCGCACCGTTTC > NC_002947/956441‑956731 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |