| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,311,359 | T→A | 20.2% | F12Y (TTT→TAT) ‡ | PP_3784 → | hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,311,359 | 0 | T | A | 20.2% | 106.1 / 21.1 | 50 | F12Y (TTT→TAT) ‡ | PP_3784 | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (20/19); new base A (6/4); total (27/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.31e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGA > NC_002947/4311213‑4311505 | acaACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGc < 2:255003/150‑1 (MQ=255) aaCCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGt > 1:1268337/1‑150 (MQ=255) cgcgTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCAAATGCGTTGc < 1:853212/150‑1 (MQ=255) gCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGc > 2:594124/1‑150 (MQ=255) aCCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCtctc > 2:1196675/1‑150 (MQ=255) aCCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCtctc < 1:594124/150‑1 (MQ=255) aCTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGtc < 2:1087909/150‑1 (MQ=255) aCTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGtc > 1:1087909/1‑150 (MQ=255) gCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATtcatca < 1:1405725/150‑1 (MQ=255) aCCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCAc < 1:463727/150‑1 (MQ=255) aCCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCAATCTCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCAc < 1:513349/150‑1 (MQ=255) gAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACgtg > 1:296107/1‑150 (MQ=255) cGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGcc > 1:1041604/1‑150 (MQ=255) cGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCAAATGCGTTg < 1:1311667/98‑1 (MQ=255) cGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCAAATGCGTTg > 2:1311667/1‑98 (MQ=255) cGCTCTGGCCCCCAAATCGATTTCGCGCCCTCACTTCATTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGTTG‑‑GCA‑ACTCATTTGCCTTGCCCTTCTCCCGCTCTCATCATCTATTCATCAGTGCCACGTGTACGCCTcaca < 2:650835/150‑1 (MQ=255) tctGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACAt > 1:119622/1‑150 (MQ=255) cccccGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGc > 1:111816/1‑150 (MQ=255) cccGAATCGATTGCGCGCACTCCCGTCAGTGAAAGGAAGATGATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAGTTCATCATTGCCACGTGTACGCCTGACACATTCAGCGa < 2:294670/150‑1 (MQ=255) gAATCGATTGCGCGCCCTCCCGTCAGTGCAAGGCTGATCATCAAGGTCTAACATCATTCCATCCGGGGGCGCAGACTCCTGTG‑GTCGCGGGTGACC‑‑CTCGCGGCATCAATTCATCAGTGCCACGTGTACGCCTGACTCATGCAGCGACTg > 1:1391700/1‑150 (MQ=255) ggcgcCCTCCCGTCAGTGAAAGGAAGATCATCAATCTCTTTCAGCATTCACTCCGGGA‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCaa < 2:337727/149‑1 (MQ=255) gcgcCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAACATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAAt > 1:600255/1‑150 (MQ=255) gcCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAACATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCa < 2:600255/150‑1 (MQ=255) cTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTg < 1:177187/150‑1 (MQ=255) tCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTgc > 1:881633/1‑150 (MQ=255) tGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCAAATGCGTTGc < 2:1316961/60‑1 (MQ=38) tGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCAAATGCGTTGc > 1:1316961/1‑60 (MQ=38) aGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGaa < 2:1041604/150‑1 (MQ=255) atcatcAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACa > 2:1038682/1‑150 (MQ=255) atcaccAATGGCTAACAGCATTCCCTCCCGGG‑‑GCA‑ACGCATTTGCGTTGCACTTCGCCGGCTCTCGTCATCAATTCATCAGTGCAACGTGTACGCCTGCCGCATTCCGAGCCTGCCTGTCGCAATCACTGACGCCATCCAGCAACGAACa < 2:558318/150‑1 (MQ=255) tcatcaATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAg > 1:140876/1‑150 (MQ=255) catcaATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGt > 2:668351/1‑80 (MQ=255) catcaATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGt < 1:668351/80‑1 (MQ=255) tcaATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGa < 2:567501/95‑1 (MQ=255) tcaATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGa > 1:567501/1‑95 (MQ=255) caATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGcac < 1:1038682/150‑1 (MQ=255) cTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCAACACTTGGCGTGCTGCGTCt > 2:1149631/1‑150 (MQ=255) ggg‑‑gCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCg > 2:283027/1‑150 (MQ=255) g‑‑gCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTTATCAATGCATCAGTGCCACGTGTACGCCTGACACATTCAGCGCCTGCCAGTCCCCACCCTTGACCCCATCCAACAACCAGCAAGGAGCAGTTGGCGTGCGGCGTCTTGTCGTg < 2:375092/150‑1 (MQ=255) g‑‑gCA‑ACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTg > 2:227388/1‑150 (MQ=255) g‑‑gCA‑ACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCATGTAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTg > 1:1085456/1‑150 (MQ=255) gCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTgg > 1:1004021/1‑150 (MQ=255) gCA‑ACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCCTCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTgg > 1:71516/1‑150 (MQ=255) gCA‑ACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTgg > 2:873479/1‑150 (MQ=255) cA‑ACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGc < 2:71516/150‑1 (MQ=255) a‑aCGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACCGCCACCACTTGGCGTGCTGCGTCTTGTCGTGGCa > 2:181278/1‑150 (MQ=255) a‑aCGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCa > 1:1199911/1‑150 (MQ=255) aCGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACt > 2:1247153/1‑78 (MQ=255) aCGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACt < 1:1247153/78‑1 (MQ=255) cATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGa < 2:1004021/150‑1 (MQ=255) | ACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGG‑‑GCA‑ACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGA > NC_002947/4311213‑4311505 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |