Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,932,638:1 | +C | 100% | intergenic (‑105/+45) | rluA ← / ← minE | bifunctional tRNA pseudouridine(32) synthase TrmQ/ribosomal large subunit pseudouridine synthase RluA/cell division topological specificity factor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,932,638 | 1 | . | C | 97.8% | 147.5 / ‑3.0 | 45 | intergenic (‑105/+45) | rluA/minE | bifunctional tRNA pseudouridine(32) synthase TrmQ/ribosomal large subunit pseudouridine synthase RluA/cell division topological specificity factor |
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (22/22); total (23/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTG‑CCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT > NC_002947/1932493‑1932781 | gcggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACCAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTG‑ccc > 5:289612/1‑148 (MQ=255) ggcgTCTTGGTAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGAAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTG‑cccgc < 2:328315/149‑4 (MQ=255) gTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCTGGGGGCGAAACAACAAAGGATATGCGTCCCCTGAAGCACCGGGC‑TTG‑ccc > 1:136412/1‑148 (MQ=255) tCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAg < 2:87869/149‑1 (MQ=255) gAATAGGATGTGCACATCCGACAGGGGCATTGCAGGTCTCACTAAAATCGCAGTAACCGCTGGGGCGGCTTTGGGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAACGGATATGCGGCCCCTGTCGGAGCGGGC‑TTGCCCGCGAAGGGCTg < 4:225976/149‑1 (MQ=25) ggATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAACC‑‑‑‑‑CCGCTCCTACAAAGCATATGCCTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAgg < 1:247718/149‑1 (MQ=38) tCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCACCTAAAAAAGATATTCGTACCCTGTATGAGCCGGCGTTC‑CCGCGAAGGGCTGCAATGCAg > 2:332530/1‑149 (MQ=255) ggACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑ccc > 5:344367/1‑149 (MQ=255) aGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTTCAAAGGATATGCGTCCCCAGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAAATGTC‑Ga > 2:90750/1‑149 (MQ=37) cAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GCCTACGCGGAt > 2:275107/1‑149 (MQ=255) tctcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGc < 3:117735/149‑1 (MQ=255) tcACTAAAATAGCAGAAAAGGGTGGGGCGGTTTGGGGGCCGTGGGGGGGCAAGC‑‑‑‑‑CGGCTCTTACAAAGTATTTGCGCCCCCTGTTGGTGGGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAg < 4:186282/146‑1 (MQ=255) tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCGCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCCCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GGCTACGCGGCTCAg > 3:186282/1‑146 (MQ=37) ctaaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCGTTCGCGGGAAAGC‑‑‑‑‑CCGCTCATACAAAGTATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATcc < 6:344367/147‑1 (MQ=38) gAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCaa > 4:73288/1‑110 (MQ=255) gAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCaa < 3:73288/110‑1 (MQ=255) aaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTTTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCa < 4:57896/149‑1 (MQ=255) cAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGAAATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGt > 6:294002/1‑148 (MQ=255) gTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGAt > 3:185545/1‑149 (MQ=255) gTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTG‑ccc < 7:65463/86‑2 (MQ=255) gTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTG‑ccc > 8:65463/1‑85 (MQ=255) tAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑cc > 7:35005/1‑106 (MQ=255) tAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTAAAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑cc < 8:35005/106‑1 (MQ=255) cggtGGGGCCCGCTTGCGGCCCTTGGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCGGTAGGCGGGGGC‑TTGCCCGCGAATGGCTGCAAGGCAG‑CCCCAAATTGTCTGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTc < 1:275107/146‑1 (MQ=255) ccgtGGGGCCGCTTTGCGGCCCTTGGCGGGCAAGC‑‑‑‑‑CCGCTCGTTCAAAGGTTATCCGTCGCCTGTAGGAGCGGGG‑TTGCCCGCGAAGGGCTGCAAGGCAGCCCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTc < 1:90750/146‑1 (MQ=255) ggggCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATACGGCAGGGTGATGTTCAGCTc > 6:322331/1‑149 (MQ=255) gggCCGCTTTTCGGCCCTTCGCTGGCAACC‑‑‑‑‑CCGCTCCTTCAAAGCATCTCCGGCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCt < 6:345410/147‑1 (MQ=37) gCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCGGCTACAAAGGATATGGGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGACGGGCTGCAAGGCAG‑CCCCAAAATGTc > 7:102026/1‑105 (MQ=25) tGCGGCCGTTCGCTGGAAAGC‑‑‑‑‑CCGTTCCTACAAAGGATATGCGTCCCCTGTAGGTGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGa < 2:136412/148‑1 (MQ=39) gCGGCCAGTAGCTGGAAAGT‑‑‑‑‑CCTCTCCTACAAAGGATTTGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGag < 6:289612/148‑1 (MQ=37) ttCGCGGGCAAGC‑‑‑‑‑CCGCTCCTAAAAACGAAATGCGTCCCCAGTAGGAGCCGGC‑TTGCCCGCCAAGGGCTGCAACGCAG‑CACCAAATTGAC‑GACTACGCGGATCAGCGCTCCGGGAGGGTGATGTTGAGATCCAGAATCGagcagcgtca > 4:307847/1‑145 (MQ=255) gagcgGGCTTGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAg > 4:328070/3‑115 (MQ=38) gagcgGGCTTGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAg < 3:328070/113‑1 (MQ=38) gggCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATc < 5:403003/98‑1 (MQ=255) gggCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATc > 6:403003/1‑98 (MQ=255) aaGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc < 5:294002/149‑1 (MQ=255) tCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGt < 5:80851/107‑1 (MQ=255) tCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGt > 6:80851/1‑107 (MQ=255) tACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat > 5:235947/1‑149 (MQ=255) aGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGa < 1:145133/141‑1 (MQ=255) aGGATATGCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGa > 2:145133/1‑141 (MQ=255) gCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATcc < 2:352124/79‑1 (MQ=255) gCGTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATcc > 1:352124/1‑79 (MQ=255) gTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg > 4:68235/1‑149 (MQ=255) gTCCCCTGTAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTg > 5:6940/1‑149 (MQ=255) tAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCa > 6:384505/1‑85 (MQ=255) tAGGAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCa < 5:384505/85‑1 (MQ=255) ggAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTca < 4:163303/149‑1 (MQ=255) gAGCGGGC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTca > 4:122430/1‑148 (MQ=255) gC‑TTGCCCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAttt > 2:385475/1‑149 (MQ=255) | GCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGC‑‑‑‑‑CCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGC‑TTG‑CCGCGAAGGGCTGCAAGGCAG‑CCCCAAATTGTC‑GACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT > NC_002947/1932493‑1932781 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |