Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,004,106 | 0 | A | C | 35.4% | 128.0 / 20.1 | 65 | intergenic (‑32/+570) | ompF/asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base A (36/6); new base C (0/23); total (36/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.09e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTG > REL606/1004031‑1004181 | aaCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:20886/1‑76 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACAcc < 2:706207/76‑1 (MQ=255) agagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCt < 1:1475923/76‑1 (MQ=255) agCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc > 2:632975/1‑76 (MQ=255) agCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc > 2:58618/1‑76 (MQ=255) gCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGcc > 1:411674/1‑76 (MQ=255) gCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGcc > 2:147467/1‑76 (MQ=255) aGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCAc > 2:315795/1‑76 (MQ=255) ggACGATCTCTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTg < 2:519595/76‑1 (MQ=255) cGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCg > 2:193513/1‑76 (MQ=255) gATCACTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACACCTGCCACTGCCGt < 1:335664/76‑1 (MQ=255) tCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCa > 2:1094772/1‑76 (MQ=255) cTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATa > 1:1232843/1‑76 (MQ=255) tGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATaa < 1:1557362/76‑1 (MQ=255) tGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATaa < 1:1522984/76‑1 (MQ=255) gCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTTTGACACCTGCCCCTTCCGGCAATAAg > 2:1075011/1‑76 (MQ=255) gCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAg > 2:132455/1‑76 (MQ=255) aGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTc > 1:1100955/1‑76 (MQ=255) aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTg > 1:665231/1‑76 (MQ=255) atTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCa < 1:699372/76‑1 (MQ=255) atTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCa > 1:1047559/1‑76 (MQ=255) tGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAAt > 2:70222/1‑76 (MQ=255) gcgcTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATa < 1:1162045/76‑1 (MQ=255) cgcTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaa < 1:1323788/76‑1 (MQ=255) cgcTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaa < 1:925275/76‑1 (MQ=255) cTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaa > 2:159131/1‑76 (MQ=255) ttCATCATTATTTATTCCCCTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaaa < 2:818378/76‑1 (MQ=255) ttCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaaa < 2:949038/76‑1 (MQ=255) ttCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaaa < 1:615863/76‑1 (MQ=255) tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAt > 1:1202622/1‑76 (MQ=255) tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAt > 2:1111549/1‑76 (MQ=255) tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAt > 1:169298/1‑76 (MQ=255) tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAt > 2:1111567/1‑76 (MQ=255) tcatTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCATCAATAAGTTCTGTCAATAAAAATTTa < 1:1631311/76‑1 (MQ=255) attattTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg < 1:1249842/76‑1 (MQ=255) attattTATTACCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg < 1:41439/76‑1 (MQ=255) ttatttatTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg > 1:1338528/1‑76 (MQ=255) tatttattCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa < 1:901921/76‑1 (MQ=255) tatttattCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa < 1:943124/76‑1 (MQ=255) tatttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa > 1:521308/1‑76 (MQ=255) tatttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa > 2:375733/1‑76 (MQ=255) atttattCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa < 1:661673/76‑1 (MQ=255) atttattCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa < 1:265331/76‑1 (MQ=255) atttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa > 1:1585089/1‑76 (MQ=255) atttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa > 1:1585093/1‑76 (MQ=255) atttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa > 1:949955/1‑76 (MQ=255) tttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc > 1:231340/1‑76 (MQ=255) ttattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACt > 1:442233/1‑76 (MQ=255) attaCCCTCAGGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt < 2:36774/76‑1 (MQ=255) ttaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt > 1:358907/1‑76 (MQ=255) aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga > 2:528657/1‑76 (MQ=255) ccTCATGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatg < 1:1295277/76‑1 (MQ=255) tCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag < 1:284447/76‑1 (MQ=255) tCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag < 1:1055042/76‑1 (MQ=255) tCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag < 2:1514717/76‑1 (MQ=255) tCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag < 2:67468/76‑1 (MQ=255) tCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag > 2:113515/1‑76 (MQ=255) aTGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag < 1:428968/76‑1 (MQ=255) aTGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag < 1:1434080/76‑1 (MQ=255) gtttttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGtt < 2:425558/74‑1 (MQ=255) gTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGttt < 2:847182/76‑1 (MQ=255) ttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgt < 1:1518454/76‑1 (MQ=255) tttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTc > 1:340283/1‑76 (MQ=255) tttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTc > 1:731368/1‑76 (MQ=255) tttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGGGTc > 1:1283667/1‑76 (MQ=255) ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCt > 1:1579362/1‑76 (MQ=255) ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCt > 1:159907/1‑76 (MQ=255) ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCt > 1:793985/1‑76 (MQ=255) tGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTa > 1:1035660/1‑76 (MQ=255) cacCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAtgt < 2:881811/76‑1 (MQ=255) acCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAtgtg > 1:1054984/1‑76 (MQ=255) | AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTG > REL606/1004031‑1004181 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |