| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,766,331 | A→G | 20.7% | K110E (AAG→GAG) | PP_1577 → | lambda family tail tape meausure protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,766,331 | 0 | A | G | 20.7% | 55.8 / 9.6 | 29 | K110E (AAG→GAG) | PP_1577 | lambda family tail tape meausure protein |
| Reads supporting (aligned to +/- strand): ref base A (11/12); new base G (3/3); total (14/15) | |||||||||||
| 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 | |||||||||||
AATGACAGTTTCCGCCATCAATTCTGCCTCGGAAATAACCAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCGGTGTCACGGCGGATCAGTTCCGTAACCTAAGTGGAAGCCAGGCCCTAGGGCTCTAC > NC_002947/1766189‑1766474 | aaTGACAGTTTCCGCCATCAATTCTGCCTCGGAAATAACCAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGAtgt < 2:640621/150‑1 (MQ=255) aGTTTCCGCCATCAATTCTGCCTCGGAAATAACCAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTTGAGGAtgtcagcaa < 2:547839/150‑7 (MQ=255) gCCATCAATTCTGCCTCGGAAATAACCAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTTGAGGAtgtc < 1:526827/138‑2 (MQ=255) gCCATCAATTCTGCCCCGGAAATAACCAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTTGAGGAtgtc > 2:526827/1‑137 (MQ=255) aTCAATTCTGCCTCGGAAATAACCAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTTGAGGAtgtca < 1:1096812/136‑3 (MQ=255) aTCAATTCTGCCTCGGAAATAACCAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTTGAGGAtgtca > 2:1096812/1‑134 (MQ=255) cAATTCTGCCTCGGAAATAACCAATCTAGCGTCGGTTGCAAACGTCAGAACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCaa < 2:347034/126‑1 (MQ=255) cAATTCTGCCTCGGAAAGAACCAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCaa > 1:347034/1‑126 (MQ=255) tCGGAAATAACCAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAac > 1:767687/1‑150 (MQ=255) tCGGAAATAACCAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACAACAAGGTAGGTGATTTCCTAAac > 1:411891/1‑150 (MQ=255) aaTAACCAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACcgg > 2:181920/1‑150 (MQ=255) ccAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAacac > 2:153577/1‑142 (MQ=255) ccAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAacac < 1:153577/142‑1 (MQ=255) cTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCggg < 1:181920/150‑1 (MQ=255) gCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATg > 2:731556/1‑150 (MQ=255) cGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGa < 2:767687/150‑1 (MQ=255) aGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGATGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAg < 1:190467/150‑1 (MQ=255) tACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAg > 1:752259/1‑76 (MQ=255) tACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAg < 2:752259/76‑1 (MQ=255) aCTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTa > 1:879632/1‑150 (MQ=255) aGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCg < 1:898740/150‑1 (MQ=255) aaGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCGgt < 1:1071634/150‑1 (MQ=255) ggggCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCGGTGTCACGGCGGAt < 1:731556/150‑1 (MQ=255) ccAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCGGTGTCACGGCGGATCAGt > 2:486716/1‑150 (MQ=255) gACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCGGTGTCACGGCGGATCAGTTCCg < 1:486716/150‑1 (MQ=255) gACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCGGTGTCACGGCGGATCAGTTCCGTAACCTAAGTGGaa < 2:759809/145‑1 (MQ=255) gACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCGGTGTCACGGCGGATCAGTTCCGTAACCTAAGTGGaa > 1:759809/1‑145 (MQ=255) gCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCGGTGTCACGGCGGATCAGTTCCGTAAc > 1:180039/1‑130 (MQ=255) gCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCGGTGTCACGGCGGATCAGTTCCGTAAc < 2:180039/130‑1 (MQ=255) ttcaCATCCTTGAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCGGTGTCACGGCGGATCAGTTCCGTAACCTAAGTGGAAGCCAGGCCCTAGGGc > 1:160110/4‑150 (MQ=255) aTCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCGGTGTCACGGCGGATCAGTTCCGTAACCTAAGTGGAAGCCAGGCCCTAGGGCTCTAc > 1:1207590/1‑150 (MQ=255) | AATGACAGTTTCCGCCATCAATTCTGCCTCGGAAATAACCAATCTAGCGTCGGTTGCAAACGTCAGCACTACAGACTTTCAGAAGCTAGCAGTTGGGGCCAAGACAGTAGGCGTAGAGCAAGACAAGCTTGCTGACATCCTCAAGGATGTGAACGACAAGGTAGGTGATTTCCTAAACACCGGCGGCGGCGGGATGGCTGACTTCTTCGAGCAGGTCGCGCCTAAGGTCGGTGTCACGGCGGATCAGTTCCGTAACCTAAGTGGAAGCCAGGCCCTAGGGCTCTAC > NC_002947/1766189‑1766474 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |