| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 337,830 | Δ1 bp | 17.5% | intergenic (‑77/+52) | PP_0279 ← / ← PP_0280 | hypothetical protein/amino acid ABC transporter permease |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 337,830 | 0 | T | . | 17.5% | 106.9 / 20.1 | 40 | intergenic (‑77/+52) | PP_0279/PP_0280 | hypothetical protein/amino acid ABC transporter permease |
| Reads supporting (aligned to +/- strand): ref base T (16/17); new base . (3/4); total (19/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 = 7.67e-01 | |||||||||||
TCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAA > NC_002947/337690‑337968 | tCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTac < 1:632194/150‑1 (MQ=255)tCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGAATGGCAAGGATGTGTGGTGGTTTTGAGTGCCTCTTAGAAAACATGCCCGCTCCAACAGGGATGTCGCA‑GCAATCAAGGTTGTCGATac > 1:186512/1‑150 (MQ=255) aCTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCtgt < 1:1016658/150‑1 (MQ=255) aCTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCtgt < 1:344239/150‑1 (MQ=255) cggggcGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGAATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCtgtg < 2:587064/145‑1 (MQ=21) cTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACATCtgtg < 2:902531/150‑1 (MQ=255) tGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTgg < 2:1190791/150‑1 (MQ=21) tCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAg < 2:406804/150‑1 (MQ=255) gAGGGCCTGGAAGCGGGTGCTGAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACTCCCGCTCCCACAGGTATAGCGCG‑ACGTTCAAGGTTGGCGCTACACCTGTGGGAGCg < 2:443690/150‑1 (MQ=255) ggCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGc < 2:288360/150‑1 (MQ=255) tggtTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCtt < 2:33429/150‑1 (MQ=255) aaGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTg > 1:129229/1‑76 (MQ=255) aaGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTg < 2:129229/76‑1 (MQ=255) gcggtttttGAGGGTCTCTTACGGGGAACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTCGGCGCTAAACTTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAgc < 2:186512/145‑1 (MQ=37) ttttGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCggg > 2:855035/1‑83 (MQ=255) ttttGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCggg < 1:855035/83‑1 (MQ=255) cccGCTCCCACAGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGaa > 2:948732/1‑136 (MQ=18) cccGCTCCCACAGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGaa < 1:948732/136‑1 (MQ=18) ccGCTCCCACAGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCaa < 2:803051/114‑1 (MQ=14) ccGCTCCCACAGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCaa > 1:803051/1‑114 (MQ=14) tCCCACAGGTGTAGCGCCAACCTTGAATG‑TGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc > 2:3844/1‑150 (MQ=21) ccACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga > 1:34736/1‑150 (MQ=255) ccACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGACAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga > 2:747597/1‑150 (MQ=255) cacaGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCAcgc < 1:321215/103‑1 (MQ=255) cacaGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCAcgc > 2:321215/1‑103 (MQ=255) cacaGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgag > 1:429291/1‑150 (MQ=255) cacaGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga > 2:422380/1‑149 (MQ=255) cacaAGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACga < 1:422380/149‑1 (MQ=255) acaGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTg < 1:343428/71‑1 (MQ=255) acaGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTg > 2:343428/1‑71 (MQ=255) acaGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACgaga > 2:1126414/1‑150 (MQ=255) gTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCAc > 1:733392/1‑150 (MQ=255) gTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCAc > 1:977391/1‑150 (MQ=255) tataGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGc > 1:177081/1‑129 (MQ=255) tataGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGc < 2:177081/129‑1 (MQ=255) c‑acaTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGaa > 2:374617/1‑150 (MQ=255) acaTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATc > 1:86398/1‑83 (MQ=255) acaTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATc < 2:86398/83‑1 (MQ=255) aTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAg < 1:1116365/141‑1 (MQ=255) aTTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAg > 2:1116365/1‑141 (MQ=255) | TCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCC‑ACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAA > NC_002947/337690‑337968 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |