Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,004,106 | 0 | A | C | 39.4% | 114.3 / 24.4 | 66 | intergenic (‑32/+570) | ompF/asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base A (37/3); new base C (0/26); total (37/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.21e-15 | |||||||||||
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 | aaCAGAGCAGGGACGATCACTGCCAGATTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca > 1:9882/1‑76 (MQ=255) aCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGacac > 2:1112204/1‑76 (MQ=255) cAGAGCAGGGACGATCACTGCCAGACTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACAcc > 1:304226/1‑76 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACAcc > 2:1487905/1‑76 (MQ=255) cAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACAcc > 2:501801/1‑76 (MQ=255) agagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCt > 2:479708/1‑76 (MQ=255) gagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTg < 2:994349/76‑1 (MQ=255) agCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGc > 2:2061699/1‑76 (MQ=255) cAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCa > 2:1960717/1‑76 (MQ=255) aGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCAc > 1:151995/1‑76 (MQ=255) gggACGATCACTGCCAAAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACt < 2:1703200/76‑1 (MQ=255) aCGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGcc > 2:1365678/1‑76 (MQ=255) aTCACTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTc < 1:1920283/76‑1 (MQ=255) aTCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTc < 1:1781945/76‑1 (MQ=255) tcccTGCCAAAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCa < 1:1477666/73‑1 (MQ=255) tCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGCCCCCTGCCACTGCCGTCa < 2:242612/76‑1 (MQ=255) tGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATaa > 1:1073583/1‑76 (MQ=255) aGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTc > 1:447953/1‑76 (MQ=255) gAATTTTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCt < 1:817884/76‑1 (MQ=255) gAATTTTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCt < 2:679658/76‑1 (MQ=255) gAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCt < 1:1442110/76‑1 (MQ=255) aaaTATTGCGCTTCATCATTATTTTTTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCt < 1:838469/75‑1 (MQ=255) atatTGCGCTTCATCATTATTTATTCCCCTCAGGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGt < 1:1705718/76‑1 (MQ=255) tatTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTc > 2:779164/1‑76 (MQ=255) tGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAAt > 1:1127621/1‑76 (MQ=255) gcgcTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATa < 1:370330/76‑1 (MQ=255) gcgcTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATa < 2:1199635/76‑1 (MQ=255) cgcTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaa < 2:297603/76‑1 (MQ=255) gcTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaaa < 2:60648/76‑1 (MQ=255) gcTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATaaa < 2:1052868/76‑1 (MQ=255) ttCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCATTaaaaa > 2:377439/1‑76 (MQ=255) tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAt > 2:1348540/1‑76 (MQ=255) catcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAtt > 1:1643647/1‑76 (MQ=255) tcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTa > 1:839287/1‑76 (MQ=255) catTATTTATTCCCCTCATGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTAc < 1:1433333/76‑1 (MQ=255) ttatttatTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACgg > 1:1851410/1‑76 (MQ=255) tttttattCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa < 2:1416591/74‑1 (MQ=255) tatttattCCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa < 2:1536121/76‑1 (MQ=255) tatttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa > 1:1739362/1‑76 (MQ=255) atttattCCCCTCATGGTTTTTTTTATGCCACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa < 1:792697/76‑1 (MQ=255) atttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa > 2:1831763/1‑76 (MQ=255) atttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGaa > 2:933436/1‑76 (MQ=255) tattcCCCTCATGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa < 1:1511987/76‑1 (MQ=255) tattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa > 2:184848/1‑76 (MQ=255) tattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa > 2:1188724/1‑76 (MQ=255) tattaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTa > 2:1584284/1‑76 (MQ=255) attaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAt < 1:345151/76‑1 (MQ=255) ttccccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt < 2:895242/73‑1 (MQ=255) ttaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt > 1:1550704/1‑76 (MQ=255) tccccTCATGTTTTTTTTTATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg < 1:1594708/74‑1 (MQ=255) tccccTCATGGTTTTTTTTAAGCCGCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg < 1:448229/74‑1 (MQ=255) taCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg > 2:1046515/1‑76 (MQ=255) aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga > 1:249128/1‑76 (MQ=255) aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga > 2:401115/1‑76 (MQ=255) cccTCATGGTTTTTTTTATAACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgat < 1:926991/76‑1 (MQ=255) cTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatga < 2:1834135/76‑1 (MQ=255) cATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga < 1:254927/76‑1 (MQ=255) cATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga > 1:2028818/1‑76 (MQ=255) aTGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTGTGTCAATAAAAATTTACGGAACTATTGATGagag > 1:1932797/1‑76 (MQ=255) aTGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag > 1:750330/1‑76 (MQ=255) aTGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGagag > 1:1932819/1‑76 (MQ=255) ggTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGtt > 1:2012575/1‑76 (MQ=255) tttttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGttt < 2:186208/75‑1 (MQ=255) ttttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTg < 1:9419/76‑1 (MQ=255) ttttttttATCACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTg < 1:894769/76‑1 (MQ=255) ttttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgt < 1:952657/76‑1 (MQ=255) ttttATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtgt < 2:1600544/76‑1 (MQ=255) ttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGgtgt > 2:1242550/1‑76 (MQ=255) tttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTc > 1:1381656/1‑76 (MQ=255) ttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCt > 2:663912/1‑76 (MQ=255) tATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCtt < 2:1982198/76‑1 (MQ=255) tATCACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCtt < 2:151995/76‑1 (MQ=255) aTGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCttt > 2:1688737/1‑76 (MQ=255) gACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAt > 1:1317232/1‑76 (MQ=255) acacCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAtg > 2:1652154/1‑76 (MQ=255) cacCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAtgt > 2:353494/1‑76 (MQ=255) | AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGT > REL606/1004031‑1004180 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |