| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,812,386 | 0 | C | A | 27.3% | 43.5 / 11.0 | 22 | L58F (TTG→TTT) | yicR | RadC‑like JAB domain‑containing protein YicR |
| Reads supporting (aligned to +/- strand): ref base C (16/0); new base A (0/6); total (16/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.34e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCG > NC_000913/3812322‑3812449 | aCTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAAc > 1:1216027/1‑65 (MQ=255)aCTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAAc > 1:1421009/1‑65 (MQ=255)aCTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAAc > 1:1582250/1‑65 (MQ=255)aCTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAAc > 1:638425/1‑65 (MQ=255)aCTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAAc > 1:1725105/1‑65 (MQ=255)aCTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAAc > 1:4190/1‑65 (MQ=255) gCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGcc > 1:3102135/1‑65 (MQ=255) ccaTTTCCATGACCCCCAAAAAATTGTCCATATTCAGAGTTAAAAAAGCCATAAAGAGAGCCGaa < 1:3090128/65‑1 (MQ=255) ccaATCCCATGAACCCCAATAGTTTGTTAATATTCAGAGGCAAAAAAGCCATAAAGAGAGCCGaa < 1:3093010/65‑1 (MQ=255) tCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATtct > 1:932606/1‑65 (MQ=255) tCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATtct > 1:2578154/1‑65 (MQ=255) gCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCAtct > 1:1638616/1‑65 (MQ=255) ccATTAATTTGTTAATTTGCAGGGGAAAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCAtctc < 1:2292460/65‑1 (MQ=255) cacaAAATTGTTAATATTCAGAGTAAAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATctct < 1:975587/65‑1 (MQ=255) aaTTGTTCATATTCAGAGGTAAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGc < 1:1671840/65‑1 (MQ=255) aTTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGcc > 1:1263298/1‑65 (MQ=255) aTTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGcc > 1:2782386/1‑65 (MQ=255) aTTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGcc > 1:1819102/1‑65 (MQ=255) tGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAg > 1:1824683/1‑65 (MQ=255) ttaatatTAAGAGATAAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAggg < 1:745789/62‑1 (MQ=255) ttaatatTAAGAGATAAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAggg < 1:2707766/62‑1 (MQ=255) ttCATATTAAGAGTTAAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAggg < 1:191501/65‑1 (MQ=255) ttCATATTAAGAGTTAAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAggg < 1:1741756/65‑1 (MQ=255) tatatatGCAGAGTAAAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAggg < 1:2034799/62‑1 (MQ=255) atCATATTCAGAGGAAAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAggg < 1:2230707/64‑1 (MQ=255) tCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGt > 1:2100316/1‑65 (MQ=255) tCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGt > 1:785434/1‑65 (MQ=255) tCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGt > 1:777430/1‑65 (MQ=255) ttCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATa < 1:2984655/65‑1 (MQ=255) ttCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATa > 1:2193539/1‑65 (MQ=255) cAGAGGCAAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACa < 1:20414/65‑1 (MQ=255) aagagTTAAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACa < 1:2851627/64‑1 (MQ=255) aagagTTAAAAAAGCCAAAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACa < 1:2131694/64‑1 (MQ=255) gtggTAAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATc < 1:210201/63‑1 (MQ=255) ggggTTAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATc < 1:3233575/63‑1 (MQ=255) gagTTTAAAAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATc < 1:175149/65‑1 (MQ=255) ggTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCtt > 1:2760264/1‑65 (MQ=255) ggTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCtt > 1:522261/1‑65 (MQ=255) ggTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCtt > 1:2309631/1‑65 (MQ=255) aCAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTAcag > 1:1554622/1‑63 (MQ=255) aCAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACcg > 1:639294/1‑65 (MQ=255) | ACTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTCATATTCAGAGGTTAACAAGCCATAAAGAGAGCCGAAATTCTCCAGCATCTCTTTTGCCAGGGTTAATACATCTTTACCG > NC_000913/3812322‑3812449 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |