Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 529,678 | C→A | 13.1% | intergenic (+14/‑106) | PP_23SC → / → PP_5SC | 23S ribosomal RNA/5S ribosomal RNA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 529,678 | 0 | C | A | 13.1% | 157.9 / 15.0 | 61 | intergenic (+14/‑106) | PP_23SC/PP_5SC | 23S ribosomal RNA/5S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base C (24/29); new base A (6/2); total (30/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.47e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-02 |
CTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATC > NC_002947/529575‑529824 | cTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGc < 8:177327/149‑1 (MQ=33) tACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGc > 1:254652/1‑148 (MQ=33) cgacgTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACAAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGAAGAAAGAACCGAAAATTCGTAAGGCca > 6:137253/1‑148 (MQ=255) ctAGCTTGGGTGTGTCAGCGCTGTGAGGCGTTGAGCTACCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATc < 6:305347/147‑1 (MQ=255) tAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGCGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCg > 4:264673/1‑148 (MQ=255) aGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCa > 3:21511/1‑149 (MQ=255) gTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCata > 8:279208/1‑149 (MQ=255) tgtgtAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGaa < 3:71689/149‑1 (MQ=255) aGCGCTTTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACTCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCt < 4:198346/149‑1 (MQ=255) aGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCt > 7:20057/1‑149 (MQ=255) gcgcTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTg < 7:279208/149‑1 (MQ=255) cgcTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTgg < 2:254652/149‑1 (MQ=255) gAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGt < 8:20057/148‑1 (MQ=255) gCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCa > 3:169564/1‑148 (MQ=255) gTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCt < 5:85168/149‑1 (MQ=255) gTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCt > 4:18542/1‑149 (MQ=255) gCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACa > 4:142569/1‑149 (MQ=255) aaCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGc < 3:266182/113‑1 (MQ=255) aaCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGc > 4:266182/1‑113 (MQ=255) aaCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTc > 3:228575/1‑149 (MQ=255) aCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCt < 7:131952/149‑1 (MQ=255) aCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCt > 3:127728/1‑149 (MQ=255) aCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCt > 8:131952/1‑149 (MQ=255) ccATTACTTATTGCCCGTGAGGGTTTACGATATAACACCCAAGCAATCTGCTACGCAGATTGTGTTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGACACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTg < 4:21503/149‑1 (MQ=18) cAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTgg < 4:228575/149‑1 (MQ=255) cccGTCACGCGTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTc < 8:60004/149‑1 (MQ=255) ccGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCt < 3:82311/149‑1 (MQ=255) ccGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCt > 5:75158/1‑149 (MQ=255) gTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGCGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTg < 3:264673/149‑1 (MQ=255) tGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTga > 8:150004/1‑149 (MQ=255) gCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTgacgac < 6:275810/149‑1 (MQ=255) gCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTgacgac < 3:52159/149‑1 (MQ=255) gCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAACTCCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTgacga < 4:241951/148‑1 (MQ=255) ttGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTg < 7:74219/117‑1 (MQ=255) ttGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTg > 8:74219/1‑117 (MQ=255) ttGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCa > 7:202029/1‑149 (MQ=255) aCCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATag < 8:202029/149‑1 (MQ=255) aCCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATa < 3:142569/148‑1 (MQ=255) atatAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGc < 2:240064/149‑1 (MQ=255) taACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACATAATTTCTTGACGACCATAGAGCAtt < 4:127728/149‑1 (MQ=255) aaCACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTcc > 1:17519/1‑102 (MQ=255) aaCACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTcc < 2:17519/102‑1 (MQ=255) aaCACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTg > 4:107313/1‑149 (MQ=255) acCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGaa > 2:101443/1‑149 (MQ=255) cccAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAac < 3:18542/149‑1 (MQ=255) ccAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACAtt < 1:287461/108‑1 (MQ=255) ccAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACAtt > 2:287461/1‑108 (MQ=255) caccacAATCTGCGTAGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAacc > 5:109619/6‑149 (MQ=255) accacAATCTGCGTAGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATc < 1:243973/95‑1 (MQ=25) accacAATCTGCGTAGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATc > 2:243973/5‑99 (MQ=25) accacAATCTGCGTAGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTAGGACAACAGAATTTCTTGACGACCATAGAGCATTGGAacca > 1:256658/5‑149 (MQ=18) ccacAATCTGCGTAGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAaccac > 6:72296/4‑149 (MQ=255) aaGCAATCTGCTTCGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAAATCGTAAGGCCACAAATATCGCATATACGAATTCGATGGGCTGTCGATATGGACATTCTGGCAACAGAATTTCTTGACGACCATAGATCATTGGAaccac > 5:161029/1‑149 (MQ=18) cacAATCTGCGTAGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGc < 6:233550/87‑1 (MQ=25) cacAATCTGCGTAGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGc > 5:233550/3‑89 (MQ=25) aaTCTGCGTAGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAACCCCCTGa > 5:40480/1‑149 (MQ=255) aTCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGat < 6:75158/149‑1 (MQ=255) cTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGatcccat > 3:133267/1‑149 (MQ=255) tACGGAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGATCATTGGAACCACCTGatcccatc < 4:21511/149‑1 (MQ=255) tACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCAt > 6:142556/1‑122 (MQ=255) tACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCAt < 5:142556/122‑1 (MQ=255) | CTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGTGGTGGTGAAGACGAAAGACCCGAAAATTCGTAAGGCCACAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATC > NC_002947/529575‑529824 |
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: ‑ |