| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 3,675,925 | C→A | 8.0% | A18A (GCC→GCA) | PP_3240 → | hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 3,675,925 | 0 | C | A | 8.0% | 158.1 / 2.5 | 50 | A18A (GCC→GCA) | PP_3240 | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (25/21); new base A (2/2); total (27/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01 | |||||||||||
GCAGTGTCGGGCGCAATTGGCCTCGGTTGCAGATTAATCAAACTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAGTAAACCACGGCCTTGATTACCGGTCTCGTCGCCGGCAAGC > NC_002947/3675780‑3676073 | gcAGTGTCGGGCGCAATTGGCCTCGGTTGCAGATTAATCAAACTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATa > 1:841611/1‑150 (MQ=255)gcAGTGTCGGGCGCAATTGGCCTCGGTTGCAGATTAATCAAACTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATa > 2:725015/1‑150 (MQ=255) gtCGGGCGCAATTGGCCTCGGTTGCAGATTAATCAAACTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAg < 2:565899/150‑1 (MQ=255) tCGGGCGCAATTGGCCTCGGTTGCAGATTAATCAAACTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAgg > 1:592472/1‑150 (MQ=255) tCGGGCGCAATTGGCCTCGGTTGCAGATTAATCAAACTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAgg < 1:355777/150‑1 (MQ=255) gcgcAATTGGCCTCGGTTGCAGATTAATCAAACTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAGCTCAGTGGGCCG‑AAGAATATCAGCCACCACACAAAGAGGGGCCCCACTTCGCGTAGCAGGCCGATAAAAAGGCCGa > 1:928577/1‑150 (MQ=255) cgcAATTGGCCTCGGTTGCAGATTAATCAAACTTTTCCTGGGGAGTGGTTTCGTCCTTGGACATTTATTAGAG‑ACTACGCCATGCTTGGATCTCAACCACCACACAAAGAGGAGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGAc < 2:401719/150‑1 (MQ=255) tCGGTTGCAGATTAATCAAACTTTTTCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGa < 2:183562/150‑1 (MQ=255) tCGGTTGCAGATTAATCAAACTTTTCCTGGCGAATCGTTTCTGCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGa < 2:1001205/150‑1 (MQ=255) aTCAAACTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAAGCCGACCTTGGTTTCGACCCCGACTCCCCgg > 1:1034888/1‑150 (MQ=255) tCAAACTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGa > 2:485831/1‑150 (MQ=255) tCAAACTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGa < 2:841611/150‑1 (MQ=255) aaaCTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATg > 2:851290/1‑150 (MQ=255) aaCTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGc > 1:212012/1‑150 (MQ=255) cTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGAcc > 1:710542/1‑150 (MQ=255) gCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTc > 1:404793/1‑150 (MQ=255) gTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGa > 2:1000219/1‑150 (MQ=255) gTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGa > 1:127821/1‑150 (MQ=255) tcttcCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGAccc < 1:1000219/150‑1 (MQ=255) ttAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAggg > 2:925913/1‑150 (MQ=255) gACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGggcc < 2:710542/150‑1 (MQ=255) cATCGGAG‑ACAACGCCATGCCTGGAGCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCaa < 1:485831/150‑1 (MQ=255) tCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGa < 2:220282/100‑1 (MQ=255) tCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGa > 1:220282/1‑100 (MQ=255) tCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCaa > 2:1117578/1‑148 (MQ=255) tCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCaa < 1:1117578/148‑1 (MQ=255) aCAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTa < 1:925913/150‑1 (MQ=255) aCAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTa > 1:175159/1‑150 (MQ=255) gCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGAc > 2:851307/1‑150 (MQ=255) tGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTg < 1:725015/150‑1 (MQ=255) tcAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGa < 2:212012/150‑1 (MQ=255) cAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGAcc < 2:227379/94‑1 (MQ=255) cAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGAcc > 1:227379/1‑94 (MQ=255) aaCCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAg < 2:1042961/149‑1 (MQ=255) aaCCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAg > 1:1042961/1‑149 (MQ=255) accacACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTa > 1:1153731/1‑125 (MQ=255) accacACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTa < 2:1153731/125‑1 (MQ=255) accacACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAGTAAAc < 2:928577/150‑1 (MQ=255) cacaAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAGTAAACCACg < 1:851290/150‑1 (MQ=255) aaaGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGa < 2:397328/92‑1 (MQ=255) aaaGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGa > 1:397328/1‑92 (MQ=255) aaaGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAGTAAACCACGGcc < 1:851307/150‑1 (MQ=255) gCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAGTAAACCACGGCCTTGATTAc > 2:1175581/1‑150 (MQ=255) tCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAGCAAACCACGGCCTTGATTACCgg > 2:1058594/1‑150 (MQ=255) ggCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGCCCAAACCACGGCCTTGATTACCGGTCtcgt > 2:1056807/1‑150 (MQ=255) cGAGCAGGCAGATAAAAAGGCCGACCTTGGTTTCGAcc > 2:1190064/1‑38 (MQ=255) cGAGCAGGCAGATAAAAAGGCCGACCTTGGTTTCGAcc < 1:1190064/38‑1 (MQ=255) aGGCAGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGAc < 2:232780/98‑1 (MQ=255) aGGCAGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGAc > 1:232780/1‑98 (MQ=255) gCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAGTAAACCACGGCCTTGATTACCGGTCTCGTCGCCGGCAAg > 1:413280/1‑150 (MQ=255) ccGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAGTAAACCACGGCCTTGATTACCGGTCTCGTCGCCGGCAAGc < 2:404793/150‑1 (MQ=255) cGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAGTAAACCACg > 1:1160867/1‑118 (MQ=255) cGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAGTAAACCACg < 2:1160867/118‑1 (MQ=255) | GCAGTGTCGGGCGCAATTGGCCTCGGTTGCAGATTAATCAAACTTTTCCTGGCGAGTGGTTTCTTCCTTAGACACTCATCGGAG‑ACAACGCCATGCCTGGATCTCAACCACCACACAAAGAGGGGCCTCACTCGGCCGAGCAGGCCGATAAAAAGGCCGACCTTGGTTTCGACCCCGACTCCCCGGATGCCGAAGACCCTCAAGTCGACCCACAGGGGCCGGCCAAGCCGCCTAGAGACATTGAGAAAAAGGAGTAAACCACGGCCTTGATTACCGGTCTCGTCGCCGGCAAGC > NC_002947/3675780‑3676073 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |