Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,004,106 | 0 | A | C | 28.1% | 147.0 / 12.4 | 64 | intergenic (‑32/+570) | ompF/asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base A (38/8); new base C (0/18); total (38/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.34e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAT > REL606/1004031‑1004178 | aaCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 2:96742/1‑76 (MQ=255) aaCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:114715/1‑76 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACAcc > 2:375316/1‑76 (MQ=255) agagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCt > 2:611194/1‑76 (MQ=255) gagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTg < 2:1005088/76‑1 (MQ=255) agCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCNTGGTTTTTTTTATGACACCTGc < 2:236338/76‑1 (MQ=255) gCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGcc < 2:596391/76‑1 (MQ=255) cAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCa > 2:19308/1‑76 (MQ=255) gACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGc < 2:1174518/76‑1 (MQ=255) gATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGt > 2:709489/1‑76 (MQ=255) tCACTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCa < 1:1208456/76‑1 (MQ=255) tCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCa > 2:198618/1‑76 (MQ=255) tCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCa > 2:106139/1‑76 (MQ=255) aCTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAAt > 1:179902/1‑76 (MQ=255) gCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAg > 2:85646/1‑76 (MQ=255) ccAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGt > 1:502544/1‑76 (MQ=255) aGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTc > 2:810360/1‑76 (MQ=255) gAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGCCACCTGCCACTGCCGTCAATAAGTTCt < 2:1736491/76‑1 (MQ=255) aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTg > 2:1043321/1‑76 (MQ=255) aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTg > 1:1242354/1‑76 (MQ=255) atatTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGt < 1:1526658/76‑1 (MQ=255) atatTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGt > 1:363297/1‑76 (MQ=255) atatTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGt > 1:1476050/1‑76 (MQ=255) atTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCa < 2:1126126/76‑1 (MQ=255) ttGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCaa < 1:1204023/76‑1 (MQ=255) ttGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCaa < 1:659175/76‑1 (MQ=255) tGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAAt < 1:815064/76‑1 (MQ=255) gcgcTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAGTa < 2:1297463/76‑1 (MQ=255) gcgcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATa > 1:638952/1‑76 (MQ=255) cTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaa < 2:653630/76‑1 (MQ=255) ttCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaaa > 1:186024/1‑76 (MQ=255) ttCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaaa > 2:191799/1‑76 (MQ=255) catcatTTTTTATTCCCCCCATGTTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAtt < 1:1639217/76‑1 (MQ=255) atcatTATTTATTCCCCTCATGTTTTTTTTTATAACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAttt < 2:1502394/76‑1 (MQ=255) catTATTTTTTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTAc < 2:319102/76‑1 (MQ=255) catTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTAc > 1:1691889/1‑76 (MQ=255) attattTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg > 2:712870/1‑76 (MQ=255) attattTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg > 2:712879/1‑76 (MQ=255) attattTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg > 2:712884/1‑76 (MQ=255) ttatttatTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg > 1:475515/1‑76 (MQ=255) tatttattCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa < 2:1262307/76‑1 (MQ=255) tttattACCCTCATGGTTTTTTNTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc > 2:238061/1‑76 (MQ=255) ttattcCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACt < 1:1378960/76‑1 (MQ=255) ttattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACt > 2:1430703/1‑76 (MQ=255) ttattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACt > 1:926303/1‑76 (MQ=255) tattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa > 2:1734846/1‑76 (MQ=255) tccccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg < 1:793937/74‑1 (MQ=255) taCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg > 2:3416/1‑76 (MQ=255) taCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg > 1:1304245/1‑76 (MQ=255) ccccTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga < 1:448973/75‑1 (MQ=255) cccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgat < 1:1582189/76‑1 (MQ=255) ccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatg < 1:1675004/76‑1 (MQ=255) cTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatga < 1:211750/76‑1 (MQ=255) cTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatga < 1:1712846/76‑1 (MQ=255) tCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag < 1:562233/76‑1 (MQ=255) aTGTTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag > 1:1452256/1‑76 (MQ=255) aTGATTTTTTTTATGACACCTGCCACTGCCGCCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag > 1:1452260/1‑76 (MQ=255) tttttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATGGATGAGAGttt < 2:1265793/75‑1 (MQ=255) ttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgt < 1:1033514/76‑1 (MQ=255) ttttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgt > 1:924454/1‑76 (MQ=255) tttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtg < 2:297285/76‑1 (MQ=255) tttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtg > 2:566302/1‑76 (MQ=255) tttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTc > 2:330779/1‑76 (MQ=255) tttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTc > 1:164802/1‑76 (MQ=255) ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCt > 1:1268195/1‑76 (MQ=255) ttATGACACCTGCCACTGCCGTCAATAAATTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCt > 1:1645865/1‑76 (MQ=255) tcggACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAGCTATTGATGAGAGTTTGGTGTCtt < 1:1114494/73‑1 (MQ=255) tATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCtt < 2:378252/76‑1 (MQ=255) tATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCtt > 2:167267/1‑76 (MQ=255) gACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAt < 2:789568/76‑1 (MQ=255) gACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAt > 2:1296674/1‑76 (MQ=255) | AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAT > REL606/1004031‑1004178 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |