| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 4,040,661 | G→A | 12.8% | intergenic (+25/+245) | rrfA → / ← mobB | 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,040,661 | 0 | G | A | 12.8% | 102.4 / 7.8 | 39 | intergenic (+25/+245) | rrfA/mobB | 5S ribosomal RNA of rrnA operon/molybdopterin‑guanine dinucleotide biosynthesis protein B |
| Reads supporting (aligned to +/- strand): ref base G (17/17); new base A (3/2); total (20/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 = 7.62e-01 | |||||||||||
TTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAA > NC_000913/4040515‑4040808 | tttGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGaa < 4:290199/149‑1 (MQ=34) tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTgg < 4:250739/149‑1 (MQ=35) ggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCt < 2:758588/149‑1 (MQ=35) cTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTc < 1:177440/149‑1 (MQ=35) cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAAc < 1:282965/149‑1 (MQ=35) aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAg > 4:330564/1‑149 (MQ=34) gcgcCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAT > 4:787175‑M1/1‑101 (MQ=255) gcgcCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAT < 3:787175‑M1/101‑1 (MQ=255) gTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATcc < 4:224419/124‑1 (MQ=35) gTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATcc > 3:224419/1‑124 (MQ=35) gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTg < 1:358129/128‑1 (MQ=32) gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTg > 2:358129/1‑128 (MQ=33) aaCTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGc < 3:330564/149‑1 (MQ=32) aTCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAgg > 4:449589/1‑133 (MQ=33) aTCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAgg < 1:161864/133‑1 (MQ=33) aTCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAgg < 3:449589/133‑1 (MQ=32) aTCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAgg > 2:161864/1‑133 (MQ=32) aaaTAAAACAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATttgatgtctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccac > 4:703260‑M1/1‑46 (MQ=255) aaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCa > 1:139240/1‑134 (MQ=33) aaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCa < 2:139240/134‑1 (MQ=33) aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcactt < 4:70420‑M1/111‑73 (MQ=255) aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcactt > 3:70420‑M1/1‑39 (MQ=255) gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctact > 3:330994‑M1/1‑37 (MQ=255) gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctga < 2:236212‑M1/113‑77 (MQ=255) gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctga > 1:236212‑M1/1‑37 (MQ=255) aaaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactg > 2:41971‑M1/1‑36 (MQ=255) aaGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgc > 3:391681‑M1/1‑35 (MQ=255) aGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgg < 3:91803‑M1/131‑98 (MQ=255) aGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgg > 4:91803‑M1/1‑34 (MQ=255) gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgccgc > 1:514169‑M1/1‑32 (MQ=255) gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatgaggagaccccacactaccatcggcgctacggcgtttcacttctgagttcggcatggggtcaggtgggaccaccgcgctactgccgc > 1:483371‑M1/1‑32 (MQ=255) cTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactacc > 2:124015‑M1/1‑31 (MQ=255) cTCAGTCGAAAGACTGGGCCTTTCGTTTTATttgatgcctggcagttccctactctcgcatggggagaccccacactacc < 1:124015‑M1/80‑50 (MQ=255) gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgccg > 1:361777/1‑73 (MQ=34) gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgccg < 4:393594/73‑1 (MQ=35) gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgccg > 3:393594/1‑73 (MQ=34) gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgccg < 2:361777/73‑1 (MQ=34) gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGa > 1:315601/1‑95 (MQ=33) gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGa < 2:315601/95‑1 (MQ=34) ggAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaaa > 1:436045/1‑149 (MQ=31) | TTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAA > NC_000913/4040515‑4040808 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |