| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 3,866,243 | C→A | 43.5% | E227* (GAG→TAG) | PP_3415 ← | LacI family transcriptional regulator |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 3,866,243 | 0 | C | A | 43.5% | 10.8 / 39.8 | 37 | E227* (GAG→TAG) | PP_3415 | LacI family transcriptional regulator |
| Reads supporting (aligned to +/- strand): ref base C (11/10); new base A (10/6); total (21/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.39e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01 | |||||||||||
TTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTG > NC_002947/3866115‑3866389 | cccttgacgcaaggcctgcagcaccgcgccctgagccaggccgtcACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCTGGGCCAATACCTCCAGCTc < 1:315041‑M2/105‑1 (MQ=255) cttgacgcaaggcctgcagcaccgcgccctgagccaggtcgtcACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCgg < 1:361299‑M2/87‑1 (MQ=255) cttgacgcaaggcctgcagcaccgcgccctgagccaggtcgtcACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCgg > 2:361299‑M2/44‑130 (MQ=255) gCAAGGCCTGCAGCACCGCGCCCTGAGCACGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTACGGGGCCAATACCTCCAGCTCCGCAg > 1:240372/1‑150 (MQ=255) gcaaggcctgcagcaccgcgccctgagccaggtcgtcACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCt > 2:403876‑M2/38‑150 (MQ=255) ccGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTACGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCct > 1:635256/1‑150 (MQ=255) ccTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCACTCGCGAATGGCTACGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCa > 1:8972/1‑150 (MQ=255) tcacccACGTCGTCTTTACAGAAGAAGATACCTTCGACATCAGGCGCCTGCTGCATCATCTGGCTGAACAGCACACCCCCCAGGTCAATCGACGATCTCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCaaa < 1:406376/146‑1 (MQ=255) gaagaagaTACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGc > 2:593048/1‑150 (MQ=255) aTACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGc < 2:695253/150‑1 (MQ=255) tACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTACGGGGCCAATACCt > 2:618946/1‑87 (MQ=255) tACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTACGGGGCCAATACCt < 1:618946/87‑1 (MQ=255) tACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTACGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCa < 2:13723/150‑1 (MQ=255) tACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTACGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCa > 1:509851/1‑150 (MQ=255) aCCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCAt < 2:319107/150‑1 (MQ=255) ccGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTACGGGGCCAATACCTCCAGCTCCGCAGCCTGc > 1:479467/1‑102 (MQ=255) ccGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTACGGGGCCAATACCTCCAGCTCCGCAGCCTGc < 2:479467/102‑1 (MQ=255) cGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTACGGGGCCAATAACTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCa > 1:755889/1‑150 (MQ=255) gACATCCGGCGCCTGCTGCCTCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCt < 2:123145/150‑1 (MQ=255) gACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGACAACGCCTGACGGAAACCCTCGGCACGCTGCATCACCCt > 2:707570/1‑150 (MQ=255) aTCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTACGGGGCCAATAACTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCgg < 2:755889/150‑1 (MQ=255) ctgcATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTgcgc > 2:590912/1‑150 (MQ=255) gctggctgAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCt > 2:156019/1‑96 (MQ=255) gctggctgAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCt < 1:156019/96‑1 (MQ=255) cacaCTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCt < 2:592025/95‑1 (MQ=255) cacaCTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCt > 1:592025/1‑95 (MQ=255) cacaCTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGcgcg < 1:804257/130‑1 (MQ=255) cacaCTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGcgcg > 2:804257/1‑130 (MQ=255) cacaCTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCtg < 1:590912/150‑2 (MQ=255) cccAGGGCAATCGACGATGGCTACGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCcg > 1:458317/1‑150 (MQ=255) aGGGCAATCGACGATGGCTACGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAATGCGCTTGCGCCCgcgc < 2:635256/150‑1 (MQ=255) gggCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCt > 1:275940/1‑150 (MQ=255) aaTCGACGATGGCTACGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGAt > 2:696029/1‑150 (MQ=255) tCGACGATGGCTACGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCa > 2:578545/1‑150 (MQ=255) acgaTGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATGAgcg > 1:266403/1‑148 (MQ=255) aTGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGcc > 2:376819/1‑150 (MQ=255) gCTACGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCCGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTg < 1:347416/150‑1 (MQ=255) | TTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGCTGCATCAGCTGGCTGAACAGCACACTGCCCAGGGCAATCGACGATGGCTCCGGGGCCAATACCTCCAGCTCCGCAGCCTGCAAACCAGCCTCTGCCAAAGCCTGACGGAAACCCTCGGCACGCTGCATCACCCTCGGGTCCAGCTGCGCGGCAATGAATGCAAGGCGCTTGCGCCCGCGCTCGATCAGGTGCCGTG > NC_002947/3866115‑3866389 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |