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 | 98.2% | 207.7 / NA | 55 | K662I (AAA→ATA) | spoT | bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base T (31/23); total (32/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATG > REL606/3762668‑3762811 | ggTGCGCGGGCAAACCTGACGGAGGAAATTAACACCGCGACTTTGAATATTCAAAGTTTGAATACGGAAGAGATAg < 1:1295264/76‑1 (MQ=255) gTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGa < 1:648597/76‑1 (MQ=255) tGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGagatagat < 1:398833/76‑1 (MQ=255) ggcgcTTGCAAACCTGACGGCGGCAATTAACCCCACGACTTCGAATATTCAAAGTTTGAATACGGAAGagatagat < 1:398832/75‑1 (MQ=255) gcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGt < 2:414558/76‑1 (MQ=255) cTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTc > 1:768990/1‑76 (MQ=255) gCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTcgcg < 2:1499464/76‑1 (MQ=255) cAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGt > 1:1241850/1‑76 (MQ=255) cAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGt > 1:1241827/1‑76 (MQ=255) aaCCTGACGGCGGCTATTAACACCGCGACTTNGAATATTCGAAGTTNGAATANGGAAGAGATAGATGGTCGCGTCt < 2:203502/76‑1 (MQ=255) aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa < 1:1338242/76‑1 (MQ=255) aCCTGACGGCGACAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa < 1:568918/76‑1 (MQ=255) ccTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTAc < 1:557041/76‑1 (MQ=255) cTGACGGCTGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa < 1:1370494/76‑1 (MQ=255) tGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAg < 1:355990/76‑1 (MQ=255) gACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgc > 1:458763/1‑76 (MQ=255) aCGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg < 2:744898/76‑1 (MQ=255) cggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt < 2:868450/76‑1 (MQ=255) cggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCtt < 1:1124265/76‑1 (MQ=255) gcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTa > 1:519671/1‑76 (MQ=255) cAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAt > 2:1164935/1‑76 (MQ=255) aaTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt < 1:699829/76‑1 (MQ=255) aTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTc > 2:960921/1‑76 (MQ=255) ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg < 2:439135/76‑1 (MQ=255) ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg > 1:699042/1‑76 (MQ=255) ccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAc > 2:118046/1‑76 (MQ=255) aCTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTc > 1:373485/1‑76 (MQ=255) cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCg > 2:1488639/1‑76 (MQ=255) cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCg > 2:1070762/1‑76 (MQ=255) ttCGAATATTCAAAGTTTGNATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTCTTCGTCTGACCGCTCGt > 1:1432706/1‑76 (MQ=255) tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg < 2:97713/76‑1 (MQ=255) cGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGGGa > 1:1008135/1‑76 (MQ=255) cGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGATCGTGa > 1:1390686/1‑76 (MQ=255) tatTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgt > 1:544469/1‑76 (MQ=255) nttCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGANCGCTCGTGACCgtg > 2:219758/2‑76 (MQ=255) atTCAAAGTTTGAATACGGAAGAGATAGTTGGTCCCATCTACAGCGCNTTAATTCGTCTGACCGCTCGTGACCgtg > 1:142200/1‑76 (MQ=255) atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtg > 1:574445/1‑76 (MQ=255) ttCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtgt > 1:524912/1‑76 (MQ=255) tCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGtgtg < 1:559312/76‑1 (MQ=255) aaaGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCa > 1:495958/1‑76 (MQ=255) aGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATc > 1:285751/1‑76 (MQ=255) ttGAATACGGAAGAGATAGATGGTCGCGTCTACNGCGCCTTTATTCGTCTGACCGCNCGTGACCGAGNGCATCTgg > 2:1278581/1‑76 (MQ=255) ttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTgg > 1:835910/1‑76 (MQ=255) tGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGc < 1:983826/76‑1 (MQ=255) gAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCg < 1:735331/76‑1 (MQ=255) gAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCg < 2:771604/76‑1 (MQ=255) gAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCg < 1:594855/76‑1 (MQ=255) aaTACGGAAGAGATAGATGGTCGCGTCTACTGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa > 2:66729/1‑76 (MQ=255) aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa > 2:581060/1‑76 (MQ=255) aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa > 2:522191/1‑76 (MQ=255) aaTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGa > 2:289411/1‑76 (MQ=255) aTACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGaa > 2:811593/1‑76 (MQ=255) tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat < 1:127731/76‑1 (MQ=255) tACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAat > 2:1338144/1‑76 (MQ=255) gAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCa < 2:731898/76‑1 (MQ=255) aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACGNCTNGTGACCGTGTGCATNTGGCGAATATCAt > 2:1384508/1‑76 (MQ=255) aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt > 2:524356/1‑76 (MQ=255) agagATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATg > 2:849387/1‑76 (MQ=255) | GGTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATG > REL606/3762668‑3762811 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |