| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,004,106 | 0 | A | C | 32.1% | 115.1 / 12.7 | 54 | intergenic (‑32/+570) | ompF/asnS | outer membrane porin 1a (Ia;b;F)/asparaginyl‑tRNA synthetase |
| Reads supporting (aligned to +/- strand): ref base A (34/2); new base C (0/17); total (34/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.29e-12 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATG > REL606/1004031‑1004179 | aaCAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGaca < 1:114655/76‑1 (MQ=255) gagCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTg > 1:906961/1‑76 (MQ=255) cAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACTCCTGCCa < 2:161771/76‑1 (MQ=255) gATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGt > 1:390714/1‑76 (MQ=255) gATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGt > 1:390708/1‑76 (MQ=255) gATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGt > 1:390689/1‑76 (MQ=255) aTCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTc > 2:382031/1‑76 (MQ=255) tCACTGCCATAATATTGCGCTTCATCATTATTTCTTACCCTCATGGTTTTTTTTATGCCACCTGCCACTGCCGTCa > 2:1538309/1‑76 (MQ=255) tCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCa > 2:1538328/1‑76 (MQ=255) tCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCa > 2:546515/1‑76 (MQ=255) cccTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCCCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCaa < 2:1527963/74‑1 (MQ=255) ccTGCCAGAATATTGCGCTTCATCATTATTTATTCCCCTCATGTTTTTTTTTATGACACCTGCCACTGCCGTCAAt < 1:388103/75‑1 (MQ=255) ccTGCCAAAATATTGCGCTTCATCATTATTTATTCCCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAAt < 2:522114/75‑1 (MQ=255) cTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATa < 2:369141/76‑1 (MQ=255) cTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATa > 2:1094409/1‑76 (MQ=255) cTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATa > 1:722362/1‑76 (MQ=255) gCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAg > 2:504447/1‑76 (MQ=255) gCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAg > 2:805838/1‑76 (MQ=255) gCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAg > 2:407503/1‑76 (MQ=255) cAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGtt > 2:721500/1‑76 (MQ=255) cAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGtt > 2:721487/1‑76 (MQ=255) aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTg > 1:1085996/1‑76 (MQ=255) aaTATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTg > 1:1053548/1‑76 (MQ=255) ttGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCaa < 1:1379593/76‑1 (MQ=255) gcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaa > 2:325309/1‑76 (MQ=255) ttCGTCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaat > 1:1201979/1‑75 (MQ=255) ttCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATaaaaa > 1:1201954/1‑76 (MQ=255) tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAt > 2:761880/1‑76 (MQ=255) tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAt > 2:796443/1‑76 (MQ=255) tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAt > 2:109757/1‑76 (MQ=255) cataatTATTTATTCCCCCCATGGTTTTTTTTAGGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAAtt < 2:1264222/76‑1 (MQ=255) tcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTa > 1:260988/1‑76 (MQ=255) tatttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa > 1:1664601/1‑76 (MQ=255) tatttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGa > 1:1562320/1‑76 (MQ=255) tttattCCCCCAATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc < 1:198711/76‑1 (MQ=255) tttattACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAAc > 2:812245/1‑76 (MQ=255) ttattcCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACt < 1:242354/76‑1 (MQ=255) ttattcCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACt < 2:135255/76‑1 (MQ=255) ttattcCCCTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTACTGTCAATAAAAATTTACGGAACt < 1:1618095/76‑1 (MQ=255) ttccccTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt < 1:836009/73‑1 (MQ=255) ttaCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTAtt > 1:538895/1‑76 (MQ=255) tccccTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg < 1:1783496/74‑1 (MQ=255) tccccTCATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtg < 1:1177727/74‑1 (MQ=255) aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga > 1:968896/1‑76 (MQ=255) aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga > 1:968883/1‑76 (MQ=255) aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga > 1:473454/1‑76 (MQ=255) aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga > 1:209372/1‑76 (MQ=255) aCCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtga > 1:1548273/1‑76 (MQ=255) cTCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATtgatga < 1:1002704/76‑1 (MQ=255) tCATGTTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgag < 1:947900/76‑1 (MQ=255) cATGGTTTTTTTTATGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATgaga < 2:1481997/76‑1 (MQ=255) tGGTTTTTTTAATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGt < 1:1049833/76‑1 (MQ=255) ggTTTTTTTTAGGCCGCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGtt < 1:821761/76‑1 (MQ=255) tttttttATGCCCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTgg < 2:53983/76‑1 (MQ=255) tATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCtt > 1:915707/1‑76 (MQ=255) tGACCCCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTa < 2:1452454/76‑1 (MQ=255) acacCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTAtg > 2:411107/1‑76 (MQ=255) | AACAGAGCAGGGACGATCACTGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATG > REL606/1004031‑1004179 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |