Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,004,106 | 0 | A | C | 29.3% | 163.2 / 11.6 | 76 | intergenic (‑32/+570) | ompF/asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base A (41/12); new base C (0/22); total (41/35) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.06e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGT > REL606/1004032‑1004180 | aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac > 1:1331754/1‑76 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACAcc > 2:1114710/1‑76 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACAcc < 1:868789/76‑1 (MQ=255) agagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCt < 2:186323/76‑1 (MQ=255) agagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCt < 1:769219/76‑1 (MQ=255) gCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGcc < 1:1550488/76‑1 (MQ=255) gCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGcc < 1:985629/76‑1 (MQ=255) aGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCAc > 2:178807/1‑76 (MQ=255) aGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCAc > 2:1364415/1‑76 (MQ=255) aGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCAc > 2:798081/1‑76 (MQ=255) aCGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGcc > 1:1114319/1‑76 (MQ=255) aCGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGcc > 2:1267210/1‑76 (MQ=255) ggATCCCTCCCAGAATTTTGCGTTTCATCTTTTTTTTTTCCCCTCATGTTTTTTTTTATGACACCTGCCACTGCCg < 2:741122/75‑1 (MQ=255) cGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGTTTTTTTTTATGACACCTGCCACTGCCg < 2:10779/76‑1 (MQ=255) cGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCg > 2:70530/1‑76 (MQ=255) cGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCg > 2:292312/1‑76 (MQ=255) cGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCg > 2:320202/1‑76 (MQ=255) gATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGt < 1:760396/76‑1 (MQ=255) aTCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTc > 1:348071/1‑76 (MQ=255) aCTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAAt > 1:551549/1‑76 (MQ=255) tGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATaa < 1:1104672/76‑1 (MQ=255) ccAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGt > 2:70313/1‑76 (MQ=255) ccAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGt > 2:1439183/1‑76 (MQ=255) cAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGtt < 1:744038/76‑1 (MQ=255) cAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGtt < 1:1513709/76‑1 (MQ=255) cAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGtt < 1:408745/76‑1 (MQ=255) cAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGtt < 1:1139041/76‑1 (MQ=255) aGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTc > 2:192548/1‑76 (MQ=255) aGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTc > 2:397379/1‑76 (MQ=255) gAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCt < 1:1160119/76‑1 (MQ=255) aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTg > 2:1029033/1‑76 (MQ=255) aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTg > 1:1332139/1‑76 (MQ=255) atatTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGt < 2:56164/76‑1 (MQ=255) tatTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTc > 2:648010/1‑76 (MQ=255) tatTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTc > 2:648022/1‑76 (MQ=255) tatTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTc > 2:822233/1‑76 (MQ=255) gcgcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATa > 1:1031325/1‑76 (MQ=255) gcTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaaa < 1:401023/76‑1 (MQ=255) gcTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaaa < 1:567082/76‑1 (MQ=255) cTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaa > 2:1135010/1‑76 (MQ=255) ttCATCTTTTTTTATTCCCCTCAGGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaaa < 2:973945/76‑1 (MQ=255) ttCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaaa < 2:1610133/76‑1 (MQ=255) atcatTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAttt < 1:397686/76‑1 (MQ=255) atcatTATTTATTACCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAttt < 1:1614646/76‑1 (MQ=255) catTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTAc < 1:1349748/76‑1 (MQ=255) catTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTAc > 1:71541/1‑76 (MQ=255) tttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc > 1:66808/1‑76 (MQ=255) ttattcCCCTCATGTTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACt < 1:403162/76‑1 (MQ=255) ttattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACt > 2:1058638/1‑76 (MQ=255) tattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa > 1:16363/1‑76 (MQ=255) ttaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt > 2:208240/1‑76 (MQ=255) tccccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg < 1:342432/74‑1 (MQ=255) tccccTCATGTTTTTTTTTATCACGCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATGTACGGAACTATtg < 1:574177/74‑1 (MQ=255) taCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg > 1:361934/1‑76 (MQ=255) cccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgat < 1:626296/76‑1 (MQ=255) ccTCTTGGTTTTTTTTGTGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatg < 2:1097938/76‑1 (MQ=255) ccTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatg < 1:1554376/76‑1 (MQ=255) cTCATTGTTTTTTTTATCACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatga < 2:339430/76‑1 (MQ=255) cTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatga < 2:1190094/76‑1 (MQ=255) tgctGGTTTTTTTTATGGTACCTGCCCCTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag < 1:1496994/73‑1 (MQ=255) tCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag < 1:272051/76‑1 (MQ=255) tCATGTTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag < 2:1331754/76‑1 (MQ=255) tCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag < 1:1009925/76‑1 (MQ=255) tCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag > 2:220098/1‑76 (MQ=255) cATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga > 2:26145/1‑76 (MQ=255) cATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga > 2:1181716/1‑76 (MQ=255) cATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga > 2:1231966/1‑76 (MQ=255) aTGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag > 1:1039069/1‑76 (MQ=255) aTGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag > 1:608840/1‑76 (MQ=255) tGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGt < 1:290328/76‑1 (MQ=255) gtttttttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGtt < 1:1306024/74‑1 (MQ=255) gtttttttttAAAGCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGtt < 1:579049/74‑1 (MQ=255) ggTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGtt < 1:1295003/76‑1 (MQ=255) ggTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGtt < 1:1219917/76‑1 (MQ=255) ggTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGtt > 2:11524/1‑76 (MQ=255) ggTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGtt > 1:407367/1‑76 (MQ=255) ttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtgt < 1:1633910/76‑1 (MQ=255) ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCt > 1:1625689/1‑76 (MQ=255) ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCt > 1:1643176/1‑76 (MQ=255) tATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCtt > 2:576938/1‑76 (MQ=255) tATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCtt > 2:944733/1‑76 (MQ=255) aTGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg < 2:86593/57‑1 (MQ=255) aTGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg > 1:86593/1‑57 (MQ=255) cacCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAtgt > 2:1510478/1‑76 (MQ=255) | ACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGT > REL606/1004032‑1004180 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |