Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,762,741 | A→T | K662I (AAA→ATA) | spoT → | bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,762,741 | 0 | A | T | 97.3% | 142.1 / ‑6.2 | 37 | K662I (AAA→ATA) | spoT | bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase |
Reads supporting (aligned to +/- strand): ref base A (0/0); major base T (13/23); minor base C (1/0); total (14/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.78e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCGCA > REL606/3762672‑3762815 | cgcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATgg > 1:252818/1‑76 (MQ=255) tGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcg > 1:370847/1‑76 (MQ=255) aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa < 2:916639/76‑1 (MQ=255) aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa < 2:846139/76‑1 (MQ=255) ccTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc < 1:874491/76‑1 (MQ=255) cTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa > 1:947965/1‑76 (MQ=255) tGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAg < 2:650218/76‑1 (MQ=255) cggcggCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgcgc < 1:48868/76‑1 (MQ=255) gcggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCt > 2:374570/1‑76 (MQ=255) cggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt < 2:1191060/76‑1 (MQ=255) cggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt < 2:1293119/76‑1 (MQ=255) ggcAATTAACACCACGACTTCGAATATTCAAAGTTTTAAGACGGAAGANATAGATGGTCGCGTCTACAGCGCCttt < 1:1168877/76‑1 (MQ=255) aaTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt < 2:1110931/76‑1 (MQ=255) aTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTc < 1:702831/76‑1 (MQ=255) ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg < 1:765496/76‑1 (MQ=255) tAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGt < 1:709582/76‑1 (MQ=255) acacCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCt < 1:695459/76‑1 (MQ=255) cacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAcc < 2:1115428/76‑1 (MQ=255) cGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGc > 2:420195/1‑76 (MQ=255) gACTTCGAATATTCAAAGGTTGAATACGGAAGGGACGGATGGTCGCGTCTACAGCGCCTTTATGCGTCTGACCGCt > 2:29019/1‑76 (MQ=255) ttCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGt > 2:1110088/1‑76 (MQ=255) tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg < 2:966041/76‑1 (MQ=255) tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg < 1:734623/76‑1 (MQ=255) tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg < 1:1026517/76‑1 (MQ=255) gAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAc > 1:485040/1‑76 (MQ=255) aaTATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGAcc > 2:100241/1‑76 (MQ=255) ttCAAAGTTTGAATACGGGAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtgt > 1:1265744/1‑76 (MQ=255) aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCa < 1:526393/76‑1 (MQ=255) aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCa < 2:629073/76‑1 (MQ=255) gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCt < 2:636933/76‑1 (MQ=255) gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCt < 1:476268/76‑1 (MQ=255) tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTg > 1:431965/1‑76 (MQ=255) tACGGAAGAGATAGATGGTCTCGTCTACAGCGACTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat > 2:107019/1‑76 (MQ=255) tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat < 2:379911/76‑1 (MQ=255) cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAatat < 2:174765/76‑1 (MQ=255) cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAatat > 2:958417/1‑76 (MQ=255) cGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAatat > 2:1176768/1‑76 (MQ=255) aaaGATCGTCNCAACTACAGCNGCTTTATTCGTCTGACCGCTCGTGACCGCGTGCATCTGGCGTATATCATCCGCa > 2:1319655/1‑76 (MQ=255) | CGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATGCGCA > REL606/3762672‑3762815 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |