Predicted mutation
evidence seq id position mutation annotation gene description
JC AE015451 4,980,585 +GGC intergenic (‑24/+44) PP_4387 ← / ← flgE conserved exported protein of unknown function/Flagellar hook protein FlgE

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AE015451 = 49805850 (0.000)8 (0.360)
+GGC
8/310 1.3 100% intergenic (‑24/+44) PP_4387/flgE conserved exported protein of unknown function/Flagellar hook protein FlgE
?AE015451 4980586 = 0 (0.000)intergenic (‑25/+43) PP_4387/flgE conserved exported protein of unknown function/Flagellar hook protein FlgE

GTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGG  >  AE015451/4980394‑4980750
                                                                                                                                                                                               |                                                                                                                                                                     
gtggtaggGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaac                                                                                                                                            <  1:284463‑M1/216‑28 (MQ=255)
              agTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattaccaacaccaaagc                                                                                                                                    >  1:336487‑M1/1‑178 (MQ=255)
                                                                                           aaCGGGGCGCGGGCAGCGGGGGGAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATGCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaat                                                                                                                           <  2:296831‑M1/145‑45 (MQ=255)
                                                                                                                                           aaCAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtattggcattggcctggtagttgctctgcgccttgatcaggtttaccagctcgccggtcaggtcgacgt  >  1:251179‑M1/1‑53 (MQ=255)
                                                                                                                                                                                              ggggccgctttgcggccccaattacaaacaccaaagccaatacaatcatca                                                                                                                      <  1:155725‑M1/51‑50 (MQ=255)
                                                                                                                                                                                              ggggccgctttgcggccccaattacaaacaccaaagccaatacaatcatca                                                                                                                      >  2:155724‑M1/1‑2 (MQ=255)
                                                                                                                                                                                               gggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgcct                                      <  1:224996‑M1/130‑130 (MQ=255)
                                                                                                                                                                                               gggccgctttgcggccccaattacaaacaccaaagccaatacaatcatcaggtcatctgaatgatggtctgcatgatggtgctttcagtggaaatggtcttggcattggcctggtagttgctctgcgcct                                      >  2:224995‑M1/1‑1 (MQ=255)
                                                                                                                                                                                               |                                                                                                                                                                     
GTTGTAGGGGCCAGAGTGGCGCACCCAGCCTTTGCCCGGGTTGTGGGCCGAACACAGGTAGCGGCCGGTGACCAGGCTTTGCCATTGATAGAACGGGGCGGGGGCAGCGGGGGCAGGCAGAACCAGCACCAGCAAGGGGAACAAGGCCAGTAACGACAGCCGCTTCATTCGCGTGAACCTTCACATTCTTGGCGCTTTGCGGCCCCAATTACAAACACCAAAGCCAATACAATCATCAGGTCATCTGAATGATGGTCTGCATGATGGTGCTTTCAGTGGAAATGGTCTTGGCATTGGCCTGGTAGTTGCTCTGCGCCTTGATCAGGTTTACCAGCTCGCCGGTCAGGTCGACGTTGG  >  AE015451/4980394‑4980750

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.