| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 514,015 | A→C | 19.1% | intergenic (+46/+31) | yhfA → / ← PP_0426 | hypothetical protein/membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 514,015 | 0 | A | C | 19.1% | 31.4 / 6.8 | 21 | intergenic (+46/+31) | yhfA/PP_0426 | hypothetical protein/membrane protein |
| Reads supporting (aligned to +/- strand): ref base A (10/7); new base C (2/2); total (12/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01 | |||||||||||
TATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGC > NC_002947/513892‑514148 | tATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCg > 1:424080/1‑150 (MQ=255) gCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGcc > 2:412929/1‑150 (MQ=255) tCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCgg > 1:390964/1‑150 (MQ=255) gAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGaa < 1:414262/150‑1 (MQ=255) cGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATc > 2:318716/1‑150 (MQ=255) ccATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGa < 2:424080/150‑1 (MQ=255) gAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGcttc > 1:99616/1‑150 (MQ=255) gCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAgcgc > 1:204137/1‑135 (MQ=255) gCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAgcgc < 2:204137/135‑1 (MQ=255) tGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAAACGGCTACCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTccc > 1:88782/1‑150 (MQ=255) tGAGAGTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTccc > 2:149027/1‑150 (MQ=255) ccGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCtt < 2:390964/150‑1 (MQ=255) gCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCGCGGCGCTTTTGGTaaa > 1:422372/1‑150 (MQ=255) gCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTaaa < 2:99616/150‑1 (MQ=255) gCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCTCTCTGGTCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCGCGGCGCTTTTGGTAAACCAgg < 2:422372/150‑1 (MQ=255) gCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAgg > 1:211696/1‑150 (MQ=255) cTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGc < 1:318716/150‑1 (MQ=255) tgggaGCGGGTTCACCCGCTTCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAgc < 1:107738/97‑1 (MQ=255) tgGGAGCGGGTTCACCCGCTTCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAgc > 2:107738/1‑101 (MQ=255) ggaGCGGGTTCACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTg < 1:362670/126‑1 (MQ=255) ggaGCGGGTTCACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTg > 2:362670/3‑128 (MQ=255) | TATTGCTCGGCTTCGATCATGCTCGAGCGTGCCGGGGTTGAGATTACCCATGGGTATGAAATCGTAGAGCTGGGTTAAGGCTGAGATTCCGGGGCCGCTTTGCAGCCCTTGTGGAAGCGGGTGAACCCGCTCCCACCCGCGGCCCGGGCGAGTTTCAGTTGCGCAGGGCCGGAAGGGCCATCCGCAACCATTGCGGAAGACGGCTTCTTAGCGCCTGGCTCAACCGCTCCCGGCGCTTTTGGTAAACCAGGCCAAGC > NC_002947/513892‑514148 |
| 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: ‑ |