Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,051,619 | G→C | 12.1% | P315A (CCG→GCG) | PP_4452 ← | NAD/NADP octopine/nopaline dehydrogenase family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,051,619 | 0 | G | C | 12.1% | 90.0 / 6.0 | 33 | P315A (CCG→GCG) | PP_4452 | NAD/NADP octopine/nopaline dehydrogenase family protein |
Reads supporting (aligned to +/- strand): ref base G (13/16); new base C (2/2); total (15/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01 |
GTATCAGCTTACGCTGCAGCTACACGCGCCAGTAGCCCTTCCATAGACTCGGACTCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCGGGTGTCAACTTCAATAGGCCCGGGCGGCCCACCTCTCTTCTGGTGTAGTTCCTCCGCAATTTCTTT > NC_002947/5051474‑5051760 | gTATCAGCTTACGCTGCAGCTACACGCGCCAGTAGCCCTTCCATAGACTCGGACTCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATc < 2:167468/150‑1 (MQ=255) ttACGCTGCAGCTACACGCGCCAGTAGCCCTTCCATAGACTCGGACTCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACAtt < 1:686399/150‑1 (MQ=255) cGCTGCAGCTACACGCGCCAGTAGCCCTTCCATAGACTCGGACTCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGc < 2:326226/150‑1 (MQ=255) cTACACGCGCCAGTAGCCCTTCCATAGACTCGGACTCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGcc < 2:1135085/150‑1 (MQ=255) tACACGCGCCAGTAGCCCTTCCATAGACTCGGACTCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGccc > 2:689515/1‑150 (MQ=255) acGCGCCAGTAGCCCTTCCATAGACTCGGACTCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCaa > 1:191816/1‑150 (MQ=255) gcCAGTAGCCCTTCCATAGACTCGGACTCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGcc > 2:888005/1‑150 (MQ=255) agacgCCCTTCCATCGATCCGGACTCAACGTTCAGTCCGTCGACGATAGCGGTTTCGTTAGAGAAACCGTCTGGCACCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCaaa < 2:524432/146‑1 (MQ=255) aGTAGCCCTTCCATAGACTCGGACTCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCaaa < 1:888005/150‑1 (MQ=255) ccATAGACTCGGACTCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGaa > 1:889860/1‑145 (MQ=255) ccATAGACTCGGACTCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGaa < 2:889860/145‑1 (MQ=255) tCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACAtt < 1:767560/104‑1 (MQ=255) tCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACAtt > 2:767560/1‑104 (MQ=255) tcgaCCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGTTGTAt > 1:1204270/4‑150 (MQ=255) aCCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACTCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCgg < 2:5664/150‑1 (MQ=255) aCCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCgg > 1:339279/1‑150 (MQ=255) aGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCt > 2:1090510/1‑87 (MQ=255) aGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCt < 1:1090510/87‑1 (MQ=255) tttCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCGGGTGTCAACtt < 2:191816/150‑1 (MQ=255) gTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACCTCGGCTGTACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCGGGTGTCAACTTCAATAGGCCCGggcgg < 2:36441/150‑1 (MQ=255) ggCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCGGGTGTCAACTTCAATAGGCCCGGGCGGCCCa > 2:825887/1‑150 (MQ=255) gcagcaGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCGGGTGTCAACTTCAATAGGCCCGGGCGGCCCAc < 1:1133559/150‑1 (MQ=255) gggCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCAGGTGTCAACTTCAATAGGCACCTGCGGCCCACCTCtcttct > 1:901762/1‑150 (MQ=255) gCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTAACGTGGATCAATTTCAATAGGCCCTGGGTGCCCACCTCTCTTCTgg > 2:968213/1‑150 (MQ=255) aacaaTACTTTGGGTGACCGCCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGt > 2:1203197/1‑92 (MQ=255) aacaaTACTTTGGGTGACCGCCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGt > 2:598241/1‑92 (MQ=255) aacaaTACTTTGGGTGACCGCCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGt < 1:598241/92‑1 (MQ=255) aacaaTACTTTGGGTGACCGCCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGt < 1:1203197/92‑1 (MQ=255) acaaTACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCGGGTGTCAACTTCAATAGGCCCGGGCGGCCCACCTCTCTTCTGGTGTAGTtcctcc < 1:825887/150‑1 (MQ=255) caaTACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCGGGTGTCAACTTCAATAGGCCCGGGCGGCCCACCTCTCTTCTGGTGTAGTTCCTCCg > 2:1125495/1‑150 (MQ=255) ggTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCGGGTGTCAACTTCAATAGGCCCGGGCGGCCCACCTCTCTTCTGGTGTAGTTCCTCCGCAAtttcttt < 1:358645/150‑1 (MQ=255) cGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCGGGTGTCAACTTCAATAGGCCCGGGCGGCCCACCTCtcttct < 2:587251/119‑1 (MQ=255) cGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCGGGTGTCAACTTCAATAGGCCCGGGCGGCCCACCTCtcttct > 1:587251/1‑119 (MQ=255) | GTATCAGCTTACGCTGCAGCTACACGCGCCAGTAGCCCTTCCATAGACTCGGACTCAAGGTTCAGTCCCTCAACCAGATCGTTTTCGTTAGAGAAACCGTCTGGCAGCAGTAGGGCCGACATTGCAACAATACTTTGGGTGACCGGCATCGGTACATTGGCTATACGCCCCAACGCCAAAGTGAATTCCAAACCGTATGGCACGTCTTCAGAGATGTATCGGGTGTCAACTTCAATAGGCCCGGGCGGCCCACCTCTCTTCTGGTGTAGTTCCTCCGCAATTTCTTT > NC_002947/5051474‑5051760 |
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: ‑ |