Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,489,830 | Δ1 bp | 17.1% | coding (264/597 nt) | PP_4825 → | MarC family membrane protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,489,830 | 0 | G | . | 17.1% | 92.5 / 17.3 | 35 | coding (264/597 nt) | PP_4825 | MarC family membrane protein |
Reads supporting (aligned to +/- strand): ref base G (11/18); new base . (4/2); total (15/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01 |
ACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATG > NC_002947/5489683‑5489969 | aCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTTCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGc < 1:210270/150‑1 (MQ=255) ccTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAgg < 1:139832/150‑1 (MQ=255) cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCAT‑GCCCAggg < 1:43512/150‑1 (MQ=255) gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCaa < 2:143511/150‑1 (MQ=255) gcgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCcggcg < 1:374953/150‑1 (MQ=255) gcgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCcggcg < 1:487892/150‑1 (MQ=255) cgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCcggcgg < 2:382853/150‑1 (MQ=255) gTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATAGATCCGACATCGGTGCTGGCGGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGt > 1:261531/1‑150 (MQ=255) cGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACa < 2:225194/150‑1 (MQ=255) ccTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTACTTCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGT‑CAGGGCCCAGGGCAAGCCGGCGGTACAGGCAg > 1:348826/1‑150 (MQ=255) gCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCCTCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCCGGGCAAGCCGGCGGTACAGGCAGACAACGtt < 2:467476/150‑1 (MQ=255) tGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGACGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCCGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGa < 1:539494/150‑1 (MQ=255) aCCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAt < 2:530842/143‑1 (MQ=255) aCCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAt > 1:530842/1‑143 (MQ=255) cTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACg > 2:580531/1‑130 (MQ=255) cTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACg < 1:580531/130‑1 (MQ=255) gATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCAt > 2:278951/1‑150 (MQ=255) gCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCAT‑GCCCAGGGcc < 1:438556/75‑2 (MQ=38) gCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCAT‑GCCCAGGGcc > 2:438556/1‑74 (MQ=38) gCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATc < 1:395150/119‑1 (MQ=255) gCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATc > 2:395150/1‑119 (MQ=255) ccTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTg < 2:261531/150‑1 (MQ=255) cGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctg > 1:312318/1‑150 (MQ=255) gggTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGC‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGt < 2:348826/149‑1 (MQ=255) gTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTc < 1:284023/75‑1 (MQ=255) gTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTc > 2:284023/1‑75 (MQ=255) tGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCt < 2:347804/100‑1 (MQ=255) tGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCt > 1:347804/1‑100 (MQ=255) ttCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAtgt < 2:578794/68‑1 (MQ=255) ttCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAtgt > 1:578794/1‑68 (MQ=255) ttgCCCTGGGCCAT‑GCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTggg > 2:573526/3‑150 (MQ=255) gCCCTGGGCCAT‑GCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGAt > 1:548357/1‑150 (MQ=255) cccTGGGCCAT‑GCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATg > 2:369824/1‑150 (MQ=255) tGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGc < 2:317163/75‑1 (MQ=255) tGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGc > 1:317163/1‑75 (MQ=255) | ACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATG > NC_002947/5489683‑5489969 |
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: ‑ |