Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,004,106 | 0 | A | C | 31.2% | 136.0 / 14.1 | 64 | intergenic (‑32/+570) | ompF/asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base A (36/8); new base C (0/20); total (36/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.58e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAT > REL606/1004031‑1004178 | aaCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGaca < 2:923945/76‑1 (MQ=255) aaCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGaca < 2:1611481/76‑1 (MQ=255) aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGacac < 1:913090/76‑1 (MQ=255) aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGacac < 1:1801776/76‑1 (MQ=255) aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGacac < 1:266577/76‑1 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACAcc > 1:420581/1‑76 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACAcc < 2:863480/76‑1 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACAcc < 2:419697/76‑1 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACAcc > 2:364631/1‑76 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACAcc < 1:294327/76‑1 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACAcc > 1:1702436/1‑76 (MQ=255) gagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTg < 2:706487/76‑1 (MQ=255) gagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTg > 2:643063/1‑76 (MQ=255) gagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTg > 2:1686605/1‑76 (MQ=255) gagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTg < 2:1290249/76‑1 (MQ=255) gCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGcc < 1:10959/76‑1 (MQ=255) aGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTTATGACACCTGcaa > 2:1659507/1‑74 (MQ=255) ggACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGT‑TTTTTTTTATGACCCCTGCCACTg < 2:162991/76‑1 (MQ=255) ggACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTg > 2:929404/1‑76 (MQ=255) cGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCg < 2:92643/76‑1 (MQ=255) aTCACTGCCAGAATATTGCGCTTCTTCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTc > 2:708594/1‑76 (MQ=255) tCACTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGG‑TTTTTTTTATGACCCCTGCCACTGCCGTCa < 1:1390366/76‑1 (MQ=255) ctcccAAAATATTGCGCTTCATCATTATTTTTTACCCTCATGT‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATa < 2:592213/73‑1 (MQ=255) cTGCCAAAATATTGCGCTTCATCATTATTTATTCCCCTCATGG‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATa < 2:1639890/76‑1 (MQ=255) gCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAg > 2:441003/1‑76 (MQ=255) gCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAg > 2:691906/1‑76 (MQ=255) cAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGT‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGtc < 1:686956/76‑2 (MQ=255) aGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTCTGACACCTGCCACTGCCGTCAATAAGTTc > 2:312563/1‑76 (MQ=255) aGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTc > 1:38373/1‑76 (MQ=255) aGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTc > 2:1686577/1‑76 (MQ=255) aGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTc > 1:1761994/1‑76 (MQ=255) aGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTc > 2:1104067/1‑76 (MQ=255) gAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCt > 2:451469/1‑76 (MQ=255) gAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCt > 2:451475/1‑76 (MQ=255) aaTATTGCGCTTCATCATTATTTATTCCCCTCATGG‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTg < 1:179867/76‑1 (MQ=255) aaTATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTg > 2:1157010/1‑76 (MQ=255) tatTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTc > 1:108937/1‑76 (MQ=255) gcgcTTCATCATTATTTATTCCCCTCATGT‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATa < 1:1427709/76‑1 (MQ=255) gcgcTTCATCATTATTTATTCCCCTCATGG‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATa < 1:499492/76‑1 (MQ=255) gcgcTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATa > 1:1621143/1‑76 (MQ=255) cgcTTCATCATTATTTATTCCCCTCATGG‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaa < 1:527303/76‑1 (MQ=255) cgcTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaa > 2:515763/1‑76 (MQ=255) gcTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaa > 1:1257188/1‑76 (MQ=255) catTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTAc > 1:790511/1‑76 (MQ=255) catTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTAc > 1:790492/1‑76 (MQ=255) attattTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg > 2:894086/1‑76 (MQ=255) attattTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg > 2:490076/1‑76 (MQ=255) attattTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg > 2:374126/1‑76 (MQ=255) attattTATTACCCTCATGG‑TTNTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACg > 2:245032/1‑76 (MQ=255) ttttttCCCCTCATGT‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc < 2:668201/76‑1 (MQ=255) tttattCCCCTCATGT‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc < 2:1162573/76‑1 (MQ=255) tttattCCCCTCATGT‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc < 1:383093/76‑1 (MQ=255) tattaCCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa > 1:1240844/1‑76 (MQ=255) attcCCCTCATGT‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt < 1:1280722/76‑1 (MQ=255) ttccccTCATGT‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt < 1:1328266/73‑1 (MQ=255) ttaCCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt > 1:399882/1‑76 (MQ=255) tccccTCATGT‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg < 1:28840/74‑1 (MQ=255) aCCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga > 1:1360041/1‑76 (MQ=255) cccTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgat > 1:592209/1‑76 (MQ=255) cccTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgat > 1:592239/1‑76 (MQ=255) cTCATGG‑TCTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatga > 2:1707251/1‑76 (MQ=255) tCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag > 1:937302/1‑76 (MQ=255) cATGT‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga < 1:1202037/76‑1 (MQ=255) cATGT‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga < 1:814834/76‑1 (MQ=255) g‑ttgttTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGttt < 2:1134232/76‑1 (MQ=255) g‑TTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGttt < 2:1739781/76‑1 (MQ=255) tttttttttACGACCCCTGCCAATGCCGGCAATAAGTTCTGTTAATAAAAATTTACGGAACTATTGATGAGAGttt < 1:737976/75‑1 (MQ=255) ttttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTg < 1:662070/76‑1 (MQ=255) ttttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgt < 1:961421/76‑1 (MQ=255) ttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtgt < 2:1497916/76‑1 (MQ=255) ttttATCACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtgt < 2:333666/76‑1 (MQ=255) ttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCt < 2:177695/76‑1 (MQ=255) ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCt < 2:230305/76‑1 (MQ=255) tATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCtt < 2:437479/76‑1 (MQ=255) aTGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCttt > 1:765976/1‑76 (MQ=255) tGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTa > 2:461716/1‑76 (MQ=255) gACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAt < 2:197044/76‑1 (MQ=255) | AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGG‑TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAT > REL606/1004031‑1004178 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |