| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 2,939,284 | A→G | 26.5% | intergenic (‑230/+210) | PP_2570 ← / ← PP_2572 | transposase/hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 2,939,284 | 0 | A | G | 26.5% | 38.6 / 18.6 | 34 | intergenic (‑230/+210) | PP_2570/PP_2572 | transposase/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (12/13); new base G (3/6); total (15/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.97e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01 | |||||||||||
CTCTGTCCCGCGCACACGAACCCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCCGCGAAAGGATCAACA > NC_002947/2939226‑2939432 | agttccgcagcgacgctgacattttgaaaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGc < 8:184646‑M1/117‑1 (MQ=255) acgctgacattttgaaaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTcc > 7:184646‑M1/21‑149 (MQ=255) atctattttgaaaaagAAAAGGCGGGTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACa > 7:155340‑M1/17‑149 (MQ=255) ttttgaaaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTcctacgggc < 4:288278‑M1/136‑5 (MQ=255) ttttgaaaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACa < 8:155340‑M1/133‑1 (MQ=255) ccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCGACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTAc < 2:162785/149‑1 (MQ=17) ccTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTAc > 1:102554/1‑149 (MQ=32) tAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActg > 4:209068/1‑149 (MQ=33) tAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActg > 3:125595/1‑149 (MQ=33) aaCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATTTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActgc < 1:175087/149‑1 (MQ=21) aCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActggt > 1:177754/1‑147 (MQ=18) aCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActgct > 5:236535/1‑149 (MQ=35) aCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGCACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActgct < 7:72265/149‑1 (MQ=25) aCAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTAAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActgct > 2:284803/1‑149 (MQ=25) aaaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTg > 5:357777‑M1/7‑125 (MQ=255) aaaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTg < 6:357777‑M1/119‑1 (MQ=255) aaaagAAAAGGCGGTTGGCCCGGTAGGGCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTActgctg < 6:113444‑M1/143‑1 (MQ=255) aaaGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGt > 8:304477/1‑149 (MQ=37) aaGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTa < 4:68346/149‑1 (MQ=38) aaaaGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGaa > 1:112921/1‑149 (MQ=40) aaaGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAAc < 7:325382/149‑1 (MQ=255) aGGCGGCTGGCCCTACCGGGCAAACCGCTTTTTTTTTTTAAAAATGTCAGCGTCGCTGCGGAACTGAAGCCTTACACCTCATTTTTGCTGCGCTACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACgg < 2:351773/149‑1 (MQ=14) ggcggcTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACggg > 8:133047/1‑149 (MQ=255) ggcTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTc < 7:65007/149‑1 (MQ=255) gcTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCa > 4:89695/1‑149 (MQ=255) cccTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGGCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCGGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCcgc < 5:121593/149‑1 (MQ=255) cTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCcgcg > 6:142090/1‑148 (MQ=255) ccGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCCGCGaaa < 7:304477/148‑1 (MQ=255) ccGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCCGCGaa < 5:142090/147‑1 (MQ=255) ccGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCCGCGAAAg > 5:162211/1‑149 (MQ=255) cGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCCGCGAAAgg < 6:21824/149‑1 (MQ=255) ggCCAGCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCCGCGAAAGGAt < 5:113444/149‑1 (MQ=255) gCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCCGCGAAAGGATc < 7:296106/149‑1 (MQ=255) gCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCCGCGAAAGGATc < 7:279994/149‑1 (MQ=255) aaCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTa < 6:359929/120‑1 (MQ=40) aaCCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTa > 5:359929/1‑120 (MQ=40) aCCGCCTTTTCTTTTTCAAACTGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCCGCGAAAGGATCAaca < 7:102964/149‑1 (MQ=255) | CTCTGTCCCGCGCACACGAACCCTAACAAAGCAAAAGGCGGCTGGCCCTACCGGGCCAACCGCCTTTTCTTTTTCAAAATGTCAGCGTCGCTGCGGAACTGAAGCCGTACACCTCATGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGTACTGCTGTAGGAACGGGTTCACCCGCGAAAGGATCAACA > NC_002947/2939226‑2939432 |
| 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: ‑ |