Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,423,383 | A→G | 23.5% | intergenic (+189/+40) | yhdZ → / ← rrfF | putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,423,383 | 0 | A | G | 23.5% | 137.7 / 40.1 | 68 | intergenic (+189/+40) | yhdZ/rrfF | putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon |
Reads supporting (aligned to +/- strand): ref base A (26/26); new base G (9/7); total (35/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.78e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01 |
GATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACT > NC_000913/3423249‑3423531 | aatttgtctggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCtt < 2:399315‑M1/22‑1 (MQ=255) atttgtctggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTc < 4:457692‑M1/23‑1 (MQ=255) tttgtctggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTcc < 1:697432‑M1/24‑1 (MQ=255) ctggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTg < 3:359411‑M1/29‑1 (MQ=255) ggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAg > 3:457692‑M1/119‑149 (MQ=255) gcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGc < 4:641411‑M1/32‑1 (MQ=255) ggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCt > 4:359411‑M1/116‑149 (MQ=255) ggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCt > 3:1331794‑M1/116‑149 (MQ=255) ggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt > 1:399315‑M1/116‑149 (MQ=255) gcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCtt > 1:456454‑M1/115‑149 (MQ=255) gcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt > 2:697432‑M1/115‑149 (MQ=255) tagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGt > 3:641411‑M1/111‑149 (MQ=255) tagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGc > 3:1378499‑M1/111‑142 (MQ=255) tagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGc < 4:1378499‑M1/32‑1 (MQ=255) tcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt < 2:752766‑M1/37‑1 (MQ=255) tcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt > 1:752766‑M1/100‑136 (MQ=255) acctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGc < 4:123090‑M1/14‑1 (MQ=255) acctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGc > 3:123090‑M1/96‑109 (MQ=255) ctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCt < 3:327972‑M1/34‑1 (MQ=255) ctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCt > 4:327972‑M1/94‑127 (MQ=255) cccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt > 4:1428409‑M1/89‑122 (MQ=255) cccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt < 3:1428409‑M1/34‑1 (MQ=255) ccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGc > 4:429407‑M1/88‑119 (MQ=255) ccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGc < 3:429407‑M1/32‑1 (MQ=255) catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt > 1:234009‑M1/87‑120 (MQ=255) catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt < 2:234009‑M1/34‑1 (MQ=255) cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc < 3:1029105/106‑1 (MQ=255) cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc > 4:1029105/1‑106 (MQ=40) gaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGcc > 4:168154‑M1/81‑113 (MQ=255) gaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGcc < 3:168154‑M1/33‑1 (MQ=255) tcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAGAGGCCCAGTCTTTCGACTGAGc < 4:178072‑M1/32‑1 (MQ=255) tcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAGAGGCCCAGTCTTTCGACTGAGc > 3:178072‑M1/77‑108 (MQ=255) agaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt < 2:747899‑M1/37‑1 (MQ=255) agaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt > 1:747899‑M1/75‑111 (MQ=255) ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt < 2:1431621/113‑1 (MQ=255) ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt > 1:1431621/1‑113 (MQ=255) ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc < 3:977342/92‑1 (MQ=255) ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc > 4:977342/1‑92 (MQ=255) aaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGc > 2:1116491‑M1/68‑99 (MQ=255) aaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGc < 1:1116491‑M1/32‑1 (MQ=255) acgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCtt > 3:1410685‑M1/66‑100 (MQ=255) acgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCtt < 4:1410685‑M1/35‑1 (MQ=255) gTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt > 4:791340/1‑97 (MQ=255) gTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt < 3:791340/97‑1 (MQ=255) gcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc > 1:865828‑M1/58‑90 (MQ=255) gcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc < 2:865828‑M1/33‑1 (MQ=255) aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc < 3:1036186/83‑1 (MQ=255) aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc > 4:1036186/1‑83 (MQ=255) aTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc > 4:260999/1‑82 (MQ=255) aTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc < 3:260999/82‑1 (MQ=255) ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAt > 2:992421/1‑89 (MQ=255) ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAt < 1:992421/89‑1 (MQ=255) ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc < 2:253651/107‑1 (MQ=255) ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc > 1:253651/1‑107 (MQ=255) gcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACg > 1:1442283/1‑149 (MQ=255) cccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt < 1:157606‑M1/35‑1 (MQ=255) cccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt > 2:157606‑M1/34‑68 (MQ=255) tCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa < 2:1442283/149‑1 (MQ=255) agagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt < 2:1372437‑M1/37‑1 (MQ=255) agagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt > 1:1372437‑M1/26‑62 (MQ=255) tCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc > 2:734708/1‑44 (MQ=40) tCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc < 1:734708/44‑1 (MQ=40) acaacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg > 4:1346152/1‑144 (MQ=40) acaacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg < 3:1346152/144‑1 (MQ=40) tAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg < 2:539010/76‑1 (MQ=255) tAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg > 1:539010/1‑76 (MQ=255) aaaaCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg < 4:1303368/129‑1 (MQ=255) aaaaCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg > 3:1303368/1‑129 (MQ=255) aaaaGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACt > 1:239662/1‑149 (MQ=38) | GATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACT > NC_000913/3423249‑3423531 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |