Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 132,463 | G→A | 10.4% | intergenic (+37/‑115) | cc → / → dsbA | cytochrome c4/thiol:disulfide interchange protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 132,463 | 0 | G | A | 10.4% | 126.5 / 4.7 | 48 | intergenic (+37/‑115) | cc/dsbA | cytochrome c4/thiol:disulfide interchange protein |
Reads supporting (aligned to +/- strand): ref base G (24/19); new base A (3/2); total (27/21) | |||||||||||
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.74e-01 |
AACCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGG > NC_002947/132322‑132608 | aaCCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGc > 5:19104/1‑149 (MQ=255) aCCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCt > 7:23212/1‑149 (MQ=255) aaCGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTCGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCtttt < 6:345486/149‑1 (MQ=255) aTGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtttcatca < 5:19718/149‑6 (MQ=255) ggCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgt > 5:203887/1‑149 (MQ=255) ggCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgt > 8:145716/1‑149 (MQ=255) gCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgtc < 4:352391/149‑1 (MQ=255) cATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGGTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGACCTTTTTTCCGTCCGCAGcc > 4:341044/1‑149 (MQ=255) cATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGGTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCGTTTTTTCCGTCCGCAGcc > 4:42645/1‑149 (MQ=255) ccATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTac < 6:6708/149‑1 (MQ=255) cATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTaca > 7:214508/1‑149 (MQ=255) aCCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTACCGTCCGCAGCCGTTACACTACAg > 7:319466/1‑148 (MQ=255) ccATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGaa < 6:284239/149‑1 (MQ=255) gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgtc > 2:338962/1‑119 (MQ=255) gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTtccgtc < 1:338962/119‑1 (MQ=255) gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGc < 8:264408/147‑1 (MQ=255) gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGc < 6:47684/147‑1 (MQ=255) gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtttc < 7:325336/112‑2 (MQ=255) gCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCAGCCCGGCCGCCCtttc > 8:325336/1‑111 (MQ=255) aGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCg < 6:19104/147‑1 (MQ=255) cAACCATGACATCGAAGCGTAGGCCAGCTACATCCAGGGCCTGCAATAACTTTCGGTTAAAGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTCCCCTTTTTTCCATACGCAGCCGTTAAACTACAGAACACCAGCCCTTCCAGGcat > 3:5076/1‑147 (MQ=255) aTGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGtct > 1:99721/1‑148 (MQ=255) tGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGa > 2:201197/1‑129 (MQ=255) tGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCCTCCGCAGCCGTTACACTACAGAACTCGa < 1:201197/129‑1 (MQ=255) aCATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGt > 1:94804/1‑149 (MQ=255) aTCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGc < 3:243940/149‑1 (MQ=255) tCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTg > 3:333093/1‑147 (MQ=255) gAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAAGTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCAGCAGCCGTTACACTACAGAACTCGCGCCATTCCCGCCCTTTCTCAGTGCAgg > 3:264066/1‑149 (MQ=255) gAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAgg > 6:164306/1‑149 (MQ=255) aGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTc < 6:203887/149‑1 (MQ=255) ggaaGCTTCATCCAGGGGCTGCATTAACTTTCGTTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGcgtcgt < 4:5076/146‑1 (MQ=255) cATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGAcc > 1:216409/1‑148 (MQ=255) tCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCAt < 2:94804/149‑1 (MQ=255) cATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAg > 7:107245/1‑148 (MQ=255) cATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGATACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAg > 3:142440/1‑148 (MQ=255) aaCTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTACCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGt < 8:319466/148‑1 (MQ=255) cTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTaaa < 8:214508/149‑1 (MQ=255) tttCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCGTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAg > 5:265503/1‑149 (MQ=255) cGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCa < 7:145716/148‑1 (MQ=255) aaTGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGt > 4:264569/1‑148 (MQ=255) aTGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTaa < 6:265503/149‑1 (MQ=255) aTGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTa > 1:296341/1‑148 (MQ=255) tAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGAt > 3:274620/1‑149 (MQ=255) tCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATtctc < 2:99721/149‑1 (MQ=255) aaaaaaGGGCAGCCCGGCCGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgct > 4:143020/4‑149 (MQ=255) aaaaaGGGCAGCCCGGCCGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgctg > 4:44626/3‑149 (MQ=255) gaAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAgcgctgc > 7:40126/1‑149 (MQ=255) gggCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGcgctgcgctg < 8:40126/149‑1 (MQ=255) ggcggcCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTgg < 2:296341/149‑1 (MQ=255) ggcCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAAc < 3:240061/125‑1 (MQ=255) ggcCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAAc > 4:240061/1‑125 (MQ=255) | AACCGCACCAACGATGGCGATGCCATGACCATGCGCACCATCGCTGGCAAGCTGAGCAACCATGACATCGAAGCGTTGGCCAGCTACATCCAGGGCCTGCATTAACTTTCGGTTAATGTTGTAGGTCAATGATGAAAGGGCGGCCGGGCTGCCCTTTTTTCCGTCCGCAGCCGTTACACTACAGAACTCGAGCCCTTCCCGGCCTGTCTCAGTGCAGGTCGCGTCGTGGCGACCCATGACTGCTCAGGAGTAAAGCATGCGTAAACTGATTCTCAGCGCTGCGCTGG > NC_002947/132322‑132608 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |