Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,004,106 | 0 | A | C | 28.8% | 134.3 / 11.4 | 59 | intergenic (‑32/+570) | ompF/asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base A (34/8); new base C (0/17); total (34/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.90e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTA > REL606/1004031‑1004177 | aaCAGAGCAGGGGCGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca < 1:1366994/76‑1 (MQ=255) aaCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 2:1332585/1‑76 (MQ=255) aaCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca < 1:312193/76‑1 (MQ=255) aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac < 2:949888/76‑1 (MQ=255) aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac > 1:644526/1‑76 (MQ=255) aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac < 2:738417/76‑1 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACAcc > 2:1107545/1‑76 (MQ=255) gagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTg < 2:1440018/76‑1 (MQ=255) gagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTg < 1:1537278/76‑1 (MQ=255) agCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc > 1:670654/1‑76 (MQ=255) agCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc < 2:197289/76‑1 (MQ=255) gCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGcc < 1:1227373/76‑1 (MQ=255) gCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGcc < 1:528932/76‑1 (MQ=255) gCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGcc < 2:791986/76‑1 (MQ=255) aGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCAc > 1:810273/1‑76 (MQ=255) aTCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTc > 2:330700/1‑76 (MQ=255) tGCCAGAATATTGCGCTTCATCTTTATTTATTACCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATaa < 2:1286014/76‑1 (MQ=255) tGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATaa < 2:1368955/76‑1 (MQ=255) tGCCAAAATATTCCGCTTCATCATTATTTATTCCCCTCAGGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATaa < 2:887134/76‑1 (MQ=255) gCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAg < 2:719844/76‑1 (MQ=255) gCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAg < 2:1491833/76‑1 (MQ=255) aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTg < 1:811188/76‑1 (MQ=255) atatTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGt > 1:1185337/1‑76 (MQ=255) atTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCa > 1:1335295/1‑76 (MQ=255) ttGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCaa > 1:809188/1‑76 (MQ=255) ttGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCaa > 1:1451172/1‑76 (MQ=255) gcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaa > 1:256295/1‑76 (MQ=255) tcatcaTTATTTATTCCCCTCAGGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAt < 1:276867/76‑1 (MQ=255) tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAt > 2:167630/1‑76 (MQ=255) atcatTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAttt < 1:822980/76‑1 (MQ=255) tcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTa > 2:896987/1‑76 (MQ=255) nattattTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTAc > 1:1413773/2‑76 (MQ=255) catTTTTTATTACCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTAc < 1:1534757/76‑1 (MQ=255) catTATTTTTTCCCCTCATGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTAc < 2:811104/76‑1 (MQ=255) ttatttatTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg > 2:787474/1‑76 (MQ=255) ttatttatTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg > 2:1046084/1‑76 (MQ=255) tatttattCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa < 1:607962/76‑1 (MQ=255) tatttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa > 2:881161/1‑76 (MQ=255) atttattACCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa < 2:1002246/76‑1 (MQ=255) tttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc > 2:1215400/1‑76 (MQ=255) tttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc > 1:100520/1‑76 (MQ=255) tattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa > 1:242473/1‑76 (MQ=255) tattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa > 1:1448650/1‑76 (MQ=255) tattaCCCTCATGGNTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa > 1:868883/1‑76 (MQ=255) tattaCCCTCATGGNTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa > 1:868896/1‑76 (MQ=255) attaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt > 1:1098589/1‑76 (MQ=255) attaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt > 1:307625/1‑76 (MQ=255) ttaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTCtt > 2:289515/1‑76 (MQ=255) ttaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt > 1:359083/1‑76 (MQ=255) taCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg > 2:48619/1‑76 (MQ=255) cccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgat < 1:744174/76‑1 (MQ=255) ccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatg < 1:561326/76‑1 (MQ=255) cTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatga < 2:91986/76‑1 (MQ=255) cTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatga > 2:414956/1‑76 (MQ=255) tCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag < 2:242473/76‑1 (MQ=255) tCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatgan < 1:1390972/76‑2 (MQ=255) cATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga < 2:644526/76‑1 (MQ=255) cATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga > 2:566233/1‑76 (MQ=255) aTGTTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag < 1:579816/76‑1 (MQ=255) aTGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag > 1:809433/1‑76 (MQ=255) tttttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTgg > 1:547879/1‑76 (MQ=255) tATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCtt > 1:455657/1‑76 (MQ=255) tATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCtt > 2:1475829/1‑76 (MQ=255) aTGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCttt > 2:1440795/1‑76 (MQ=255) tGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTa < 1:330700/76‑1 (MQ=255) tGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTa > 2:1173391/1‑76 (MQ=255) | AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTA > REL606/1004031‑1004177 |
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: ‑ |