Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,004,106 | 0 | A | C | 42.1% | 71.0 / 26.4 | 57 | intergenic (‑32/+570) | ompF/asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base A (20/13); new base C (0/24); total (20/37) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.33e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTT > REL606/1004034‑1004175 | agagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACAcc < 1:1286839/75‑1 (MQ=255) agagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACAcc < 1:1528737/75‑1 (MQ=255) gagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCt < 1:158895/75‑1 (MQ=255) gagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCt < 1:1296676/75‑1 (MQ=255) gCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc < 1:1347317/75‑1 (MQ=255) aGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCa > 2:1707238/1‑75 (MQ=255) aGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCa > 2:48763/1‑75 (MQ=255) gggACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCTTGTTTTTTTTTATGCCCCCTGCCAc < 2:79325/75‑1 (MQ=255) aCGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGc > 1:1349771/1‑75 (MQ=255) gATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCg < 2:97370/75‑1 (MQ=255) ccTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCaa < 2:505978/74‑1 (MQ=255) ccTGCAAGAATTTTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTCaa < 1:585758/74‑1 (MQ=255) cTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACACCTGCCACTGCCGTCAAt < 2:1782254/75‑1 (MQ=255) cTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAAt < 1:595156/75‑1 (MQ=255) cTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGGCAAt > 2:2006221/1‑75 (MQ=255) tGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACACCTGCCACTGCCGTCAATa < 2:832022/75‑1 (MQ=255) ccAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAg > 1:1205539/1‑75 (MQ=255) ccAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAg > 1:1463331/1‑75 (MQ=255) gAATATTGCTCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTc < 1:967134/75‑1 (MQ=255) gAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTc < 1:562739/75‑1 (MQ=255) aaaTTTTGCGTTTCATCATTATTTATTCCCCTCATGTTTTTTTTTAGGCCCCCTGCCACTGCCGTCAATAAGTTc < 1:757658/74‑1 (MQ=255) aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCt > 2:997658/1‑75 (MQ=255) tatTGCGCTTCATCATTATTTATTACCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGt < 1:801685/75‑1 (MQ=255) atTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTc > 2:209484/1‑75 (MQ=255) ttGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCa > 2:1104017/1‑75 (MQ=255) tGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCaa < 1:719894/75‑1 (MQ=255) tGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCaa < 1:1336899/75‑1 (MQ=255) gcgcTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAAt < 1:342790/75‑1 (MQ=255) gcgcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAAt > 1:524462/1‑75 (MQ=255) ttCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaa < 1:849222/75‑1 (MQ=255) ttCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaa > 1:796021/1‑75 (MQ=255) tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaaa > 1:2109964/1‑75 (MQ=255) tatcatTTTTTTTTCCCCCCATGGTTTTTTTTTTGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAt < 2:252057/74‑1 (MQ=255) atcatTATTTATTCCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAtt < 2:289281/75‑1 (MQ=255) atcatTATTTATTCCCCTCAGGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAtt < 1:1711900/75‑1 (MQ=255) tctttatttatTCCCCTCATGTTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAttt < 2:484672/72‑1 (MQ=255) tcatTATTTATTCCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAttt < 2:484694/75‑1 (MQ=255) tcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAttt > 2:1244964/1‑75 (MQ=255) ttatttatTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg < 2:1021881/75‑1 (MQ=255) tatttattCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg < 2:1963328/75‑1 (MQ=255) tatttattCCCCTCATGGTTTTTTTTAGGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg < 1:622965/75‑1 (MQ=255) tatttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg > 2:2083300/1‑75 (MQ=255) tatttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg > 2:2083306/1‑75 (MQ=255) atttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa > 1:1434725/1‑75 (MQ=255) tttattCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa < 1:1707238/75‑1 (MQ=255) tttattACCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa < 2:215503/75‑1 (MQ=255) ttattaCCCTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc < 2:1784006/75‑1 (MQ=255) attaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa > 2:994156/1‑75 (MQ=255) ttgcccTCATGGTTTTTTTTATCACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt < 2:1757830/72‑1 (MQ=255) ttaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt > 1:1389625/1‑75 (MQ=255) tccccTCCTGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt < 2:880043/73‑1 (MQ=255) tccccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt < 2:1168387/73‑1 (MQ=255) aCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg < 2:1768079/75‑1 (MQ=255) cTCATGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatg < 2:342959/75‑1 (MQ=255) cTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatg < 1:1287092/75‑1 (MQ=255) cTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatg < 2:2098730/75‑1 (MQ=255) tCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatga > 2:1523285/1‑75 (MQ=255) ggTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGt < 2:484047/75‑1 (MQ=255) gTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGtt < 2:97543/75‑1 (MQ=255) tttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgt < 2:1897180/75‑1 (MQ=255) ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTc > 2:1786006/1‑75 (MQ=255) tATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCt < 1:92376/75‑1 (MQ=255) aTGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCtt > 1:2081868/1‑75 (MQ=255) | AGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTT > REL606/1004034‑1004175 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |