Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 1,317,152:1 | +GT | 100% | intergenic (+23/+5) | phrA → / ← yjpA | secreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 1,317,152 | 1 | . | G | 94.7% | 129.8 / 1.0 | 38 | intergenic (+23/+5) | phrA/yjpA | secreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base G (20/16); total (21/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
* | BS168_reference_genome | 1,317,152 | 2 | . | T | 94.7% | 106.2 / ‑0.2 | 38 | intergenic (+23/+5) | phrA/yjpA | secreted inhibitor of the activity of phosphatase RapA (quorum sensing)/putative enzyme |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base T (19/16); total (21/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GGTTTGTTGCTCGTTGCGGTCGGGTTCAGCTTTACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGG‑‑C‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACATACGCT > BS168_reference_genome/1317016‑1317274 || ggTTTGTTGCTCGTTGCGGTCGGGTTCAGCTTTACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAggggt < 6:265981/139‑3 (MQ=255) gttgCTCGTTGCGGTCGGGTTCAGCTTTACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGG‑‑G‑tcttt < 5:215194/138‑1 (MQ=255) cGTTGCGGGCGGGTTCAGCTTTACTCAGGTGATGGGTCCTTCAGGTGAAACCGCAAACACAGAAGGGAAAACATTTAATATTGCGGCCCGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGG‑‑G‑TCttttttattt > 8:187731/1‑135 (MQ=255) gTTCAGCTTTACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCa < 3:132922/139‑1 (MQ=255) tACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTAt > 7:217612/1‑139 (MQ=255) aGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGt < 8:14815/137‑1 (MQ=255) ggTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGACTTTTTTATTTCTTAAGCTTCGATTCTTTTATCGTAAg > 7:286902/1‑139 (MQ=255) gTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGc > 1:196504/1‑139 (MQ=255) gTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGc > 3:142324/1‑139 (MQ=255) tGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCt < 2:196504/139‑1 (MQ=255) gATGGTTCATGCAGGAGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTc < 5:10100/139‑1 (MQ=255) aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc > 5:191511/1‑139 (MQ=255) aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc > 7:41682/1‑139 (MQ=255) aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc > 7:287362/1‑139 (MQ=255) aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc > 6:139742/1‑139 (MQ=255) aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc > 5:186488/1‑139 (MQ=255) aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc > 3:300277/1‑139 (MQ=255) aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATc > 2:238473/1‑139 (MQ=255) aTGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTAAGCTCAGAAGATc > 1:189283/1‑139 (MQ=255) gTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTg > 3:91987/1‑139 (MQ=255) cAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGcc > 7:234020/1‑138 (MQ=255) gAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCa < 5:139742/139‑1 (MQ=255) aTTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTACATTCTTTTATCGTACGCTCAGAAGATCCACTTGCAACCAGCCGATCCGCATGGGCGATTTc > 2:213323/1‑139 (MQ=255) ttCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCg > 3:223014/1‑139 (MQ=255) tCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCg < 4:223014/138‑1 (MQ=255) cATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTc > 8:263825/1‑73 (MQ=255) cATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTc < 7:263825/73‑1 (MQ=255) tGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCg < 4:191086/132‑1 (MQ=255) tGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCg > 3:191086/1‑132 (MQ=255) cGGCACGAAATCAAACATGATGAATAAAAATAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGTATCCGCATGGCCGATTTCCGCTGCctct < 1:277492/138‑1 (MQ=255) cGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGt > 7:229739/1‑75 (MQ=255) cGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGt < 8:229739/75‑1 (MQ=255) aTCAAACATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAAt < 2:270571/139‑1 (MQ=255) cATGATGCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCtt < 8:287362/139‑1 (MQ=255) gatgCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGa > 6:26379/1‑139 (MQ=255) atgCATAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGGCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGTCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGAt < 1:213323/139‑1 (MQ=255) aTAAAAAAAGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACa < 6:186488/139‑1 (MQ=255) aaaaGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACATACGc < 8:41682/139‑1 (MQ=255) acaGACCCTTAGGGGTCTTTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACATACGCt < 7:187731/137‑1 (MQ=255) || GGTTTGTTGCTCGTTGCGGTCGGGTTCAGCTTTACTCAGGTGATGGTTCATGCAGGTGAAACAGCAAACACAGAAGGGAAAACATTTCATATTGCGGCACGCAATCAAACATGATGCATAAAAAAAGACCCTTAGGG‑‑C‑TTTTTATTTCTTCAGCTTCCATTCTTTTATCGTCAGCTCAGAAGATCCACTTGCCACCAGCGGATCCGCATGGCCGATTTCCGCTGCCTCTTCCAGTGAATCTGCTTCGATGACATACGCT > BS168_reference_genome/1317016‑1317274 |
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: ‑ |