New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 823993	35 (0.540)	3 (0.050)	3/248	NT	8.0%	coding (505/762 nt)	ybhP	endo/exonuclease/phosphatase family protein
	?	NC_000913	= 824001	36 (0.580)					coding (497/762 nt)	ybhP	endo/exonuclease/phosphatase family protein

GTAGAGGAAATTGCACCGGAAACGTGCGCGCCGGGCGTCCGTGGGCGCGGGTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGCTTTTCCGGTCATCGGCGGCACAATGCGGCAGTAGAGCACGCCGCGCTTT  >  NC_000913/823855‑824144                                                                                                                                                   |                                                                                                                                                gtagAGGAAATTGCACCGGAAACGTGCGCGCCGGGCGTCCGTGGGCGCGGGTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTcc                                                                                                                                               >  2:420481‑M2/1‑147 (MQ=255) gtagAGGAAATTGCACCGGAAACGTGCGCGCCGGGCGTCCGTGGGCGCGGGTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTcc                                                                                                                                               <  1:440422‑M2/149‑3 (MQ=255) gtagAGGAAATTGCACCGGAAACGTGCGCGCCGGGCGTCCGTGGGCGCGGGTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTcc                                                                                                                                               <  2:277150‑M2/149‑3 (MQ=255)    gagGAAATTGCACCGGAAACGTGCGCGCCGGGCGTCCGTGGGCGCGGGTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTccggt                                                                                                                                            >  2:440422‑M2/1‑144 (MQ=255)     agGAAATTGCACCGGAAACGTGCGCGCCGGGCGTCCGTGGGCGCGGGTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTccggta                                                                                                                                           <  1:39332‑M2/149‑7 (MQ=255)             gCACCGGAAACGTGCGCGCCGGGCGTCCGTGGGCGCGGGTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTc                                                                                                                                   <  1:277150/149‑1 (MQ=255)                 cGGAAACGTGCGCGCCGGGCGTCCGTGGGCGCGGGTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCAccc                                                                                                                               <  2:39332/149‑1 (MQ=255)                        gTGCGCGCCGGGCGTCCGTGGGCGCGGGTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGc                                                                                                                        <  1:572722/149‑1 (MQ=255)                        gTGCGCGCCGGGCGTCCGTGGGCGCGGGTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGc                                                                                                                        <  3:292968/149‑1 (MQ=255)                                 gggCGTCCGTGGGCGCGGGTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGc                                                                                                               <  2:196858/149‑1 (MQ=255)                                                   gTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTg                                                                                             <  3:382975/149‑1 (MQ=255)                                                          tctcATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCa                                                                                      >  1:283189/1‑149 (MQ=255)                                                            tcATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACg                                                                                    >  1:244042/1‑149 (MQ=255)                                                               tCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAg                                                                                 >  4:299421/1‑149 (MQ=255)                                                                ccAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAgg                                                                                <  2:75511/149‑1 (MQ=255)                                                                 cAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGc                                                                               >  3:415003/1‑149 (MQ=255)                                                                    tCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTgcgc                                                                                          <  3:234815/135‑1 (MQ=255)                                                                    tCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTgcgc                                                                                          >  4:234815/1‑135 (MQ=255)                                                                    tCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAgg                                                                                <  3:572832/145‑1 (MQ=255)                                                                    tCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAgg                                                                                >  4:572832/1‑145 (MQ=255)                                                                    tCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCa                                                                            >  2:367807/1‑149 (MQ=255)                                                                      cGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGcc                                                                                                     <  2:336686/122‑1 (MQ=255)                                                                      cGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGcc                                                                                                     >  1:336686/1‑122 (MQ=255)                                                                            gCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACg                                                                    >  4:35071/1‑149 (MQ=255)                                                                             cACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGc                                                                   <  1:90173/149‑1 (MQ=255)                                                                                    aaCGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGcacatcac                                                            <  3:401477/149‑1 (MQ=255)                                                                                             ttAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGc                                                   >  3:337520/1‑149 (MQ=255)                                                                                              tAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCt                                                                                        >  3:240833/1‑111 (MQ=255)                                                                                              tAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCt                                                                                        <  4:240833/111‑1 (MQ=255)                                                                                                  tttttGCCGCCAGTCATTGAAATCACCCGCCACCAATTCCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGCTTTTc                                              <  3:6469/149‑1 (MQ=255)                                                                                                  tttttGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGCTTTTc                                              >  2:66112/1‑149 (MQ=255)                                                                                                  tttttGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGCTTTTc                                              >  2:180283/1‑149 (MQ=255)                                                                                                                     aaaTCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGCTTTTCCGGTCATCGGCGGCACAAt                           >  2:199357/1‑149 (MQ=255)                                                                                                                          aCCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCAtt                                                                                                                            >  4:39685/1‑47 (MQ=255)                                                                                                                          aCCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCAtt                                                                                                                            <  3:39685/47‑1 (MQ=255)                                                                                                                            ccGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGCTTTTCCGGTCATCGGCGGCACAATGCGGCAg                    >  2:454848/1‑149 (MQ=255)                                                                                                                               ccaccaATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGCTTTTCCGGTCATCGGCGGCACAATGCGGCAGTag                 <  2:63733/149‑1 (MQ=255)                                                                                                                                caccaATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGCTTTTCCGGTCATCGGCGGCACAATGCGGCAGTaga                <  3:60060/149‑1 (MQ=255)                                                                                                                                caccaATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGCTTTTCCGGTCATCGGCGGCACAATGCGGCAGTaga                <  1:454848/149‑1 (MQ=255)                                                                                                                                  ccaATACCGGTTCGCCGTccggtattggtggcgggtgatttcaatgactggcggcaaaaagctaatcatccgttaaaagtgcaggccggactggatgagatttttacccgcgcccacggacgcccggcgcgcacgtttccggtgcaatt              <  2:311495‑M2/149‑132 (MQ=255)                                                                                                                                   caATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGCTTTTCCGGTCATCGGCGGCACAATGCGGCAGTAGAGCa             >  1:187211/1‑149 (MQ=255)                                                                                                                                    aaTACCGGTTCGCCGTccggtattggtggcgggtgatttcaatgactggcggcaaaaagctaat                                                                                                 <  4:236559‑M2/64‑49 (MQ=255)                                                                                                                                    aaTACCGGTTCGCCGTccggtattggtggcgggtgatttcaatgactggcggcaaaaagctaat                                                                                                 >  3:236559‑M2/1‑16 (MQ=255)                                                                                                                                              cGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGCTTTTCCGGTCATCGGCGGCACAATGCGGCAGTAGAGCACGCCGCGCttt  >  2:110403/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                GTAGAGGAAATTGCACCGGAAACGTGCGCGCCGGGCGTCCGTGGGCGCGGGTAAAAATCTCATCCAGTCCGGCCTGCACTTTTAACGGATGATTAGCTTTTTGCCGCCAGTCATTGAAATCACCCGCCACCAATACCGGTTCGCCGTCCGGTAGCTCATTCACCCATTCGGCGAGCATCGCAAGCTGCGCCTGACGGTGCGCCTCACGCAGGCCCAGATGTACGCACATCACATGAATCGCTTTTCCGGTCATCGGCGGCACAATGCGGCAGTAGAGCACGCCGCGCTTT  >  NC_000913/823855‑824144

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