Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 303,631 | 0 | A | G | 22.9% | 63.1 / 15.0 | 35 | intergenic (+26/+222) | yagU/ykgJ | DUF1440 family inner membrane acid resistance protein/UPF0153 cysteine cluster protein |
Reads supporting (aligned to +/- strand): ref base A (13/14); new base G (4/4); total (17/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGGACATTTAGTCTGGTTCTGGTCTATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAA > NC_000913/303500‑303774 | tttgccagcaacaatgcctctcttctgtagcgttcttcggatctacatcatcctgaagtagcgcgatatactcactgactatctgcatcagaatataaaaagcaatgtttttaacctATAAAAATGGCGCAAATACAGCGCCAttt > 4:106926‑M1/118‑146 (MQ=255) tctggttctggtCTATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCTGTATTTgcgccatttttatcattcaatg < 6:32381‑M2/148‑23 (MQ=255) tctggttctggtCTATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTaaaaa < 7:128092/149‑1 (MQ=255) gctggttctcgtCTATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCAAGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTaaaaa < 7:126269/148‑1 (MQ=255) tATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAgcgccattttta > 1:141655‑M2/1‑113 (MQ=255) tATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAgcgccattttta < 2:141655‑M2/125‑13 (MQ=255) aaaTTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCt > 5:708/1‑149 (MQ=255) attattCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCTGTATTTgcgccatttttatcattcaatgcattatctgtttgagcctaaa > 5:32381‑M2/1‑106 (MQ=255) attattCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCTGTATTTgcgccatttttatcattcaatgcattatctgtttgagccta < 7:46879‑M2/147‑42 (MQ=255) cGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCTGTATTTgcgccatttttatcattcaatgcattatctgtttgagcctaaagggatctcagggtctggc > 8:46879‑M2/1‑88 (MQ=255) gAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGta > 3:60160/1‑149 (MQ=255) aaTTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAgcgccattttt < 7:34396‑M2/91‑12 (MQ=255) aaTTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAgcgccattttt > 8:34396‑M2/1‑80 (MQ=255) cTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTAc < 6:708/148‑1 (MQ=255) cTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACt > 3:28946/1‑149 (MQ=255) aTGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCa < 4:60160/149‑1 (MQ=255) aGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGgatgat < 2:168275/149‑1 (MQ=255) aGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGgatgat > 1:168275/1‑149 (MQ=255) cTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCgaagaa < 7:23144/149‑1 (MQ=255) tCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGt < 2:37096/104‑1 (MQ=255) tCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGt > 1:37096/1‑104 (MQ=255) tCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAg < 8:32665/121‑1 (MQ=255) tCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAg > 7:32665/1‑121 (MQ=255) aaaaagcaatgtttttaacctATAAAAATGGCGCTGTATTTGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGgatga < 3:84202‑M1/103‑1 (MQ=255) aaaaagcaatgtttttaacctATAAAAATGGCGCTGTATTTGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGgatga > 4:84202‑M1/22‑124 (MQ=255) gATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGc > 3:8749/1‑137 (MQ=255) gATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGc < 4:8749/137‑1 (MQ=255) tAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGt < 1:125760/120‑1 (MQ=255) tAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGt > 2:125760/1‑120 (MQ=255) atgtttttaacctATAAAAATGGCGCTGTATTTGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCa < 7:81588‑M1/97‑1 (MQ=255) atgtttttaacctATAAAAATGGCGCTGTATTTGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCa > 8:81588‑M1/14‑110 (MQ=255) ttGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCAttgttg > 8:128979/1‑149 (MQ=255) aTGGCGCAAATACAgcgcaatttttatcattaaatgcattatctgtttgagcctaaagggatctcagggtctatctcatgagtaattctgtttagtaaatatctgagaataatatcaatagaacagaac < 1:174576‑M2/129‑116 (MQ=255) TGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATag < 2:123008‑M1/146‑1 (MQ=255) GCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAgaagaa > 1:123008‑M1/1‑149 (MQ=255) | TGGACATTTAGTCTGGTTCTGGTCTATTGAAATTATTCGCAGAGATTTACGAAACAGAATTACTCATGAGCCAGACCCTGAGATCCCTTTAGGCTCAAACAGATAATGCATTGAATGATAAAAATGGCGCAAATACAGCGCCATTTTTATAGGTTAAAAACATTGCTTTTTATATTCTGATGCAGATAGTCAGTGAGTATATCGCGCTACTTCAGGATGATGTAGATCCGAAGAACGCTACAGAAGAGAGGCATTGTTGCTGGCAAATAGAAGAA > NC_000913/303500‑303774 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |