| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP009273 | 352,394 | A→G | 22.4% | E171E (GAA→GAG) | codA → | cytosine/isoguanine deaminase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP009273 | 352,394 | 0 | A | G | 22.4% | 23.4 / 2.0 | 18 | E171E (GAA→GAG) | codA | cytosine/isoguanine deaminase |
| Reads supporting (aligned to +/- strand): ref base A (4/9); new base G (2/2); total (6/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.84e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 | |||||||||||
CATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGTGAAGCAGGAAGTCGCGCCGTGGATTGATCTGCAAATCGTCGCCTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCG > CP009273/352248‑352510 | cATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGTGAAGCAGGAAGTCGCGCCGTGGATTGATCTGCAAATCGTCGCCTTCCCTTAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAaga < 2:214028/149‑1 (MQ=255) tCGATGTTTCGGATGCAACGCTAACCGCGCTGAAAGCAATGCTGGAAGTGAAGCAGGAAGTCGCGCCGTGGATTGATCTGCAAATCGTCGCCTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAag < 1:96251/144‑1 (MQ=255) tttCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGTGAAGCAGGAAGTCGCGCCGTGGATTGATCTGCAAATCGTCGCCTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGt < 1:160224/144‑1 (MQ=255) ggATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGTGAAGCAGGAAGTCGCGCCGTGGATTGATCTGCAAATCGTCGCCTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACg > 1:231529/1‑144 (MQ=255) aTGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGTGAAGCAGGAAGTCGCGCCGTGGATTGATCTGCAAATCGTCGCCTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAGGAGGCGTTACGCTTAggg > 2:42744/1‑149 (MQ=255) aTGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGTGAAGCAGGAAGTCGCGCCGTGGATTGATCTGCAAATCGTCGCCTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAGGAGGCGTTACGCTTAggg > 2:179735/1‑149 (MQ=255) gcggAAAGCAATGCTGGAAGTGAAGCACGAAGTCGCGCCGTGGATTGATCTGCAAATCGTCGCCTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTggggg < 2:186354/146‑1 (MQ=255) ttgAAGTGAAGCAGGAAGTCGCGCCGTGGATGGATCGGCAAGACGTTACCTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCAtttt < 2:57951/147‑1 (MQ=255) aaGTGAAGCAGGAAGTCGCGCCGTGGATTGCTCTGCCATCCTTTGCCTTCCCTTAGGAAGGGATTTTGTCGTATCCCAACGGTTAAGCTTTGCTGGAAGAGGCGTTACGCTTAGGGGTAGATGTAGTGGGGGCGATTCCGCATTTTGaa < 2:45590/149‑1 (MQ=255) aGCAGGAAGTCGCGCCGTGGATTGATCTGCAAATCGTCGCCTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGAGGTAGTGGGGGCGATTCCGCATTATGaag > 1:98120/1‑143 (MQ=255) aGGAAGTCGCGCCGTGGATTGATCTGCAAATCGTCGCCTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTATGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGt < 2:216003/149‑1 (MQ=255) tGGATTGATCTGCAAATCGTCGCCTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGTAGATGTAGTGGTGGCGATTTCGCATTTTGAATTTACCCGGGAATACGGCGTGGa < 2:17994/149‑1 (MQ=255) gcttggCAAATCGTCGGGTTCCCTCAGGAAAGGATTTTGTCGTATCCCACCGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAg < 1:42723/139‑1 (MQ=255) aTCTGCAAATCGTCGCCTCCCCTCAGGAAGGAATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACTTTTAGGGGTAGTTTTATTGGTTTCTATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTg < 2:95901/149‑1 (MQ=255) cTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGATGAAGCATTACTGAAAGAGGCGATACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGccc > 2:17185/1‑149 (MQ=255) ggAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGACGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGTGTTCGCTGCATAAGACCTTCGCCCTgga < 1:144076/144‑2 (MQ=255) gAAGGTATTTTGTCGTATCCCAGCGGTGAAGCGTTGCTGGAAGAGGCGGTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGgcg < 1:152136/144‑1 (MQ=255) tCGTATCCCAACGGTGAAGCGTTGCTGGAGGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGAc < 1:162523/144‑1 (MQ=255) tCGGATCCCCACGGTGAAGCGTTGCTGGAGGAGGGGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGAc < 1:39241/144‑1 (MQ=255) gTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAGAAATACGACCg > 1:160983/1‑144 (MQ=255) | CATGTCGATGTTTCGGATGCAACGCTAACTGCGCTGAAAGCAATGCTGGAAGTGAAGCAGGAAGTCGCGCCGTGGATTGATCTGCAAATCGTCGCCTTCCCTCAGGAAGGGATTTTGTCGTATCCCAACGGTGAAGCGTTGCTGGAAGAGGCGTTACGCTTAGGGGCAGATGTAGTGGGGGCGATTCCGCATTTTGAATTTACCCGTGAATACGGCGTGGAGTCGCTGCATAAAACCTTCGCCCTGGCGCAAAAATACGACCG > CP009273/352248‑352510 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |