Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,777,456:1 | +T | 7.8% | intergenic (+19/+24) | betA‑II → / ← PP_5065 | choline dehydrogenase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,777,456 | 1 | . | T | 7.8% | 187.8 / 11.1 | 51 | intergenic (+19/+24) | betA‑II/PP_5065 | choline dehydrogenase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base . (23/24); new base T (2/2); total (25/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01 |
TGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGCCTGGCGC > NC_002947/5777310‑5777605 | tGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGAGGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCa > 1:1157642/1‑150 (MQ=255) gagaAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCTGC‑GCGGcccc < 2:1113212/150‑1 (MQ=255) aagaTTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑TGCGCGGCCCCtga < 2:1005955/150‑3 (MQ=255) tCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGtt < 2:1300924/150‑1 (MQ=255) cGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGttt < 2:1157642/150‑1 (MQ=255) cGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGttt < 1:144893/150‑1 (MQ=255) ggACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCa > 1:1074545/1‑131 (MQ=255) ggACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCa < 2:1074545/131‑1 (MQ=255) aTCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGc > 2:1324004/1‑150 (MQ=255) ccGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGttt < 2:967711/140‑1 (MQ=255) ccGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGttt > 1:967711/1‑140 (MQ=255) cGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACAAGTGAAAGACAAGTCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGTGGCCGC‑GCAGCTGCCTCTTTCACGTTTGCTCAGGCggg < 2:162557/150‑1 (MQ=255) gCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCa > 2:539539/1‑150 (MQ=255) cTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCGCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtgt < 1:1151052/150‑1 (MQ=255) tGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtgtg > 2:296932/1‑150 (MQ=255) gCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTgc < 2:735359/150‑1 (MQ=255) aCCGCCAAGTGCTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGa < 2:1140846/150‑1 (MQ=255) ccGCCACGTAGTACTTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTTAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCCTCAACGaa > 1:484086/1‑150 (MQ=255) ccGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGaa < 2:848084/150‑1 (MQ=255) ccGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCACCTTCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGaa < 1:175392/150‑1 (MQ=255) tactacGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTc > 2:1034335/1‑150 (MQ=255) tGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACAccc > 1:817214/1‑150 (MQ=255) cGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGt < 1:296932/150‑1 (MQ=255) gATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGccc < 2:787570/150‑1 (MQ=255) aTGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGcccc < 1:70635/150‑1 (MQ=255) gCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGc < 1:1143169/72‑1 (MQ=255) gCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGc > 2:1143169/1‑72 (MQ=255) aCCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTgcg > 2:959945/1‑112 (MQ=255) aCCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTgcg < 1:959945/112‑1 (MQ=255) aCCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTGTCACGTTTGCTCAGGCGGGCAGCGCAAAAACCGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCa > 2:888872/1‑150 (MQ=255) gAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCt < 2:1096121/150‑1 (MQ=255) aaGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAg > 2:1049664/1‑81 (MQ=255) aaGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAg < 1:1049664/81‑1 (MQ=255) aaGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATc > 1:755099/1‑150 (MQ=255) ccGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtgt > 2:47319/1‑93 (MQ=255) ccGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtgt < 1:47319/93‑1 (MQ=255) gATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCa < 2:674428/96‑1 (MQ=255) gATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGAGCa > 1:674428/1‑96 (MQ=255) tGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGa > 2:1141732/1‑141 (MQ=255) tGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGa < 1:1141732/141‑1 (MQ=255) gtgaagtgaaACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGACCGAACAACtct > 2:2103/1‑148 (MQ=255) cTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGc > 1:558106/1‑150 (MQ=255) ttGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGcc > 1:321327/1‑150 (MQ=255) aTCAGGGGCCGC‑GCAGCGGCCCCTTTCACTTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCACCAACGAATCATGCTCGACACCCAGTCGCCCCGCCAGAGCTTGGCGATCGAACACTGGTCCTGTGCGCATTTCACGCACGGCCTTc < 1:888872/150‑1 (MQ=255) aTCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAACAAATATACGCATCAACGAAACGTCTGCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTc > 1:174277/1‑150 (MQ=255) aGGGGCCGCTGC‑GCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGAt < 2:1112921/109‑1 (MQ=255) aGGGGCCGCTGC‑GCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGAt > 1:1112921/1‑109 (MQ=255) aGGGGCCGCTGC‑GCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCg > 1:658078/1‑150 (MQ=255) aGGGGCCGC‑TGCGCGGCCCCTTGCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCg > 1:405099/1‑150 (MQ=255) gggCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTg < 2:1117479/102‑1 (MQ=255) gggCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTg > 1:1117479/1‑102 (MQ=255) gc‑gcAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCt > 2:1218084/1‑69 (MQ=255) gc‑gcAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCt < 1:1218084/69‑1 (MQ=255) c‑gcAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGCCTGgcgc > 2:745806/1‑150 (MQ=255) | TGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGCCTGGCGC > NC_002947/5777310‑5777605 |
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: ‑ |