| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,634,023 | 0 | A | T | 59.5% | 37.9 / 44.4 | 37 | G123G (GGT→GGA) | nohA | NohA |
| Reads supporting (aligned to +/- strand): ref base A (6/9); new base T (8/14); total (14/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATC‑CCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTG > NC_000913/1633948‑1634115 | aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg > 1:1703294/4‑93 (MQ=0) tccgTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATC‑CCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGaa > 1:2723532/4‑93 (MQ=1) agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt > 1:636218/4‑93 (MQ=0) agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt < 1:6308474/90‑1 (MQ=0) agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt > 1:1936998/4‑93 (MQ=0) agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt > 2:1522747/4‑93 (MQ=0) gggaaaCGCCGCTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGtt > 2:3752460/4‑93 (MQ=0) gaagcgccgcTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc < 2:356664/89‑1 (MQ=0) gaagcgccgcTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc < 1:5087865/89‑1 (MQ=0) gaagcgccgcTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc > 2:5709120/5‑93 (MQ=0) gaagcgccgcTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc < 2:3699466/89‑1 (MQ=0) gaagcgccgcTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTGCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc < 2:527712/89‑1 (MQ=0) aagcgccgcTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc < 2:2768562/90‑1 (MQ=0) aagcgccgcTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc < 1:1922894/90‑1 (MQ=0) aagcgccgcTGAACCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTGCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc < 1:2167919/90‑1 (MQ=0) cccacTGCACCGACAGGGGGATC‑CCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTccacca < 1:3863793/89‑1 (MQ=1) tgaaCCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCa < 2:5138648/90‑1 (MQ=0) tgaaCCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCa < 2:429946/90‑1 (MQ=0) tgaaCCGACAGAGGAACT‑CCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCa > 2:3817845/4‑93 (MQ=0) tgaaCCGACA‑GAGGATCTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCa < 1:5171265/90‑1 (MQ=0) tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt < 2:6089771/90‑2 (MQ=0) tcaagaATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCAtttttt < 2:3936211/90‑2 (MQ=0) atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc < 2:604941/90‑1 (MQ=0) atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc > 2:2976294/4‑93 (MQ=0) gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg > 1:5924767/1‑93 (MQ=31) gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg > 2:3638288/1‑93 (MQ=31) cAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGt < 1:6360284/93‑1 (MQ=32) cAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGt < 2:5130681/93‑1 (MQ=32) aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc < 2:4984056/93‑1 (MQ=33) aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc < 1:1540697/93‑1 (MQ=33) aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg < 2:1390536/93‑1 (MQ=33) aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg < 1:542489/93‑1 (MQ=33) tttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCGGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCgg > 1:3636526/1‑93 (MQ=18) tttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCgg < 2:5119015/93‑1 (MQ=33) ttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGc < 1:2930288/93‑1 (MQ=33) tCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGcc > 1:5425564/1‑93 (MQ=33) tCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGcc < 2:5704645/93‑1 (MQ=34) tCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGcc < 1:204833/93‑1 (MQ=33) cACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCt > 2:2832563/1‑93 (MQ=33) aCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTg > 2:2130831/1‑93 (MQ=33) aCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTg > 2:4600362/1‑93 (MQ=33) aCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTg < 1:4420383/93‑1 (MQ=33) aCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTg > 1:3537400/1‑93 (MQ=33) aCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTg < 1:2549570/93‑1 (MQ=33) aCCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTg > 1:1206249/1‑93 (MQ=33) | AACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATC‑CCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGGCCTG > NC_000913/1633948‑1634115 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |