Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,327,127 | C→T | 54.6% | 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 | 54.6% | 1.6 / 65.0 | 53 | intergenic (+97/‑196) | PP_16SE/PP_23SE | 16S ribosomal RNA/23S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base C (14/10); new base T (19/10); total (33/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.76e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01 |
TCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAA > NC_002947/1326983‑1327230 | tCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTC‑‑TCACCTACCACCAAACAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGa > 2:271594/1‑150 (MQ=11) tAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATg > 2:1099143/1‑150 (MQ=17) aCAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGtt < 2:937035/150‑1 (MQ=32) aaGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGa > 1:768245/1‑150 (MQ=17) tAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGCAAATGAGCATTCGCACCGAATGTTGATTTc > 2:417123/1‑150 (MQ=17) ggggAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCtttt < 1:1004121/150‑1 (MQ=32) gggAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTAGCTTTTg < 2:1128214/150‑1 (MQ=17) gAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTTTTTTTTTTTACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGtt > 2:77118/1‑149 (MQ=11) gCTGGATCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTc > 2:636240/1‑150 (MQ=17) gATCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTa < 1:1099143/150‑1 (MQ=17) aTCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTaa > 2:113276/1‑150 (MQ=17) aTCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTg < 1:1053113/127‑1 (MQ=33) aTCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTg > 2:1053113/1‑127 (MQ=32) tCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCTTTCTTAGAACCTGGAAATGAGCATTCGTATCGAATGCTGATTTCTGGCTGTCGGCAGAGCTTTCGTTaaa > 1:801302/1‑150 (MQ=11) tCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTaaa > 2:77614/1‑150 (MQ=17) cACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTaaaa > 2:219410/1‑150 (MQ=32) aCCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTaaaaa > 2:633302/1‑150 (MQ=32) tcctTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTc > 2:158319/1‑150 (MQ=17) tcctTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTc > 2:209779/1‑150 (MQ=32) cctTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTTTTTTTTTTTTCTTGTAAATGAGCATTCGTATATAATGTTGATTTCTGGCTTTTGTCAGATCTTTCTTTAAAAATTct > 2:519050/1‑149 (MQ=11) ctTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCg‑g < 2:947456/150‑1 (MQ=17) aTCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCGTGATTCATTGTCGAAGACGTTTTTCCTGGGCAGAGCCTGGAAAGGCGCATTCGCGTCGCATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTACAAATTCg‑ggtct > 1:406025/1‑146 (MQ=11) tCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATAtg > 1:686321/1‑150 (MQ=17) tCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATAtg > 2:756751/1‑150 (MQ=17) cgacgaC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCTTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATAtgt < 1:413347/150‑1 (MQ=17) gac‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTgata < 1:113276/150‑1 (MQ=17) gac‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTgata < 1:235175/150‑1 (MQ=17) c‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTgataga < 2:152911/150‑1 (MQ=17) c‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTgataga > 1:99886/1‑150 (MQ=32) aTCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCGTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTgatagat > 1:341311/1‑150 (MQ=17) gCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAg > 1:826701/1‑150 (MQ=17) ccTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGACGCCGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGa > 1:1041697/1‑150 (MQ=12) ctgctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGACTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGaa > 2:411187/1‑149 (MQ=17) ctgctgATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGAc < 2:706409/150‑1 (MQ=32) tgatgaGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGaa < 2:763008/150‑1 (MQ=32) cTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGt > 1:1124046/1‑150 (MQ=32) tCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAAt < 2:325230/77‑1 (MQ=35) tCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAAt > 1:325230/1‑77 (MQ=35) ccACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTc < 1:778659/115‑1 (MQ=17) ccACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTc > 2:778659/1‑115 (MQ=17) cacaCGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCa < 1:625622/150‑1 (MQ=17) acacGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCAc > 2:13669/1‑150 (MQ=17) cacGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCAct < 1:271594/150‑1 (MQ=17) ttGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggt < 1:807562/150‑1 (MQ=32) gCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtgg < 1:848414/150‑1 (MQ=32) ttGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGAt > 1:349029/1‑150 (MQ=32) ttCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg > 2:707201/1‑150 (MQ=17) ttCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg < 2:99886/150‑1 (MQ=32) ttCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg > 2:1124177/1‑150 (MQ=32) cATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGCCCAATACTGCGCAAACCGTGCGCGATGA‑AGATGTAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCt > 1:407126/1‑150 (MQ=2) ttGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa > 2:55726/1‑150 (MQ=17) tCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCAc > 1:984077/1‑128 (MQ=17) tCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGTATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCAc < 2:984077/128‑1 (MQ=17) | TCGTAACAAGGTAGCCGTAGGGGAACCTGCGGCTGGATCACCTCCTTAATCGACGAC‑‑ATCAGCCTGCTGATGAGCTCCCACACGAATTGCTTGATTCATTGTCGAAGACGTTTTTCCTGGTCAGAACCTGGAAATGAGCATTCGCATCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCG‑GATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAA > NC_002947/1326983‑1327230 |
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: ‑ |