Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 1,290,312 | +G | coding (277/1014 nt) | hnr → | response regulator of RpoS |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,290,311 | 1 | . | G | 92.6% | 99.9 / 2.7 | 27 | coding (276/1014 nt) | hnr | response regulator of RpoS |
Reads supporting (aligned to +/- strand): ref base . (2/0); new base G (11/14); total (13/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.22e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01 |
AAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCA‑GATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGC > REL606/1290237‑1290380 | aaaCTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCA‑g > 1:456436/1‑76 (MQ=255) aCTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCA‑gga > 1:41417/1‑74 (MQ=255) ctgGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATAtt > 2:926421/1‑76 (MQ=255) cATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCaaa < 1:19489/76‑1 (MQ=255) atatACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCA‑GATATTGCCAAAGc > 1:389734/1‑76 (MQ=255) taCGTNACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGt < 2:778059/76‑1 (MQ=255) aCGTANCAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGtt < 2:779866/76‑1 (MQ=255) aCGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACCTAAAATATGGCAGGATATTGCCAAAGCGtt < 1:946416/76‑1 (MQ=255) cGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTa < 2:681045/76‑1 (MQ=255) cGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTa > 1:266938/1‑76 (MQ=255) gTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTTCCAAAACGTTAc > 1:905734/1‑76 (MQ=255) aCAGAGGCGACCAGACACCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTc > 2:918585/1‑76 (MQ=255) agagGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTg > 2:530175/1‑76 (MQ=255) agGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTgag < 2:752442/76‑3 (MQ=255) ggCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGc < 1:202635/76‑1 (MQ=255) gCGACCAGACCCCAGTTCTGGTGGTATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCg < 1:662799/76‑1 (MQ=255) aCCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGAt < 2:574000/76‑1 (MQ=255) ccAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGtt < 1:586998/76‑1 (MQ=255) cAGTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGttt < 2:456436/76‑1 (MQ=255) gTTCTGGTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTg > 2:330421/1‑76 (MQ=255) ggTGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGaa > 2:620748/1‑76 (MQ=255) tGATATCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCCGAAAc > 2:178826/1‑76 (MQ=255) atatCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCa < 1:209317/76‑1 (MQ=255) tatCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAg < 2:15799/76‑1 (MQ=255) tCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGtt > 1:341978/1‑76 (MQ=255) tCTGCCACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGANGTTTTGCTGAAACCAGtt > 1:752932/1‑76 (MQ=255) cACTGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGa < 2:704466/76‑1 (MQ=255) tGAAAATATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCt > 1:333264/1‑76 (MQ=255) aaaTATGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGaa < 1:530175/76‑1 (MQ=255) aaTATGGCA‑GATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATc > 2:187667/1‑76 (MQ=255) atGGCAGGATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGc > 2:931339/1‑76 (MQ=255) | AAACTGCTGGAGCATATACGTAACAGAGGCGACCAGACCCCAGTTCTGGTGATATCTGCCACTGAAAATATGGCA‑GATATTGCCAAAGCGTTACGTCTGGGCGTTGAAGATGTTTTGCTGAAACCAGTTAAAGATCTGAATCGC > REL606/1290237‑1290380 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |