| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 3,700,551 | C→A | 19.8% | L328F (TTG→TTT) | PP_3266 ← | hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 3,700,551 | 0 | C | A | 19.8% | 44.3 / 2.9 | 20 | L328F (TTG→TTT) | PP_3266 | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (10/6); new base A (2/2); total (12/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ACATTCAGCGTCAGCACATTGAGCCGGCGTACAGCGGTCACCTTGTCGATGATGCAGCGTGGGGCGGGCAGGGTCCGGTTCACAAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTGGTCGTTGTCGAAGTAGCAGTACACATCGCGTGTCGCACGCCGAGGCGGCGGCCCAGGCACGATCAGGCGAGCATCGTCGGCCTGGCGGCCTTCACTCCATTGCTGG > NC_002947/3700403‑3700692 | acaTTCAGCGTCAGCACATTGAGCCGGCGTACAGCGGTCACCTTGTCGATGATGCAGCGTGGGGCGGGCAGGGTCCGGTTCACAAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTgcagca > 2:504696/1‑150 (MQ=255) cGTCAGCACATTGAGCCGGCGTACAGCGGTCACCTTGTCGATGATGCAGCGTGGGGCGGGCAGGGTCCGGTTCACAAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTg > 1:275829/1‑150 (MQ=255) cGTCAGCACATTGAGCCGGCGTACAGCGGTCACCTTGTCGATGATGCAGCGTGGGGCGGGCAGGGTCCGGTTCACAAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTg < 2:275829/150‑1 (MQ=255) aTTGAGCCGGCGTACAGCGGTCACCTTGTCGATGATGCAGCGTGGGGCGGGCAGGGTCCGGTTCACAAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAg > 1:487841/1‑150 (MQ=255) ttGTCGATGATGCAGCGTGGGGCGGGCAGGGTCCGGTTCACAAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTGGTCGTTGt < 2:487841/150‑1 (MQ=255) tgCAGCGTGGGGCGGGCAGGGTCCGGTTCACAAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGC‑TTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTGGTCGTTGTCGAAGTAGCAg < 2:459025/150‑1 (MQ=255) gggcCGGGTCCGGTTCACAAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTGGTCGTTGTCGAAGTAGCAGTACACATCGCgt < 1:203420/150‑1 (MQ=255) gcAGGGTCCGGTTCACAAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTGGTCGTTGTCGAAGTAGCAGTACACATCGCgtgt > 1:278675/1‑150 (MQ=255) ggTTCACAAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTGGTCGTTGTCGAAGTAGCAGTACACATCGCGTGTCGCACGCCg > 1:294940/1‑150 (MQ=255) acaAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTg < 1:424245/103‑1 (MQ=255) acaAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTg > 2:424245/1‑103 (MQ=255) acaAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTGGTCGTTGTCGAAGTAGCAGTACACATCGCGTGTCGCACGCCGAggcg > 2:468662/1‑150 (MQ=255) cACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTGGTCGTTGTCGAAGTAGCAGTACACATCGCGTGTCGCACGCCGAGGCGGCGGccc > 1:204164/1‑150 (MQ=255) ccTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTGGTCGTTGTCGAAGTAGCAGTACACATCGCGTGTCGCACGCCGAGGCGGCGGCCCAGGCACGATCAGGCGa < 2:294940/150‑1 (MQ=255) ggTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGAAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTg > 1:389040/1‑72 (MQ=255) ggTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGAAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTg > 1:388522/1‑72 (MQ=255) ggTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGAAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTg < 2:388522/72‑1 (MQ=255) ggTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGAAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTg < 2:389040/72‑1 (MQ=255) aTGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTGGTCGTTGTCGAAGTAGCAGTACACATCGCGTGTCGCACGCCGAGGCGGCGGCCCAGGCACGATCAGGCGAGCATCGTCGGCCTGGCGGcc > 1:328863/1‑150 (MQ=255) gTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTAGTAGGGTGTGTGTACCTTGTGGTCGTTGTCGAAGTAGCAGTACACATCGCGTGTCGCACGCCGAGGCGGCGGCCCAGGCACGATCAGGCGAGCATCGTCGGCCTGGCGGCCTTc > 1:342775/1‑150 (MQ=255) ttttGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTGGTCGTTGTCGAAGTAGCAGTACACATCGCGTGTCGCACGCCGAGGCGGCGGCCCAGGCACGATCAGGCGAGCATCGTCGGCCTGGCGGCCTTCACTCCATTGCTgg > 1:339895/1‑150 (MQ=255) | ACATTCAGCGTCAGCACATTGAGCCGGCGTACAGCGGTCACCTTGTCGATGATGCAGCGTGGGGCGGGCAGGGTCCGGTTCACAAGGCCACCTCCGGCACCACGCCTGGCTCGGTCATCAGCTCATGGTCCAACGACAGCTTTTGCAGCAAGCGCCTGGCGTCGTAGGGTGCGTGAACCTTCTGGTCGTTGTCGAAGTAGCAGTACACATCGCGTGTCGCACGCCGAGGCGGCGGCCCAGGCACGATCAGGCGAGCATCGTCGGCCTGGCGGCCTTCACTCCATTGCTGG > NC_002947/3700403‑3700692 |
| 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: ‑ |