Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,570,909 | T→A | I37N (ATC→AAC) | uspA → | universal stress global response regulator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,570,909 | 0 | T | A | 97.5% | 152.5 / ‑2.5 | 40 | I37N (ATC→AAC) | uspA | universal stress global response regulator |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base A (8/31); total (8/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-01 |
CCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACG > REL606/3570837‑3570981 | ccccGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCa < 2:1545532/76‑1 (MQ=255) cccGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCAc < 2:1146831/76‑1 (MQ=255) ccGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACCATGCGAAAGTTTCTCTGAACCACg > 1:79085/1‑76 (MQ=255) aaaGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGGCCGTACAATGAGGAAGCTTCTCTGTGCCACGTAgc > 2:220663/1‑75 (MQ=255) aaGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGAt < 2:1240206/76‑1 (MQ=255) aaGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGAt < 2:396878/76‑1 (MQ=255) aaGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGAt < 1:621060/76‑1 (MQ=255) cAAAGTTCTGGGGGAAACAGCCGTTACTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTa < 2:669648/76‑1 (MQ=255) aaaGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTaa > 1:1585246/1‑76 (MQ=255) aaGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTaaa < 1:1270551/76‑1 (MQ=255) aaGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTaaa > 2:1504109/1‑76 (MQ=255) aaGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTaaa < 2:693171/76‑1 (MQ=255) aGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAac < 2:186141/76‑1 (MQ=255) aGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCCACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAac > 2:1006/1‑76 (MQ=255) ttCTGGTAGAGAAAGCAGTCTCTATGGCTCTCCCCTTCAACAAGTACGCTTTTCTGAACCACGAAGGTGTTAacta > 2:475985/1‑76 (MQ=255) cTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAAAGCGAAAATTTCTCTTAACCACGTAGATGTAAaataat > 1:784873/1‑74 (MQ=255) gTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTg < 1:1392527/76‑1 (MQ=255) tAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGa < 2:951211/76‑1 (MQ=255) tAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGa > 2:801784/1‑76 (MQ=255) tAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGa < 1:740929/76‑1 (MQ=255) tAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGa < 1:217903/76‑1 (MQ=255) aaaGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCtat < 2:683206/76‑1 (MQ=255) aaaGCAGTCTCTATGGCTCGCCCCAACAATGCGAAAGTTTCTCTCACCAAAGGAGATGGAAACTACTCTGCCCtat > 2:191227/1‑76 (MQ=255) aaGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCtata > 1:658856/1‑76 (MQ=255) gCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATacn < 1:1386087/76‑2 (MQ=255) gCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATaca < 1:208312/76‑1 (MQ=255) gTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCg < 2:141817/76‑1 (MQ=255) tctcTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCgg < 1:109557/76‑1 (MQ=255) ctATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCt < 1:566567/76‑1 (MQ=255) tATGGCTCCCCCCCACAATGCGAAAGTTTTTCTGAACCACGTAGATGTAAAATACTCTGACCTATACACCGGGcct > 2:684937/1‑74 (MQ=255) ggCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTAtt < 2:1426893/76‑1 (MQ=255) gCTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTg < 1:1568505/76‑1 (MQ=255) cTCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGa > 1:609896/1‑76 (MQ=255) tCGCCCCTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGAt < 1:1075774/76‑1 (MQ=255) gcccTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGAtgtg < 2:1023618/75‑1 (MQ=255) cccTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGa > 1:149493/1‑76 (MQ=255) cTACAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAAt < 1:1263398/76‑1 (MQ=255) aCAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCt < 2:590733/76‑1 (MQ=255) cAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTg > 1:334421/1‑76 (MQ=255) cAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTg > 1:725949/1‑76 (MQ=255) cAATGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTg < 2:1404214/76‑1 (MQ=255) aaTGCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTgg < 1:999046/76‑1 (MQ=255) gCGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTg < 2:663533/76‑1 (MQ=255) gCGAAAGTTTCTCGGAAATACGTAGATGTAAACTACTCTGACCTATTCACCGGGCTTATTGATGTGAATCAGGTTg < 2:830983/76‑1 (MQ=255) cGAAAGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGa > 1:556038/1‑76 (MQ=255) aaGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGatat > 1:1340834/1‑76 (MQ=255) aaGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGatat > 2:794959/1‑76 (MQ=255) aGTTTCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATg < 1:1148662/76‑1 (MQ=255) ttCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAg < 2:875753/76‑1 (MQ=255) ttCTCTGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAg < 2:875776/76‑1 (MQ=255) ctGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAAc < 2:1450409/76‑1 (MQ=255) tGAACCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACg > 2:1133194/1‑76 (MQ=255) | CCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACG > REL606/3570837‑3570981 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |