Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,004,106 | 0 | A | C | 44.2% | 63.6 / 26.4 | 52 | intergenic (‑32/+570) | ompF/asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base A (24/5); new base C (0/23); total (24/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.00e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGT > REL606/1004031‑1004180 | aaCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca < 2:1419727/76‑1 (MQ=255) aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac > 2:943576/1‑76 (MQ=255) aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac > 2:1065264/1‑76 (MQ=255) aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac > 2:653730/1‑76 (MQ=255) aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac < 2:515917/76‑1 (MQ=255) aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac > 2:1720915/1‑76 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACAcc < 1:717252/76‑1 (MQ=255) agCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc < 1:1839796/76‑1 (MQ=255) gCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGcc < 2:1403916/76‑1 (MQ=255) cAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCa > 2:1447449/1‑76 (MQ=255) cGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCg > 2:1096773/1‑76 (MQ=255) aTCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTc > 2:1669300/1‑76 (MQ=255) tcccTGCCAAAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCa < 1:1441968/73‑1 (MQ=255) cACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCaa < 1:1857755/76‑1 (MQ=255) cACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCaa < 2:1194863/76‑1 (MQ=255) ccTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAAt < 1:1859730/75‑1 (MQ=255) aCTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAAt > 2:1767371/1‑76 (MQ=255) cTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATa > 1:359344/1‑76 (MQ=255) ccAAAATTTTGCGCTTCATCTTTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGt < 2:753807/76‑1 (MQ=255) gAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCt < 2:819299/76‑1 (MQ=255) gAATATTGCGCTTCATCATGATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCt < 2:1292070/76‑1 (MQ=255) aaaTATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCt < 2:1153453/75‑1 (MQ=255) ttGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCaa < 1:1294205/76‑1 (MQ=255) cgcTTCATCTTTATTTATTACCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaa < 2:617792/76‑1 (MQ=255) cgcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaa > 2:957233/1‑76 (MQ=255) cgcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaa > 2:957228/1‑76 (MQ=255) gcTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaaa < 2:1713946/76‑1 (MQ=255) gcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaa > 1:1730242/1‑76 (MQ=255) cTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaa < 1:1916474/76‑1 (MQ=255) tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAt > 1:1303106/1‑76 (MQ=255) catcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAtt > 2:1603279/1‑76 (MQ=255) atcatTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAttt < 1:489481/76‑1 (MQ=255) atcatTATTTATTCCCCTCATGTTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAttt < 1:1821489/76‑1 (MQ=255) atcatTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAttt < 1:653934/76‑1 (MQ=255) catTATTTATTACCCTCATGG‑TTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg > 1:593710/1‑76 (MQ=255) ttatttatTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg > 1:124121/1‑76 (MQ=255) ttattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACt > 1:1783884/1‑76 (MQ=255) tattcCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa < 2:609378/76‑1 (MQ=255) tattaCCCTCATGG‑TTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt > 2:1718710/1‑76 (MQ=255) ttaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt > 1:780799/1‑76 (MQ=255) tccccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg < 1:669432/74‑1 (MQ=255) tccccTCAGGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg < 1:1532543/74‑1 (MQ=255) taCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg > 1:713950/1‑76 (MQ=255) tCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag < 1:740466/76‑1 (MQ=255) cATGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga < 1:1039556/76‑1 (MQ=255) cATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga > 1:555728/1‑76 (MQ=255) aTGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag > 2:1912454/1‑76 (MQ=255) tgtttttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGt < 1:407934/73‑1 (MQ=255) tgtttttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGt < 2:280368/73‑1 (MQ=255) gTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGttt < 1:1908073/76‑1 (MQ=255) ttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACNGAACTATTGATGAGAGTTTGgt < 2:262017/76‑1 (MQ=255) ttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgt < 1:217683/76‑1 (MQ=255) ttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgt < 1:380552/76‑1 (MQ=255) ttttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgt > 1:728261/1‑76 (MQ=255) tttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtg > 2:253172/1‑76 (MQ=255) ttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtgt < 2:395602/76‑1 (MQ=255) ttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtgt > 1:1929896/1‑76 (MQ=255) tttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTc < 1:1743168/76‑1 (MQ=255) tATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCtt > 1:1035251/1‑76 (MQ=255) aTGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCttt > 2:11443/1‑76 (MQ=255) gACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAt < 2:359344/76‑1 (MQ=255) gACCCCTGCCACTGCCGTCAANAAGTTCTGTCAATAAATATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAt < 2:314396/76‑1 (MQ=255) gACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAt > 2:464547/1‑76 (MQ=255) cacCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAtgt < 1:1578668/76‑1 (MQ=255) cacCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAtgt > 2:1151480/1‑76 (MQ=255) cacCTGCCACTGCCGTCAATAAGTTCTGCCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAtgt > 2:1151509/1‑76 (MQ=255) | AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGT > REL606/1004031‑1004180 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |