Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_002947	1,385,534	Δ69 bp	89.6%	coding (861‑929/1266 nt)	oprD ←	basic amino acid specific porin OprD

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 1385533	2 (0.060)	21 (0.700)	21/234	NT	89.6%	coding (930/1266 nt)	oprD	basic amino acid specific porin OprD
	?	NC_002947	1385603 =	3 (0.100)					coding (860/1266 nt)	oprD	basic amino acid specific porin OprD

GCGTACATAACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/1385407‑1385533 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggccaggGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGTTGAAGTCGAACACCAGGCCTTGATTGTCGGA  >  NC_002947/1385603‑1385751                                                                                                                                                                                                                                                                                      GCGTACATAACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTT                                                                                                                                >  2:236950/1‑150    TACATAACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTA                                                                                                                             <  1:850559/150‑1     ACATAACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAG                                                                                                                            <  1:260113/150‑1       ATAAGGCGTCATGAGTGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGC                                                                                                                          >  2:425755/1‑150          ACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTCCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCTTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCC                                                                                                                       <  2:615220/150‑1                  TGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCAAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCACCCGACAGGC                                                                                                               >  1:367656/1‑150                       GTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAAC                                                                                                          <  2:541673/150‑1                           GGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGCATTTGGTCAGTGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGC                                                                                                      <  1:320520/150‑1                                GGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTG                                                                                                 <  1:425755/150‑1                                  AATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTC                                                                                               <  2:367656/150‑1                                   ATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCC                                                                                              <  2:422248/150‑1                                     GCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAG                                                                                            <  1:236950/150‑1                                               GCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTG                                                                                                                               <  2:326405/105‑1                                               GCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTG                                                                                                                               >  1:326405/1‑105                                               GCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAA                                                                                                                                                  <  1:520243/86‑1                                               GCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAA                                                                                                                                                  >  2:520243/1‑86                                                           TCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGG                                                                      >  2:679053/1‑150                                                             GTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTA                                                                    <  1:820499/150‑1                                                                  GAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGG                                                               >  1:800559/1‑150                                                                                          CGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCAT                                       <  1:398624/150‑1                                                                                              GTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGAT                                   <  1:679053/150‑1                                                                                               TTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGGGTCCAGAGGTGATCGCTGATGTTGTAGGCGCCCGACCGCCAGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATG                                  <  1:191851/150‑1                                                                                                     TCGGAGCGGGCCACGGAGTTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGTTGAAG                            <  1:466068/150‑1                                                                                                                       TTGGCCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGTTGAAGTCGAACACCAGGCCTTGA          <  2:553427/150‑1                                                                                                                           CCAGGGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCAAGGTTGTTCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGTTGAAGTCGAACACCAGGCCTTGATTGT      >  2:292877/1‑150                                                                                                                               GGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGTTGAAGTCGAACACCAGGCCTTGATTGTCGGA  >  1:331195/1‑150                                                                                                                               GGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGTTGAAGTCGAACACCAGGCCTTGATTGTCG    <  1:567548/148‑1                                                                                                                               GGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGTTGAAGTCGAACACCAGGCCTTGATTGTCG    >  2:567548/1‑148                                                                                                                               GGTGAAGGTGTGAGCGCCGATGTTCAAGACGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGTTGAAGTCGAACACCAGGCCTTGATTGTCGGA  >  1:331165/1‑150                                                                                                                                                                                                                                                                                      GCGTACATAACGGGTCATGAAGGTCAGGCCTGGAATGCCGTATTCGGCGAAGTTCAGGTCGTAGCGAGCCTGCCAGGATTTCTCGTCTTCGGCGTTGAAGTCGGAGCGGGCCACGGAGTTGGCCAGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/1385407‑1385533 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggccaggGTGAAGGTGTGAGCGCCGATGTTGTAGGCGCCCGACAGGCTGAACGCGCGGTTGTCCAGCTTGTCGCCATCGAAGGCGCGGTACTCGGCAGAGCCTTCGCTCTTGGTGTCATAGATGTTGAAGTCGAACACCAGGCCTTGATTGTCGGA  >  NC_002947/1385603‑1385751

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ 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.