Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,812,824 | A→G | 9.8% | intergenic (+75/‑35) | rplT → / → pheS | 50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,812,824 | 0 | A | G | 9.8% | 102.6 / 3.4 | 41 | intergenic (+75/‑35) | rplT/pheS | 50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha |
Reads supporting (aligned to +/- strand): ref base A (20/17); new base G (2/2); total (22/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01 |
TGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGCGA > NC_002947/2812676‑2812968 | tgtcTGATCTGGCCAGGCACGCACCAGCCGCCTATGCTGCGATCGTCGCGAACGCTACAGCCCTACTCGCCTAACTACCCGCGGCAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAg < 2:283002/147‑1 (MQ=255) aGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCTGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTAtt > 2:182717/1‑150 (MQ=255) gAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTNAACGTCNTCGATAGGGGAANAGCCTTCAAAGCTCTTCCCCTATTTTc < 1:504826/150‑1 (MQ=255) aaCGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCg > 1:174183/1‑150 (MQ=255) cGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGaa > 2:590535/1‑150 (MQ=255) gTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCt > 1:42364/1‑128 (MQ=255) gTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCt < 2:42364/128‑1 (MQ=255) tCGAGAACGCTAAAGCCGTTCTGGCTTAAGTCCCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCCAAGCTCTTCCCCTATTTTCGTATCTGGCGTGTGTACATGGAAAACCTg > 2:13791/1‑150 (MQ=255) gagaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGa > 2:21876/1‑150 (MQ=255) cTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTg > 2:266875/1‑150 (MQ=255) tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTg < 1:356508/118‑1 (MQ=255) tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTg > 2:356508/1‑118 (MQ=255) ttCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGtctc < 1:25222/146‑1 (MQ=255) ttCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGtctc > 2:25222/1‑146 (MQ=255) ccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTg > 2:488266/1‑150 (MQ=255) ccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTg > 2:141716/1‑150 (MQ=255) cgcgACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTa > 1:198347/1‑140 (MQ=255) cgcgACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTa < 2:198347/140‑1 (MQ=255) aaTCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAcgcg > 2:40305/1‑150 (MQ=255) cccGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTc < 1:335189/76‑1 (MQ=255) cccGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTc > 2:335189/1‑76 (MQ=255) cggTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATTCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACAt < 1:13791/150‑1 (MQ=255) tGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCaa < 1:141716/150‑1 (MQ=255) gcTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTATGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTg > 1:47110/1‑150 (MQ=255) gCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGaa < 1:488266/150‑1 (MQ=255) gggCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAg > 2:462869/1‑150 (MQ=255) cATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTaga < 1:239628/107‑1 (MQ=255) cATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTaga > 2:239628/1‑107 (MQ=255) cATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATc < 1:590535/150‑1 (MQ=255) gATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTAtct < 1:40305/150‑1 (MQ=255) aTAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTAtctc < 1:462869/150‑1 (MQ=255) gggAAGAGCCTTCAAAGCTGTTTCACTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCa < 2:47110/150‑1 (MQ=255) gggAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCa > 2:448474/1‑150 (MQ=255) ggAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCt > 2:346408/1‑60 (MQ=21) ggAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCt < 1:346408/60‑1 (MQ=21) ggAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTaa > 1:430051/1‑138 (MQ=255) ggAAGAGCTTTGAAGGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTaa < 2:430051/138‑1 (MQ=255) gagCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTCTACATGGAAATCCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCaagaag < 2:174183/150‑1 (MQ=255) ttCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGCGa > 2:354148/1‑150 (MQ=255) cAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGc < 1:153223/146‑1 (MQ=255) cAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGc > 2:153223/1‑146 (MQ=255) | TGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAGACATCAACGCCCTGGAACAGATCCGGGTTAATTATCTCGGCAAGAAGGGCGA > NC_002947/2812676‑2812968 |
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: ‑ |