New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4074001 =	NA (NA)	4 (0.100)	4/244	NT	18.8%	intergenic (‑495/+129)	mdtA/PP_3586	multidrug transporter membrane fusion protein/transposase
	?	NC_002947	4302589 =	18 (0.420)					intergenic (+1035/‑265)	PP_3774/PP_5627	hypothetical protein/hypothetical protein

TGTTCACTCTTCTATACGTGGTCACAGCAAGCTGTGCCGACAGGGATTAACTCCGATGCCATTTGGACCACTCTCTCATTCGGGAGCAAAGCAGCACACCAGTGCGGTCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCTGGAAGGGCTGGGTAATCTCTCATGTGGCATTGGGATTTGG  >  NC_002947/4302441‑4302737                                                                                                                                                     |                                                                                                                                                     tcaactggtggccgcttgagtgcttagaaaacaatgggttgcattttgtttgatgagagatcacccagcccttccgggctgcgctgtcgcacagcaaacaggatgcgcatgacgtcagatctgcactgacatcgagatctagaaaagcAt                                                                                                                                                     >  1:767284‑M2/149‑150 (MQ=255)   aactggtggcctcttgagtgcttagaaaacagtgggttgcatgttgtttgatgagagatcacccagcccttccgggctgcgctgtcgcacagcaaacaggatgcgcatgacgtcagatctgcactgacatcgagatctagaaaagcATTc                                                                                                                                                   >  2:69732‑M2/147‑150 (MQ=255)   aactggtggccgcttgagtgcttagaaaacaatgggttgcattttgtttgatgagagatcacccagcccttccgggctgcgctgtcgcacagcaaacaggatgcgcatgacgtcagatctgcactgacatcgagatctagaaaagcatta                                                                                                                                                   >  1:905774‑M2/147‑149 (MQ=255)   aactggtggccgcttgagtgcttagaaaacaatgggttgcattttgtttgatgagagatcacccagcccttccgggctgcgctgtcgcacagcaaacaggatgcgcatgacgtcagatctgcactgacatcgagatctagaaaagcATTc                                                                                                                                                   >  2:751363‑M2/147‑150 (MQ=255)   aactggtggccgcttgagtgcttagaaaacaatgggttgcattttgtttgatgagagatcacccagcccttccgggctgcgctgtcgcacagcaaacaggatgcgcatgacgtcagatctgcactgacatcgagatctagaaaagcATTc                                                                                                                                                   >  1:794723‑M2/147‑150 (MQ=255)   aactggtggccgcttgagtgcttagaaaacaatgggttgcattttgaatgattagagatcacccagcccttccgggctgcgctgtcgcacagcaaacaggatgcgcatgacgtcagatctgcactgacatcgagatctagaaaagcATTc                                                                                                                                                   >  1:202773‑M2/147‑150 (MQ=255)                  gTGGTCACAGCAAGCTGTGCCGACAGGGATTAACTCCGATGCCATTTGTACCACTCTCTCATTCGGGAGCAAAGCAGCACACCAGTGCGGTCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGGCCTAcc                                                                                                                                    >  2:444755/1‑150 (MQ=25)                            cAAGCTGTGCCGACAGGGATTAACTCCGATGCCATTTGGACCACTCTCTCATTCGGGAGCAAAGCAGCACACCAGTGCGGTCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTc                                                                                                                          <  1:284353/150‑1 (MQ=35)                                 tgtgCCGACAGGGATTAACTCCGATGCCATTTGGACCACTCTCTCATTCGGGAGCAAAGCAGCACACCAGTGCGGTCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGa                                                                                                                     >  1:460560/1‑150 (MQ=35)                                  gtgCCGACAGGGATTAACTCCGATGCCATTTGGACCACTCTCTCATTCGGGAGCAAAGCAGCACACCAGTGCGGTCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGaa                                                                                                                    >  2:1054231/1‑150 (MQ=35)                                         gcattttgtttgatgagagatcacccagcccttccgggctgcgctgtcgcacagcaaacaggatgcgcatgacgtcagatctgcactgacatcgagatctagaaaagcATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCtttt                                                                                                             <  1:1005899‑M2/42‑1 (MQ=255)                                            ggATTAACTCCGATGCCATTTGGACCACTCTCTCATTCGGGAGCGAAGCAGCACACCAGTGCGGTCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTTGGCGGCTGACCTCCCCGGATCAGTCGGCGAATGCTTTTCTa                                                                                                          >  1:904243/1‑150 (MQ=255)                                              aTTAACTCCGATGCCATTTGGACCACTCTCTCATTCGGGAGCAAAGCAGCACACCAGTGCGGTCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGa                                                                                                        >  1:206832/1‑150 (MQ=35)                                                    tCCGATGCCATTTGGACCACTCTCTCATTCGGGAGCAAAGCAGCACACCAGTGCGGTCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGa                                                                                                  >  2:112672/1‑150 (MQ=33)                                                     atgagagatcacccagcccttccgggctgcgctgtcgcacagcaaacaggatgcgcatgacgtcagatctgcactgacatcgagatctagcaaagtATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGAt                                                                                                 <  2:405300‑M2/54‑1 (MQ=255)                                                            cATTTGGACCACTCTCTCATTCGGGAGCAAAGCAGCACACCAGTGCGGTCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTg                                                                                          <  2:287619/150‑1 (MQ=32)                                                             aTTTGGACCACTCTCTCATTCGGGAGCAAAGCAGCACACCAGTGCGGTCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGc                                                                                         >  1:1113963/1‑150 (MQ=32)                                                              tttGGACCACTATCTCAATCGGTAGCAAAGCAGCACACCAGTGCGGTCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCa                                                                                        <  2:527864/150‑1 (MQ=12)                                                                                           acagcaaacaggatgcgcatgacgtcagatctgcactgacatcgagatctagaaaagcATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGa                                                                                      <  1:500435‑M2/65‑1 (MQ=255)                                                                                           acagcaaacaggatgcgcatgacgtcagatctgcactgacatcgagatctagaaaagcATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGa                                                                                      >  2:500435‑M2/59‑123 (MQ=255)                                                                                                            tCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGGCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCt                                          >  2:635978/1‑150 (MQ=11)                                                                                                               aCGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCTGGa                                       <  2:481383/150‑1 (MQ=32)                                                                                                                cGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCTGGaa                                      <  2:206832/150‑1 (MQ=32)                                                                                                                       tttCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCTGGAAGGGCTgg                               <  1:658759/150‑1 (MQ=32)                                                                                                                           ccGCGCATTGGTAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCTGGAAGGGCTGGGTaa                           <  1:444755/150‑1 (MQ=21)                                                                                                                                  ttGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCTGGAAGGGCTGGGTAATCTCTCa                    <  1:455137/150‑1 (MQ=35)                                                                                                                                      gAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGACAATGACCATCTCGATGTGAGATCTTTTCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCTGGAAGGGCTGGGTAATCTCTCAtgtg                >  1:411913/1‑150 (MQ=11)                                                                                                                                            tAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATATAAACAGATACGCATCCTGTTTGCTGTGCGACAGCGCAGCCTGGAAGGGCTGGGTAATCTCTCATGTGGCAttg          >  1:1040075/1‑150 (MQ=14)                                                                                                                                                    tATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCTGGAAGGGCTGGGTAATCTCTCATGTGGCAttgggatttgg  <  2:460560/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TGTTCACTCTTCTATACGTGGTCACAGCAAGCTGTGCCGACAGGGATTAACTCCGATGCCATTTGGACCACTCTCTCATTCGGGAGCAAAGCAGCACACCAGTGCGGTCTACGGCCTCTTTCCCGCGCATTGGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCTGGAAGGGCTGGGTAATCTCTCATGTGGCATTGGGATTTGG  >  NC_002947/4302441‑4302737

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