Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,796,838 | A→T | 12.7% | intergenic (+22/‑110) | PP_t43 → / → PP_2449 | tRNA‑Val/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,796,838 | 0 | A | T | 12.7% | 150.8 / 11.8 | 55 | intergenic (+22/‑110) | PP_t43/PP_2449 | tRNA‑Val/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (30/18); new base T (3/4); total (33/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.19e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01 |
AAAGAGCTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCC > NC_002947/2796690‑2796984 | aaaGAGCTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGaa < 7:97818/149‑1 (MQ=255) gCTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGAt < 3:163814/149‑1 (MQ=255) cGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGc < 5:364412/149‑1 (MQ=255) cATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATGCAAAAAGGGCGATCATTTCTGATCGCCCtttttcgat < 3:256393/149‑5 (MQ=17) cATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCAt < 3:313345/149‑1 (MQ=255) tCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCAttt < 4:125645/149‑1 (MQ=255) gggCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATGCAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCGATTTgg < 4:226387/149‑2 (MQ=17) aCTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGccc > 4:354778/1‑148 (MQ=255) aaTATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAag > 3:126319/1‑149 (MQ=255) cacaCGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAAc > 8:273528/1‑149 (MQ=255) cacaCGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAAc > 7:35253/1‑149 (MQ=255) cacGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGt < 8:35253/149‑1 (MQ=255) cacGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGt < 5:144166/149‑1 (MQ=255) acGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCATACCAAATCGAAAAAGCGCGATCAGAAATGATCGCCCTTATTGCATTTACGCCCGAACCAGCGTGTAACGta > 4:95426/1‑148 (MQ=18) cgtcAGGCACATAGCTCAGTTGGTTAGATCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGttt < 8:90914/146‑1 (MQ=255) gCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTa > 6:332814/1‑137 (MQ=255) gCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTa < 5:332814/137‑1 (MQ=255) aTAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCAt < 6:9115/112‑1 (MQ=255) aTAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCAt > 5:9115/1‑112 (MQ=255) aGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCGATTTg > 8:71196/1‑108 (MQ=14) aGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCGATTTg > 8:66082/1‑108 (MQ=14) aGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCGATTTg < 7:66082/108‑1 (MQ=14) aGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCGATTTg < 7:71196/108‑1 (MQ=14) gTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATAcccc > 3:255305/1‑149 (MQ=255) gagCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGc < 3:95426/148‑1 (MQ=255) ccaccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa < 7:77126/149‑1 (MQ=255) ccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGtt > 4:253248/1‑114 (MQ=255) ccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGaa > 5:131179/1‑149 (MQ=255) ttGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATc > 1:275250/1‑149 (MQ=255) gACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATcc < 4:126319/148‑1 (MQ=255) ggtggGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCt > 2:156751/1‑149 (MQ=255) gtggGGGTCGTTGGTTCGAGTCCAATTGTGCATACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATAccc > 7:201293/1‑126 (MQ=255) gTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGCAACGTTTGGTCGCGGTCCg < 2:196694/111‑1 (MQ=255) gTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGCAACGTTTGGTCGCGGTCCg > 1:196694/1‑111 (MQ=255) gTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTgg < 4:255305/149‑1 (MQ=255) ttGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGa < 3:354778/149‑1 (MQ=255) cAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTc > 1:193155/1‑128 (MQ=255) cAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTc < 2:193155/128‑1 (MQ=255) cAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTg > 7:22174/1‑149 (MQ=255) gCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAAc > 5:122552/1‑149 (MQ=255) tACCAAATCGAAAAAGGGCGATCAGAAATGCTCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATACGCACCAAAGAAACAGGAGGAAAACGCGTGCGCAAAcaa > 7:232701/1‑146 (MQ=255) ccAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCt > 2:166985/1‑149 (MQ=255) cAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCtg > 3:312315/1‑149 (MQ=255) cAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCtg > 8:262263/1‑149 (MQ=255) aaTCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCtgtg < 7:273528/149‑1 (MQ=255) aaaaaGGGCGATCATTTCTGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAAt > 8:226564/1‑149 (MQ=255) aaaaaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAAt > 4:194556/1‑149 (MQ=255) aaaaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATg > 4:91086/1‑149 (MQ=255) aaaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATgg > 6:14162/1‑149 (MQ=255) aaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATgg > 7:113307/1‑148 (MQ=255) aaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATgg > 6:333147/1‑148 (MQ=255) gggCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAg > 3:64973/1‑149 (MQ=255) gCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTg > 5:181046/1‑108 (MQ=255) gCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTg < 6:181046/108‑1 (MQ=255) cAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGt > 4:20682/1‑108 (MQ=255) gAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGccc > 8:151169/1‑149 (MQ=255) | AAAGAGCTTTACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGCCC > NC_002947/2796690‑2796984 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |