| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 516,893 | 0 | G | C | 32.1% | 51.6 / 12.0 | 28 | G104A (GGA→GCA) | fetB | iron export ABC transporter permease; peroxide resistance protein |
| Reads supporting (aligned to +/- strand): ref base G (18/1); new base C (0/9); total (18/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.45e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGT > NC_000913/516831‑516950 | aCGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGCAATTa < 1:24947/68‑1 (MQ=255) gCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTAGGGGCGGGAATTAca > 2:974018/1‑67 (MQ=255) tatatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACACTGGCGGt > 2:396734/1‑68 (MQ=255) tatatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGCAATTACCCTGGCGGt < 2:681934/68‑1 (MQ=255) tatatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGCAATTACCCTGGCGGt < 2:1029447/68‑1 (MQ=255) tatTGCTAAAGCTTTTATCTCATCGTTTATTGCCTTTAAGGTCGGGGCGGGAATTACACTGGCGGTga > 2:521930/1‑67 (MQ=255) tatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGc > 2:396496/1‑68 (MQ=255) tatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGc > 2:205891/1‑68 (MQ=255) ttGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTg > 1:595095/1‑68 (MQ=255) tGCTACAGCTTTTATCTCATCGTTTATTGCCATTACGGGCGGGGCGGCAATTACCCTGGCGGTGCTGa < 2:118209/68‑1 (MQ=255) tGCTAATGCTTTTATCTCATCGTTTATTCCCATAACGGTCGGGGCGGCAATTACCCTGGCGGTGCTGa < 2:293910/68‑1 (MQ=255) gcaaaaGCTATTATCTAGTCGATTATTGCCATTACGGTCGGGGCGGGAATTACACTGGCGGTGCTGAt > 2:553327/4‑68 (MQ=255) gCTAAAGCTTTTATCTCGTCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGAt > 1:482243/1‑68 (MQ=255) gCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGAt > 2:1042465/1‑68 (MQ=255) gCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGAt > 1:621444/1‑68 (MQ=255) gCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTAACCTGGCGGTGCTGAt > 2:624866/1‑68 (MQ=255) aaGCTGTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGCAATTACCCTGGCGGTGCTGATtctc < 2:324413/68‑1 (MQ=255) gCTTTTATCTCATCGTATATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATtctctc > 2:511251/1‑68 (MQ=255) cTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCa > 2:968921/1‑68 (MQ=255) tctcATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCg > 1:250971/1‑68 (MQ=255) tctcATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCg < 1:359360/68‑1 (MQ=255) cctcATCTTTTTTTGCCTTTCCGGTGGGGGCGCCAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCg < 2:91386/67‑1 (MQ=255) cGTTTATTGCCATTACGGTCGGGGCGGCAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGaa < 1:836647/68‑1 (MQ=255) tATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTa > 1:630800/1‑68 (MQ=255) ccATTACGGTCGTGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCg > 2:506155/1‑68 (MQ=255) ccATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCg > 2:310498/1‑68 (MQ=255) agggTCGGGGCGGCAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCa < 2:227773/66‑1 (MQ=255) gTCGGGGCGGGAATTACACTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGt > 2:814172/1‑68 (MQ=255) | ACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGT > NC_000913/516831‑516950 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |