Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 699,725 | A→G | 49.0% | intergenic (+115/‑129) | PP_t06 → / → PP_23SD | tRNA‑Ala/23S ribosomal RNA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 699,725 | 0 | A | G | 49.0% | ‑5.0 / 54.6 | 43 | intergenic (+115/‑129) | PP_t06/PP_23SD | tRNA‑Ala/23S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base A (10/12); new base G (14/7); total (24/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.23e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01 |
TCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATT > NC_002947/699581‑699804 | tCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGaa < 1:29750/150‑1 (MQ=38) tCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGACGGACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGACAGATATAGACTGaa < 2:193652/150‑1 (MQ=255) aGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAaca > 2:802815/1‑150 (MQ=255) aGCGGTTCGATCCCGCTTGGCTCCACCACTTTCACGGTAAGCTGATTTTTGACTTGACTTGGGAATTGACAGACAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGCGATAGATATAGGCTGAaca > 1:328566/1‑150 (MQ=2) ggTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCa < 2:244232/150‑1 (MQ=38) ttCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGt > 1:1184982/1‑150 (MQ=255) cGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGttt > 1:1040101/1‑150 (MQ=255) aTCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCATTTTCa > 1:593705/1‑150 (MQ=255) aTCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTAAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACAACAATTTCa > 2:469036/1‑150 (MQ=255) tCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAc < 1:1576/150‑1 (MQ=38) cGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActgc > 1:542532/1‑150 (MQ=39) cGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActgc > 1:181858/1‑150 (MQ=39) ttGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtgg < 1:926522/150‑1 (MQ=38) tGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggt < 2:1080491/150‑1 (MQ=39) cTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGa > 2:121850/1‑150 (MQ=38) caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGCGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCa < 2:328566/150‑1 (MQ=255) accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAg > 1:545352/1‑150 (MQ=38) ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCActgct > 1:702476/1‑139 (MQ=255) ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCActgct < 2:702476/139‑1 (MQ=255) ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAgg > 2:610373/1‑150 (MQ=255) ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg < 2:122577/150‑1 (MQ=39) ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCGTTTG‑CAGGCCAAGGGTTTAAAAATTCGGGTATGTGATAGATATAGACTGAACACCAGTTTTCATTCTGGTGGATTAgg > 1:677195/1‑150 (MQ=11) tttCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa < 1:87361/150‑1 (MQ=255) gCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg > 1:4983/1‑142 (MQ=39) gCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg < 2:4983/142‑1 (MQ=39) gTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCGGGTGGATCAGGCTAAGGTaaaa > 2:24064/1‑150 (MQ=37) tACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAt > 2:306847/1‑150 (MQ=255) gTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa < 1:847863/150‑1 (MQ=38) tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt > 2:996105/1‑150 (MQ=35) aGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCg > 2:914826/1‑150 (MQ=35) gAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGa > 2:593620/1‑150 (MQ=35) aCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa < 1:1073774/150‑1 (MQ=35) cTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAag < 1:914826/150‑1 (MQ=35) ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGCTCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga > 1:30829/1‑150 (MQ=25) tAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAagag < 1:802815/150‑1 (MQ=25) tAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAATGCTGGTGGATCAGGCTAAGGTAACATTTTTATTGTCTGCGCGAAagag < 2:677195/150‑1 (MQ=11) aaaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa < 1:917484/131‑1 (MQ=37) aaaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa > 2:917484/1‑131 (MQ=37) aaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc < 2:542532/150‑1 (MQ=35) aaCATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTAc > 1:14579/1‑150 (MQ=25) tCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGt > 1:967704/1‑83 (MQ=37) tCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGt < 2:967704/83‑1 (MQ=255) tCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAtt > 2:247641/1‑150 (MQ=18) | TCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGT‑TCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATT > NC_002947/699581‑699804 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |