New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 3336220	26 (0.580)	3 (0.070)	3/232	NT	11.5%	coding (419/843 nt)	yghU	GSH‑dependent disulfide bond oxidoreductase
	?	NC_002947	= 3336246	22 (0.520)					coding (393/843 nt)	yghU	GSH‑dependent disulfide bond oxidoreductase

GCTTGGCCTCCATGGTGAAGCGGTTGATGGCGTACTCGAACTTTTCCGGTGCATACACGTAGAAGTGGCCAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCT  >  NC_002947/3336086‑3336358                                                                                                                                        |                                                                                                                                           gCTTGGCCTCCATGGTGAAGCGGTTGATGGCGTACTCGAACTTTTCCGGTGCATACACGTAGAAGTGGCCAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGAtt                                                                                                                             >  1:1148432/1‑150 (MQ=255)   ttGGCCTCCATGGTGAAGCGGTTGATGGCGTACTCGAACTTTTCCGGTGCATACACGTAGAAGTGGCCAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGAtttt                                                                                                                           <  2:1148432/150‑1 (MQ=255)       ccTCCATGGCGAACCGGTTAATGGCGCCCACGGAAGTTTACGAAGCAATCACGGAGAAGCGGCCAAAACGACTGCTCTGTGAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGAATTTGgca                                                                                                                       >  2:12113/1‑150 (MQ=255)         tCCATGGTGAAGCGGTTGATGGCGTACTCGAACTTTTCCGGTGCATACACGTAGAAGTGGCCAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGgcagc                                                                                                                     >  2:1111607/1‑150 (MQ=255)                        ttGATGGCGTACTCGAACTTTTCCGGTGCATACACGTAGAAGTGGCCAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCtttgctgccaaaatcgccggctgcgcgcaccgaatgcc                                                                                                      <  2:1033866‑M1/150‑39 (MQ=255)                        ttGATGGCGTACTCGAACTTTTCCGGTGCATACACGTAGAAGTGGCCAAAACCACCGCCCAGGTAACGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCCCCCCGCCGGCGATTTTGGCAGCAACGCACCGAACTtc                                                                                                      >  1:610227/1‑150 (MQ=255)                         tGATGGCGTACTCGAACTTTTCCGGTGCATACACGTAGAAGTGGCCAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTtct                                                                                                     <  1:910605/150‑1 (MQ=255)                                    tCGAACTTTTCCGGTGCATACACGTAGAAGTGGCCAAAACCAACGCCCAGGAAAGACGACACCCCCATCAGACA‑AACAAACAATCTCACTCAGACTTTGTACGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTAc                                                                                          >  2:739079/1‑150 (MQ=255)                                         cTTTTCCGGTGCATACACGTAGAAGTGGCCAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCtttgctgccaaaatcgccggctgcgcgcacc                                                                                                             >  1:1028111‑M1/1‑95 (MQ=255)                                         cTTTTCCGGTGCATACACGTAGAAGTGGCCAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCtttgctgccaaaatcgccggctgcgcgcacc                                                                                                             <  2:1028111‑M1/126‑32 (MQ=255)                                                    cATACACGTAGAAGTGGCCAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTg                                                                          <  1:636351/150‑1 (MQ=255)                                                     aTACACGTAGAAGTGGCCAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGa                                                                         >  2:1031340/1‑150 (MQ=255)                                                                      cAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTctct                                                                                                   >  2:921176/1‑107 (MQ=255)                                                                      cAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTctct                                                                                                   <  1:921176/107‑1 (MQ=255)                                                                         aaCCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACcggcg                                                     >  1:1008856/1‑150 (MQ=255)                                                                                       tAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCgg                                       >  2:796656/1‑150 (MQ=255)                                                                                           gcgcCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGa                                   >  2:630553/1‑150 (MQ=255)                                                                                                gccCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAgcagg                              <  2:1172828/150‑1 (MQ=255)                                                                                                          gCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCt                    >  1:205790/1‑150 (MQ=255)                                                                                                            cAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGACTTTGGCAGCAAAGCACCGAACTTCTGGGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCCACCGACCAGGACAGGGATCGTc                  >  2:153634/1‑150 (MQ=255)                                                                                                              gAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCATCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGa                >  2:92987/1‑150 (MQ=255)                                                                                                                  aacCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTg            >  2:914138/1‑150 (MQ=255)                                                                                                                   acCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTgg           <  1:1031340/150‑1 (MQ=255)                                                                                                                     cAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTgggg         >  2:1029124/1‑150 (MQ=255)                                                                                                                     cAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGAAGCAACGGCTGGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTgggg         <  2:31020/150‑1 (MQ=255)                                                                                                                        t‑tCAGGCATTCGGTGCGCGCAGCCGGCCATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTg      >  1:202761/1‑150 (MQ=255)                                                                                                                          tCAGGCATTCGGTGCGCGCAGCCGGCGATCTTGGCAGCAAAGCACCGAACTCCCCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGa     <  1:12113/150‑1 (MQ=255)                                                                                                                            aGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCg                                                                     >  2:557463/1‑84 (MQ=255)                                                                                                                            aGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCg                                                                     <  1:557463/84‑1 (MQ=255)                                                                                                                             ggCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAgcaccgaatgcctgaactggttgttctggcagatgggggcggcgccttacctgggcggtggttttggccacttctacgtgtatgcaccggaaaagttcgagtacgccatcaa  <  2:1033866‑M2/150‑113 (MQ=255)                                                                                                                                    ggTGCGCGCAGCCGGCGATTTTGGCAGCAAAgcaccgaatgcctgaactggttgttctggcagatgggggcggcgccttacctgggcggtggttttggccacttctacgtgtatgcaccggaaaag                   >  1:1028111‑M2/1‑31 (MQ=255)                                                                                                                                    ggTGCGCGCAGCCGGCGATTTTGGCAGCAAAgcaccgaatgcctgaactggttgttctggcagatgggggcggcgccttacctgggcggtggttttggccacttctacgtgtatgcaccggaaaag                   <  2:1028111‑M2/126‑96 (MQ=255)                                                                                                                                        |                                                                                                                                           GCTTGGCCTCCATGGTGAAGCGGTTGATGGCGTACTCGAACTTTTCCGGTGCATACACGTAGAAGTGGCCAAAACCACCGCCCAGGTAAGGCGCCGCCCCCATCTGCCAGAACAACCAGT‑TCAGGCATTCGGTGCGCGCAGCCGGCGATTTTGGCAGCAAAGCACCGAACTTCTCTGCCAGGTACAGCAGGATCGAGCCTGACTCGAACACCCGCACCGGCGGCTCCACGCTGCGGTCGAGCAGGGCAGGGATCTTCGAGTTGGGGTTGACCT  >  NC_002947/3336086‑3336358

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