New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 176362	NA (NA)	37 (0.640)	22/266	2.2	NA	intergenic (+6/‑457)	PP_5SA/PP_16SB	5S ribosomal RNA/16S ribosomal RNA
	?	NC_002947	= 176400	NA (NA)					intergenic (+44/‑419)	PP_5SA/PP_16SB	5S ribosomal RNA/16S ribosomal RNA

TTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCACCTCCCGCGACGA  >  NC_002947/176257‑176511                                                                                                          |                                                                                                                                                       ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTTATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                           <  1:320216/150‑1 (MQ=14) ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTaa                                                                                                           <  2:1163066/150‑1 (MQ=31)  tGGAACCACCTGAGCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAGTTCTCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAg                                                                                                          <  2:549506/150‑1 (MQ=12)    gAACCACCTGATACCATCGCGAACTCAGTAGTGAAACGGTGCATCGCCGATGGTAGTGTGGGGATTCCCCATGTGAGAGTAGGTCATCGTCAAGCTTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTa                                                                                                        >  2:858372/1‑150 (MQ=11)       ccaccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGaa                                                                                                     >  1:575495/1‑150 (MQ=34)            tGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTAACAGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCa                                                                                                <  1:1215778/150‑1 (MQ=255)              atcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATg                                                                                              >  1:575313/1‑150 (MQ=40)                                      aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGt                                                                      <  2:111955/150‑1 (MQ=255)                                              aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGt                                                              <  2:575313/150‑1 (MQ=255)                                                      tGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTaa                                                      >  1:479588/1‑150 (MQ=255)                                                            tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAggt                                                >  1:594776/1‑150 (MQ=255)                                                              tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAggtgg                                              >  2:140371/1‑150 (MQ=255)                                                                        ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAgtgt                                    >  1:1095548/1‑150 (MQ=255)                                                                        ccATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCT‑CCTGCCGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCAGGGTAAGCCGCAAGGTCAGGCAGTTGCAGTAAAAAGGTGGTTGACAgtgt                                    >  1:1150537/1‑150 (MQ=25)                                                                            gtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGttttt                                >  1:723598/1‑150 (MQ=255)                                                                               agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGaa                             <  1:90356/150‑1 (MQ=255)                                                                                  gTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTc                                                                                          >  1:584664/1‑86 (MQ=255)                                                                                  gTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTc                                                                                          <  2:584664/86‑1 (MQ=255)                                                                                   tAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCt                         >  2:1162932/1‑150 (MQ=255)                                                                                    aGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTg                        <  2:479588/150‑1 (MQ=255)                                                                                     ggTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGt                       >  2:59031/1‑150 (MQ=255)                                                                                      gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGta                      >  1:1304742/1‑150 (MQ=255)                                                                                         aTCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGt                                                                                                         >  2:726164/1‑64 (MQ=37)                                                                                         aTCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGt                                                                                                         <  1:726164/64‑1 (MQ=37)                                                                                          tCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGtattat                  >  1:581955/1‑150 (MQ=255)                                                                                              cAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCAc              <  1:59031/150‑1 (MQ=255)                                                                                               aaGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTattatt                 >  2:1164951/1‑146 (MQ=255)                                                                                               aaGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTattatt                 <  1:1164951/146‑1 (MQ=255)                                                                                                  aTTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAggtg                                               >  2:1094863/1‑113 (MQ=255)                                                                                                  aTTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAggtg                                               <  1:1094863/113‑1 (MQ=255)                                                                                                      aTTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTNAGTAGAAACTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTGTTTAAATGCTGTATTATTCACCTCCcgcg      >  1:34497/1‑150 (MQ=255)                                                                                                        ttCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTCGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCACCTCCCGcgac    <  2:575495/150‑1 (MQ=255)                                                                                                        atCGCAAAACTCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCACCTCCCGcgac    >  2:1208235/2‑150 (MQ=255)                                                                                                          cGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCACCTCCCGcgacga  <  1:1162932/150‑1 (MQ=255)                                                                                                          cGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCACCTCCCGcgacga  <  1:256093/150‑1 (MQ=255)                                                                                                          |                                                                                                                                                       TTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTG‑CGCATGCAGGTAGGGGTTTTGTCTTTAAGTAGAAGCTCATGCTTCGTCCTGGTAAGCCGCAAGGTCAGCCAGTTGCAGTAAAAAGGTGGTTGACAGTGTTTTTGAATGCTGTATTATTCACCTCCCGCGACGA  >  NC_002947/176257‑176511

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.