| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,888,383 | 0 | A | G | 20.4% | 106.8 / 24.4 | 56 | intergenic (‑71/+4) | cysH/cysI | phosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent/sulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding |
| Reads supporting (aligned to +/- strand): ref base A (28/15); new base G (8/3); total (37/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.33e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01 | |||||||||||
CAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTT > NC_000913/2888286‑2888481 | cAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATg > 2:648716/1‑100 (MQ=255) tttAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATc > 2:157845/1‑100 (MQ=255) tttAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATc < 1:80533/100‑1 (MQ=255) tttAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTACTTATCAGGCCTACAAGACCGGGCTGATGGTTAATc < 1:920632/100‑1 (MQ=255) tAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGGTGGTTAATcc < 2:519601/99‑1 (MQ=255) tAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATccc < 2:1028515/100‑1 (MQ=255) aTCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCcaca < 1:466263/100‑1 (MQ=255) tCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGATCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACaa > 1:450476/1‑100 (MQ=255) cGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACaaa > 2:1284617/1‑100 (MQ=255) gAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAAt > 1:664370/1‑100 (MQ=255) aGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATc < 2:1145103/100‑1 (MQ=255) tGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCAcgc > 2:874256/1‑100 (MQ=255) ggACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGGTGGTTAATCCCACAAATCAcgcg > 1:684261/1‑100 (MQ=255) ggACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCAcgcg > 2:1024393/1‑100 (MQ=255) gACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACgcgc > 2:106169/1‑100 (MQ=255) aTAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCgg > 1:727355/1‑100 (MQ=255) tAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGGTGGTTAATCCCACAAATCACGCGCCGGa < 1:990596/100‑1 (MQ=255) tAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGa > 1:236938/1‑100 (MQ=255) aGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGAt > 1:357135/1‑100 (MQ=255) cTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCa < 1:326089/100‑1 (MQ=255) cTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCa > 1:532806/1‑100 (MQ=255) cACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCAcc > 2:921946/1‑100 (MQ=255) ntttGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGGTGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGc > 2:398124/2‑100 (MQ=255) gTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGCTGGTTAATCCCCCCAATCAAGAGCCGGATAAAGCaaaggg > 1:433898/1‑94 (MQ=25) gTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGc > 2:621281/1‑100 (MQ=255) ttgccttgccTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGa > 1:1214520/1‑100 (MQ=255) cttgcctGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCACACAAATCACGCGCCGGATCGAGCACCGGGCGAatga > 2:878042/1‑100 (MQ=255) ttgcctGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAatgat < 2:499797/100‑1 (MQ=255) tgcctGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGGTGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAatgatg > 1:1062880/1‑100 (MQ=255) cGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCcgcacgc > 1:517380/1‑100 (MQ=255) gACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGGTGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCcgcacgca > 1:933150/1‑100 (MQ=255) cGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCAcc < 1:157847/100‑1 (MQ=255) cGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCAcc < 1:106171/100‑1 (MQ=255) gCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGGTGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCg < 1:1048393/100‑1 (MQ=255) cTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGt < 1:938604/100‑1 (MQ=255) cACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTaa > 1:1130298/1‑100 (MQ=255) cgTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTc > 2:632579/1‑100 (MQ=255) cgTTCTTATCAGGCCTACAAGACCGGCCTGCCGGTTAATCCCACACACCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTc < 2:451883/100‑1 (MQ=255) cTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCg > 2:237055/1‑100 (MQ=255) tCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGc < 2:1038528/100‑1 (MQ=255) ggCCTACAAGACCGGGCTGGTGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCtt > 1:628366/1‑100 (MQ=255) ggCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCtt > 2:170489/1‑100 (MQ=255) ggCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCtt > 1:578770/1‑100 (MQ=255) gCCTACAAGACCGGGCTGGTGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTc > 1:488926/1‑100 (MQ=255) gCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTc < 2:1117830/100‑1 (MQ=255) ccTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCa > 2:524216/1‑100 (MQ=255) cAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCg > 2:437165/1‑100 (MQ=255) cAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCAACCg > 1:1008538/1‑100 (MQ=255) aaGACCGGGCTGGTGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGc > 2:338732/1‑100 (MQ=255) aGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCt > 1:1236160/1‑100 (MQ=255) ccGGGCTGGTGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCTTcg > 1:609062/1‑100 (MQ=255) ccGGGCTGATGGTTAATCCAACAAATCACGCGCCGGCTCGAGCTCCGGGTAAATGATGATCGCACGATTCTTAAAGTCGTCGAATCCTTCACTCGCTTcg < 2:433894/100‑1 (MQ=255) cGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCAACCGATTcg > 2:646962/1‑99 (MQ=255) ggCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCTTCgcgc < 1:648720/100‑1 (MQ=255) gCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGct < 1:80418/100‑1 (MQ=255) gATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCtt > 2:1142151/1‑100 (MQ=255) | CAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAAGTCGCCGAAGCCTTCACCCGCTTCGCGCTCTT > NC_000913/2888286‑2888481 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |