Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,327,127 | C→T | 58.1% | intergenic (+97/‑196) | PP_16SE → / → PP_23SE | 16S ribosomal RNA/23S ribosomal RNA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,327,127 | 0 | C | T | 58.1% | ‑2.5 / 47.5 | 50 | intergenic (+97/‑196) | PP_16SE/PP_23SE | 16S ribosomal RNA/23S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base C (13/8); new base T (15/14); total (28/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.69e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAA > NC_002947/1326982‑1327235 | gTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCg > 2:537975/1‑150 (MQ=17) gTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTAATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAACTGCGCATTCGCATCGAATGTTGATTTCTGGCt > 1:402624/1‑150 (MQ=12) tAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCtt < 1:281207/150‑1 (MQ=17) ggttccccTACGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCtttt < 1:834129/144‑1 (MQ=255) ggggAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCCGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCtttt < 2:574062/150‑1 (MQ=17) ggggAACCTGCGGCGTGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGGATTGCTTGATTAATCATCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCtttt < 1:700254/150‑1 (MQ=11) gAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGATGATGTTTCCGATGCGTTA‑AACCGGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTc > 1:101342/1‑150 (MQ=11) aCCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTgg > 1:826412/1‑140 (MQ=17) aCCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTgg < 2:826412/140‑1 (MQ=17) ggCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTAt < 2:763419/117‑1 (MQ=18) ggCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGtt < 2:408511/150‑1 (MQ=32) ggATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGAt < 2:631591/142‑1 (MQ=17) ggATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGAt > 1:631591/1‑142 (MQ=17) gagcACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTTAGATCGTTCTTTa > 1:851262/4‑150 (MQ=12) gATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTa > 1:464877/1‑150 (MQ=32) cACCTCCTTACTCGACGACATCAGCCTGCTGATGAGCTCTCACACGAATTGCTTGATGCATTGTCGAAGACGTTTTTCCTG‑GTCAGATCCTGGAAATGAGAATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTaaaa < 2:86832/150‑1 (MQ=11) ctcctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAAtt > 2:426697/1‑150 (MQ=32) cctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCg < 1:426697/150‑1 (MQ=32) ctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCgg > 2:238635/1‑150 (MQ=17) ctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCa < 1:766095/129‑1 (MQ=32) ctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCa > 2:766095/1‑129 (MQ=32) ctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCtt < 2:244939/139‑1 (MQ=32) ctTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCtt > 1:244939/1‑139 (MQ=32) aaTCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGata > 1:834683/1‑150 (MQ=17) acgacATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTga > 2:580181/1‑150 (MQ=32) gCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAg < 1:2233/150‑1 (MQ=32) ccTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTg > 1:453751/1‑108 (MQ=17) ccTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTg < 2:453751/108‑1 (MQ=17) ccTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGa > 2:89139/1‑150 (MQ=32) gctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTg > 2:446468/1‑150 (MQ=17) tgaGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTg > 1:374868/1‑93 (MQ=21) tgaGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTg < 2:374868/93‑1 (MQ=21) aGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCa < 1:89139/150‑1 (MQ=32) gCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAg > 2:78252/1‑150 (MQ=32) tCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGtt > 1:734866/1‑150 (MQ=32) cccACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGttt > 1:587527/1‑150 (MQ=17) ccACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTc > 1:321682/1‑150 (MQ=32) cacaCGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCa < 1:537975/150‑1 (MQ=17) acacGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAc > 1:97131/1‑150 (MQ=17) acacGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGCATGTCGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAc > 1:530443/1‑150 (MQ=17) acGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActg > 1:513463/1‑150 (MQ=32) ttGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggt < 2:101342/150‑1 (MQ=17) ttGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGAt < 2:834683/150‑1 (MQ=17) gagtCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgat > 1:851561/4‑111 (MQ=16) gATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTgat < 2:851561/111‑1 (MQ=18) ttGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGaa > 2:398632/1‑114 (MQ=17) ttGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGAc < 1:398632/115‑1 (MQ=33) tGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg > 1:257159/1‑150 (MQ=17) tGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg < 2:97131/150‑1 (MQ=17) gAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTGGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGCTCAGGCTAAGGaaa > 1:359643/1‑150 (MQ=2) | GTCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGACATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTG‑GTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAA > NC_002947/1326982‑1327235 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |