Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 886,893 | C→A | G334G (GGC→GGA) | ECB_00825 → | putative replication protein for prophage |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 886,893 | 0 | C | A | 97.4% | 141.9 / ‑3.2 | 39 | G334G (GGC→GGA) | ECB_00825 | putative replication protein for prophage |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base A (18/20); total (18/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01 |
TGATGCATAAACGTGAGCAGCGCCGTAAGTCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGCGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCATCGCCGCAATGAGATGATG > REL606/886820‑886965 | tGATGCATAAACGTGAGCAGCGCCGTAAGTCGCTGGTTTTTTTTCGTNCTCATGAACTGGTGAACGAAGACGgaga > 1:124186/1‑76 (MQ=255) gCATAATCGTGAGCAGCGCCGTAAGTCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACg < 2:1388278/76‑1 (MQ=255) cATAAACGNGAGCAGCGCCNTAAGTCGCTGGAGTTTTTTCGTNCTCATGAACTGGTGAACGAAGACGGAGACACGc < 2:1186679/76‑1 (MQ=255) aTAAACGTTAGCAGCGCCGTAAGTCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCt < 1:1258958/76‑1 (MQ=255) aTAAACGTGCGCAGCGCCGTAAGTCGCTGGAGTGTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCt > 2:1219230/1‑76 (MQ=255) aTAAACGTGAGCAGCGCCGTAAGTCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGAAGACACGCt > 1:47863/1‑76 (MQ=255) aCGTGAGCAGCGCCGTAAGTCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGat < 2:333410/76‑1 (MQ=255) gTGAGCAGCGCCGTAAGTCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGatat < 2:968339/76‑1 (MQ=255) gTGAGCAGCGCCGTAAGTCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGatat < 2:301334/76‑1 (MQ=255) cagcGCCGTAAGTCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAg < 1:1110982/76‑1 (MQ=255) cGTAAGTCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTgg > 1:1111460/1‑76 (MQ=255) gTAAGTCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGCGACACGCTGGATATGGAAGATGTGGt < 2:884204/76‑1 (MQ=255) gTAAGTCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGt > 2:597891/1‑76 (MQ=255) aaGTCGCTGGAGTTTTTTCGGTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTaa < 2:554421/76‑1 (MQ=255) tCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACg < 2:1084580/76‑1 (MQ=255) cGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGc > 2:61942/1‑76 (MQ=255) cGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGc > 2:496826/1‑76 (MQ=255) gCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGcc < 2:671981/76‑1 (MQ=255) ggAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCcagc > 1:240852/1‑76 (MQ=255) gAGTTTCTTCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGGTGGATATGGAAGATGTGGTAAACGCcagca < 1:1302886/76‑1 (MQ=255) ttttttCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCcagcagca < 2:1391425/76‑1 (MQ=255) tttCGTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAgcagcagca < 2:462586/76‑1 (MQ=255) gTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAAccc > 1:790879/1‑76 (MQ=255) gTTCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAAccc > 2:435577/1‑76 (MQ=255) ttCTCATGAACTGGTGAACGAAGACGGAGACACGCTNGATATGGAAGCTGTGGTAAACGCCAGCAGCAGCANCCCg > 1:173816/1‑76 (MQ=255) ttCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCg > 1:1268245/1‑76 (MQ=255) ttCTCATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCg > 2:1104061/1‑76 (MQ=255) ctcATGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCgcc > 2:26757/1‑74 (MQ=255) aTGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCa > 2:876273/1‑76 (MQ=255) tGAACTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCAt < 2:892435/76‑1 (MQ=255) gAACTGGTGAACGAAGACGGAGACACGTTGGATATGGAAGATGANGTNAACGCCAGCAGCAGNAACCCGGCGCATc > 2:1282622/1‑76 (MQ=255) aaCTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCATcg < 2:563464/76‑1 (MQ=255) aaCTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCATcg < 2:942902/76‑1 (MQ=255) cTGGTGAACGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCATcgcc > 2:793870/1‑76 (MQ=255) cGAAGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCATCGCCGCAATgag < 1:933337/76‑1 (MQ=255) gaGACGGAGACACNCTGGATATGGTAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCATCGCCGCAATGAgat < 1:190990/75‑1 (MQ=255) aaGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCATCGCCGCAATGAgat < 2:1222575/76‑1 (MQ=255) aGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCATCGCCGCAATGAgatg < 2:1026506/76‑1 (MQ=255) aGACGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGAAACCCGGCGCACCGCCGCACTgggagg > 2:573829/1‑74 (MQ=255) cGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCATCGCCGCAATGAGatgatg < 2:436099/76‑1 (MQ=255) cGGAGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCATCGCCGCAATGAGatgatg > 2:964348/1‑76 (MQ=255) | TGATGCATAAACGTGAGCAGCGCCGTAAGTCGCTGGAGTTTTTTCGTTCTCATGAACTGGTGAACGAAGACGGCGACACGCTGGATATGGAAGATGTGGTAAACGCCAGCAGCAGCAACCCGGCGCATCGCCGCAATGAGATGATG > REL606/886820‑886965 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |