Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 1,004,259 | (A)5→6 | intergenic (‑185/+417) | ompF ← / ← asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,004,254 | 1 | . | A | 85.2% | 86.2 / 8.3 | 27 | intergenic (‑180/+422) | ompF/asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base . (1/3); new base A (12/11); total (13/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.96e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTC‑A‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCC > REL606/1004179‑1004317 | gtgtCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTc < 1:760959/76‑1 (MQ=255) tGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAATATAAAGATTTGGTTTC‑a‑a‑a‑a‑a > 2:378113/1‑76 (MQ=255) tGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTC‑a‑a‑a‑a‑a > 2:12962/1‑76 (MQ=255) cAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTC‑a‑a‑a‑a‑aag > 2:216485/1‑74 (MQ=255) aGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTC‑a‑a‑a‑a‑aaga > 2:61010/1‑73 (MQ=255) ccATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTa < 1:550000/76‑1 (MQ=255) ccATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTa < 2:466561/76‑1 (MQ=255) cATTAAAACTAACGCTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTC‑AGA‑A‑A‑AGAAAAATATGTaa < 1:64545/76‑1 (MQ=255) aaCTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCaaaa < 1:59195/76‑1 (MQ=255) aCTAACGTTACGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAg < 2:53600/76‑1 (MQ=255) aaCGTTTCGCTACCATGAAAGTGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAGTATGTATCCAAAAGTaa < 1:70573/76‑1 (MQ=255) gTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAAt < 2:511975/76‑1 (MQ=255) tCGCTACCGTGNTAGTGCTACAAAGATAAAGATTAGATTGC‑A‑ATN‑A‑AGAAGAATATGTAACCGAAACTAAAATTTa < 1:97604/76‑1 (MQ=255) cTACCGTGAAAGTGCTACAAAGTTAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTCACCAAAAGTAAAATTTAAgg < 2:68568/76‑1 (MQ=255) cTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAgg > 1:294156/1‑76 (MQ=255) cTACCGTGAAAGCGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATATAAgg > 2:31936/1‑76 (MQ=255) tACCGTGAAAGTNCTNCAAAGATAAAGATTTGGTNTC‑A‑A‑A‑ACAGAAAAATATGTAACCAAAAGTAAAATTTAAGGa < 1:689479/76‑1 (MQ=255) tACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGa > 1:643887/1‑76 (MQ=255) aCCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGaa > 1:633731/1‑76 (MQ=255) tGAAAGTGGTACAAAGATAAAGATTTGGTTTC‑A‑A‑ATA‑AGAAAAATACGTATCCATAAGTAGAATTTAAGGAACttt > 2:15698/1‑76 (MQ=255) aaaGTGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtgt < 1:557163/76‑1 (MQ=255) gTGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGaa > 2:366954/1‑76 (MQ=255) tGCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAac < 1:411158/76‑1 (MQ=255) gCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAaca > 1:598434/1‑76 (MQ=255) gCTACAAAGATAAAGATTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAaca > 2:560306/1‑76 (MQ=255) aaagataaagatTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGACCTTTGTGAACACCgca > 2:186240/1‑74 (MQ=255) agataaagatTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCat < 2:598434/76‑1 (MQ=255) taaagatTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATAtt > 1:648773/1‑76 (MQ=255) aaagatTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATAttt > 1:159187/1‑76 (MQ=255) aagatTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTc > 2:707895/1‑76 (MQ=255) aagatTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTc > 2:707884/1‑76 (MQ=255) agatTTGGTTTCAA‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTcc < 1:666472/76‑1 (MQ=255) | GTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTC‑A‑A‑A‑A‑AGAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCC > REL606/1004179‑1004317 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |