| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,634,023 | 0 | A | T | 70.0% | 10.7 / 14.0 | 20 | G123G (GGT→GGA) | nohA | NohA |
| Reads supporting (aligned to +/- strand): ref base A (2/4); new base T (8/6); total (10/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.28e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCC > NC_000913/1633948‑1634113 | aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCTCGACAg > 1:228532/4‑93 (MQ=0)aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg > 1:2004858/4‑93 (MQ=0)aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg > 2:383135/4‑93 (MQ=0) gtcccGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATg < 1:1652262/90‑1 (MQ=1) gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc > 1:3756090/5‑93 (MQ=0) gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc < 1:3787623/89‑1 (MQ=0) gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc > 1:1188043/5‑93 (MQ=0) gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTATCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc > 1:1784265/5‑93 (MQ=0) aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc < 1:3525661/90‑1 (MQ=0) tgaaCCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCa < 2:425980/90‑1 (MQ=0) tgaaCCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCa > 2:4154656/4‑93 (MQ=0) atcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAttttt > 2:4206311/5‑93 (MQ=0) tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt < 2:2880813/90‑2 (MQ=0) atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc < 1:4560970/90‑1 (MQ=0) atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc < 1:3367950/90‑1 (MQ=0) gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg < 2:2080269/93‑1 (MQ=31) cAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGt < 2:1292979/93‑1 (MQ=32) cAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGt < 2:1181650/93‑1 (MQ=32) tttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCgg > 2:2862938/1‑93 (MQ=33) tttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCgg > 1:3623326/1‑93 (MQ=33) tCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGcc < 2:907627/93‑1 (MQ=33) | AACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCC > NC_000913/1633948‑1634113 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |