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 | 96.0% | 189.6 / NA | 50 | 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 (20/28); minor base G (2/0); total (22/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.89e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATG > REL606/3762669‑3762811 | gTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGa < 1:852211/76‑1 (MQ=255) gTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGa < 1:284911/76‑1 (MQ=255) tGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGagatagat < 1:582531/76‑1 (MQ=255) tGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGagatagat > 1:291463/1‑76 (MQ=255) tGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGagatagat < 1:1275988/76‑1 (MQ=255) gcgcTGGCAAACCTGACGGCGGCAATTAGCACCACGTATTCGAATATNCACAGTTTGAATACGGAAGAGATAGATg > 1:137307/1‑76 (MQ=255) gcgcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATg < 2:1134096/76‑1 (MQ=255) gcgcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATg < 1:937137/76‑1 (MQ=255) gcTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGt > 2:1427413/1‑76 (MQ=255) aaCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATANATGGTCGCGTCt < 1:13619/76‑1 (MQ=255) aCCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTa < 1:1415442/76‑1 (MQ=255) cTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACa < 2:287787/76‑1 (MQ=255) aCGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAgcg < 1:1278647/76‑1 (MQ=255) ggcggcAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAGAGATGGTCGCGTCTACAGCGcc > 2:51682/1‑76 (MQ=255) aaTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCNCCTTTAtt > 1:1259977/1‑76 (MQ=255) aaTTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTAtt < 1:1326015/76‑1 (MQ=255) ttAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCg > 2:773320/1‑76 (MQ=255) tAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGt > 1:224423/1‑76 (MQ=255) agcaccacGACTTCGAATACTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTc > 1:931286/3‑76 (MQ=255) aaCACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATCCGTc < 1:82496/76‑1 (MQ=255) aaCACCACGACTTAGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTc > 1:107452/1‑76 (MQ=255) acacCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCt < 1:295317/76‑1 (MQ=255) caccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTg < 1:264946/76‑1 (MQ=255) cccacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGa < 1:875999/75‑1 (MQ=255) accacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGa < 1:336664/76‑1 (MQ=255) cacGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGAcc < 1:714298/76‑1 (MQ=255) acGACTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCg > 1:227932/1‑76 (MQ=255) acGACTTCGAATATTCAAAGGTCGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCg > 1:8065/1‑76 (MQ=255) cTTCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCg > 1:1125450/1‑76 (MQ=255) ttCGAATATTCAAAGTTTGAATANGGAAGNGATAGCTGGTCGCGNCTCCAGAGCCTTTATTAGTCTGACCGCTCGt > 2:191707/1‑76 (MQ=255) tCGAATATTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTg < 1:449286/76‑1 (MQ=255) atatTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCg > 1:547282/1‑76 (MQ=255) natTCAAAGTTTGAATACGCAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGNCCGCTCGTGACCgt > 2:198700/2‑76 (MQ=255) atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtg < 2:80520/76‑1 (MQ=255) atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtg > 1:661077/1‑76 (MQ=255) atTCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATCCGTCTGTCCGCTCGTGACCgtg > 2:43171/1‑76 (MQ=255) ttCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtgt < 1:491077/76‑1 (MQ=255) ttCAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCgtgt < 1:291572/76‑1 (MQ=255) cAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGc > 1:463713/1‑76 (MQ=255) cAAAGTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGc < 1:263765/76‑1 (MQ=255) aaaGTTTGAATACGNAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGcn < 2:204219/76‑2 (MQ=255) gTTTGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCt > 2:1093317/1‑76 (MQ=255) tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCCTCTg < 2:1429136/76‑1 (MQ=255) tttGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTg < 2:1286782/76‑1 (MQ=255) tGAATACGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGc > 1:367655/1‑76 (MQ=255) gAATACGGAAGAGATAGATGGTCGCGTCTCCAGCGCCTTTATTCGTCTGACCGCTCGNGACCGTGTGCATCTGGCg < 1:1376187/76‑1 (MQ=255) gAATACCGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCg < 1:1428677/76‑1 (MQ=255) aCGGAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAata < 1:1325666/76‑1 (MQ=255) ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATc > 2:1254070/1‑76 (MQ=255) ggAAGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATc < 1:406750/76‑1 (MQ=255) aaGAGATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt > 2:1063594/1‑76 (MQ=255) aaGAGAGAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAt > 1:89932/1‑76 (MQ=255) agagATAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCAtn < 1:617613/76‑2 (MQ=255) agatagatGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGc < 2:1122963/57‑1 (MQ=255) agatagatGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGc > 1:1122963/1‑57 (MQ=255) | GTGCGCTGGCAAACCTGACGGCGGCAATTAACACCACGACTTCGAATATTCAAAGTTTGAATACGGAAGAGAAAGATGGTCGCGTCTACAGCGCCTTTATTCGTCTGACCGCTCGTGACCGTGTGCATCTGGCGAATATCATG > REL606/3762669‑3762811 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |