| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 181,818 | G→A | 10.3% | intergenic (+32/‑74) | PP_r01 → / → PP_5SB | 5S ribosomal RNA/5S ribosomal RNA |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 181,818 | 0 | G | A | 10.3% | 150.9 / 6.9 | 58 | intergenic (+32/‑74) | PP_r01/PP_5SB | 5S ribosomal RNA/5S ribosomal RNA |
| Reads supporting (aligned to +/- strand): ref base G (25/27); new base A (4/2); total (29/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.70e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01 | |||||||||||
TTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCG > NC_002947/181670‑181960 | ttGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAgg > 4:299939/1‑149 (MQ=32) gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑ggg < 4:283735/149‑1 (MQ=32) gacCATAGAGCATTGGAAACACCTGATCCAATCCCGAACTCAGTAGTGAAACGATGCATAGCCGATGGTAGTGTGGGGTTTCACAATGTTAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑ggg > 2:154165/1‑149 (MQ=11) ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGt > 4:64722/1‑149 (MQ=32) ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGt < 6:258828/149‑1 (MQ=32) agCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTc > 2:138150/1‑149 (MQ=32) cATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCtt < 3:64722/149‑1 (MQ=32) aTTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCttt < 6:288437/149‑1 (MQ=32) ccaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTACCCATGTGAGAGTAGGTCATCGTAAAGATTCATTTCGCAAAACCCCTATCTGCACATGAAGGTA‑GGGGTTTTGTCTTTAAGTAg > 8:271538/1‑148 (MQ=12) ccTGATCCCATCCCGAAATCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACACCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAAc > 6:240991/1‑149 (MQ=25) atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTAc < 3:299939/148‑1 (MQ=38) cccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGt > 6:138482/1‑118 (MQ=32) cccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGt < 5:138482/118‑1 (MQ=32) ccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACaga < 4:162667/148‑1 (MQ=40) cGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑Gc > 7:208079/1‑149 (MQ=255) aatCATTAGTGAAACGATTCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACACCTATCTGCACATGCAGGTA‑GTGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCa < 7:271538/145‑1 (MQ=255) agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATA‑g < 2:252958/106‑1 (MQ=11) agtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTACCTGCATGTGCAGATA‑g > 1:252958/1‑106 (MQ=11) gtagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACag < 8:265089/133‑1 (MQ=40) tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTc > 4:182320/1‑149 (MQ=255) aaaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAgg > 3:49672/1‑149 (MQ=255) aTCGCCGATGGTAGTTGGGGGTTTCCCCATGTGAGAGTACGTCAGCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCACGCACGTCCTAGGGACGGGCCa < 8:287046/149‑1 (MQ=18) ggTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAg > 2:91383/1‑148 (MQ=255) ggTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAg > 6:20710/1‑148 (MQ=255) tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCAc < 3:154849/121‑1 (MQ=255) tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCAc > 4:154849/1‑121 (MQ=255) tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAAt < 1:138150/149‑1 (MQ=255) tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGATATTC‑TCAGGCACGTCCTAGGGACTGGCCAGCAAACAGAAtta > 3:4319/1‑147 (MQ=255) gggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTg < 4:49672/148‑1 (MQ=255) cccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCa > 8:129026/1‑148 (MQ=255) gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCCCGTccctagg < 3:182689/109‑6 (MQ=255) gagaGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCCCGTccctagg > 4:182689/1‑104 (MQ=255) agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTCTTCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAt < 1:91383/149‑1 (MQ=255) agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCAtt < 1:154165/149‑1 (MQ=255) tAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGaa > 6:333983/1‑149 (MQ=255) tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCGCACAGAATTTCTTGACGACCATaga > 1:40093/1‑136 (MQ=255) tCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCGCACAGAATTTCTTGACGACCATaga < 2:40093/136‑1 (MQ=255) aTCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAaccacc < 3:182320/148‑1 (MQ=255) aaGATTGATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcc < 5:240991/148‑1 (MQ=255) aaGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTga < 2:120881/118‑1 (MQ=255) aaGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTga > 1:120881/1‑118 (MQ=255) aTTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatc > 2:49864/1‑149 (MQ=255) ttaaagacAAAACCCCTACCTGCATGTGCAGATA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccatccc > 2:213133/8‑148 (MQ=37) tttCGCAAAACCCCTATCTGCACATGCAGGTAGGGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTgacga > 1:46020/1‑114 (MQ=255) tttCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTgacga < 2:46020/113‑1 (MQ=255) tttCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa > 4:321046/1‑149 (MQ=255) tttCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGaa > 7:61917/1‑149 (MQ=255) cAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag < 8:61917/149‑1 (MQ=255) aaaaCCCCTACCTGCATGTGCAGATA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag > 4:99470/1‑148 (MQ=255) aaaCCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTcag > 1:266289/1‑147 (MQ=255) ccccTACCTGCATGTGCAGATA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat > 4:227254/1‑128 (MQ=255) ccccTACCTGCATGTGCAGATA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat < 3:227254/128‑1 (MQ=255) tCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAAc < 5:20710/147‑1 (MQ=255) tCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAAc < 3:321046/147‑1 (MQ=255) cacaTGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATg < 2:266289/147‑1 (MQ=255) cacaTGCAGGTA‑GGGGTTTTGTATTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCAAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt < 4:4319/146‑1 (MQ=255) aTGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCg > 2:301889/1‑149 (MQ=255) aGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt > 2:147268/1‑139 (MQ=255) aGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGAt < 1:147268/139‑1 (MQ=255) | TTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCACATGCAGGTA‑GGGGTTTTGTCTTTAAGTAGGAACTACAGAGATTC‑GCAGGCACGTCCTAGGGACGGGCCAGCACACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCG > NC_002947/181670‑181960 |
| 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: ‑ |