Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 2,557,839 | T→G | intergenic (+108/‑43) | BW25113_RS12770 → / → BW25113_RS12775 | hypothetical protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,557,839 | 0 | T | G | 82.3% | 85.0 / 14.2 | 34 | intergenic (+108/‑43) | BW25113_RS12770/BW25113_RS12775 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (3/3); new base G (10/18); total (13/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.53e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01 |
GGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGA > NZ_CP009273/2557705‑2557973 | ggtagAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGtt < 2:284494/138‑1 (MQ=255) ggtagAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGtt < 2:451756/138‑1 (MQ=255) ggtagAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGtt < 1:603666/138‑1 (MQ=255) tagAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAAc < 2:71587/139‑1 (MQ=255) ttGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAAt > 1:587682/1‑139 (MQ=255) gAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATg > 1:243177/1‑138 (MQ=255) attgCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATgg < 1:651300/139‑1 (MQ=255) aTAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATaa < 2:133668/139‑1 (MQ=255) aacaacCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTa > 1:357118/1‑139 (MQ=255) ccGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAAttt > 2:651490/1‑139 (MQ=255) ccGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAAttt < 2:587682/139‑1 (MQ=255) ggCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCaa > 2:79719/1‑86 (MQ=255) ggCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCaa < 1:79719/86‑1 (MQ=255) ggCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGc > 1:203055/1‑139 (MQ=255) ggCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGc > 2:413655/1‑139 (MQ=255) ggCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGc > 2:339951/1‑139 (MQ=255) cccTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGAc < 2:393232/139‑1 (MQ=255) cTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACAt > 2:721442/1‑139 (MQ=255) gCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCt > 2:534685/1‑139 (MQ=255) tgATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGc < 1:326150/138‑1 (MQ=255) aTCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGa < 1:413655/138‑1 (MQ=255) tCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATg > 2:158997/1‑139 (MQ=255) cGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATgg < 2:357118/136‑1 (MQ=255) tCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa < 1:721442/139‑1 (MQ=255) tCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa < 1:339951/139‑1 (MQ=255) cAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGATATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa < 1:328205/138‑1 (MQ=255) atatTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGAc > 1:750520/1‑109 (MQ=255) atatTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGAc < 2:750520/109‑1 (MQ=255) cTTCCGTTACGCTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAgg < 2:203055/139‑1 (MQ=255) ttCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGGTTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa > 2:568485/1‑128 (MQ=255) ttCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa < 1:568485/128‑1 (MQ=255) cTTACCCTCCACCAATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGt < 2:438176/139‑1 (MQ=255) caccaATACGCGGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCaaa < 1:315894/139‑1 (MQ=255) cgcgGGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga < 1:241268/139‑1 (MQ=255) | GGTAGAGTATTTGAAAAACGTAATTGATTGCGAATAGCAACAACTAAACAACCAATACACCCGGCCCCTCGCCGGGTTTTTTGTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGA > NZ_CP009273/2557705‑2557973 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |