New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3214941 =	NA (NA)	4 (0.070)	3/250	NT	8.9%	intergenic (+53/+26)	rpoD/mug	RNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase
	?	NC_000913	3214966 =	41 (0.670)					intergenic (+78/+1)	rpoD/mug	RNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase

CCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGT  >  NC_000913/3214795‑3215089                                                                                                                                                   |                                                                                                                                                     ccGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACaaa                                                                                                                                                    >  4:229254/1‑149 (MQ=255)  cgcgAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAAc                                                                                                                                                   <  1:144719/149‑1 (MQ=255)   agtgggggaaacaaacgctcaccattggttcgacgcagatttgggtgctgccaaatcccagcggtttaagtcgcgtttcactggagaaactggttgaagcgtatcgcgagctggaccaggcgctggtagtgcgtgggcgataaaAAACg                                                                                                                                                  <  4:394024‑M1/5‑1 (MQ=255)     gAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGcc                                                                                                                                                >  1:376435/1‑149 (MQ=255)         gTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCt                                                                                                                                            >  2:571124/1‑149 (MQ=255)                tCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGAACAGGCGTTACGCATCACCCGGCACTAGGCCCTCTGCACAAACTCCACCTTTTCgat                                                                                                                                     <  4:498910/149‑3 (MQ=255)                 cAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTg                                                                                                                                    <  3:191828/149‑1 (MQ=255)                    aTCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCg                                                                                                                                 >  4:217027/1‑149 (MQ=255)                    aTCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCg                                                                                                                                 >  4:561854/1‑149 (MQ=255)                       accattggttcgacgcagatttgggtgctgccaaatcccagcggtttaagtcgcgtttcactggagaaactggttgaagcgtatcgcgagctggaccaggcgctggtagtgcgtgggcgataaaAAACGCCACCTTTTCGGTGGCGttt                                                                                                                              >  1:182690‑M1/125‑149 (MQ=255)                               ggCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGc                                                                                                                      <  3:148948/149‑1 (MQ=255)                                          aaCTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAc                                                                                                           <  1:34062/149‑1 (MQ=255)                                          aaCTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAc                                                                                                           >  2:106345/1‑149 (MQ=255)                                                gTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAgcgc                                                                                                     <  1:106345/149‑1 (MQ=255)                                                      cGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTc                                                                                               <  2:5016/149‑1 (MQ=255)                                                            gTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTc                                                                                         >  4:79785/1‑149 (MQ=255)                                                              tCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTcgc                                                                                       <  4:556136/149‑1 (MQ=255)                                                              tCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTcgc                                                                                       <  2:405972/149‑1 (MQ=255)                                                                  aaGTGCCGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATa                                                                                   >  2:366691/1‑149 (MQ=255)                                                                           gTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAAc                                                                          >  2:576738/1‑149 (MQ=255)                                                                           gTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAAc                                                                          <  1:315497/149‑1 (MQ=255)                                                                            tAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAAcc                                                                         <  3:60512/149‑1 (MQ=255)                                                                                ttCCTGGACTATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGtt                                                                     >  4:459753/1‑149 (MQ=255)                                                                                ttCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGtt                                                                     >  1:586030/1‑149 (MQ=255)                                                                                 tCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGttt                                                                    <  3:133755/149‑1 (MQ=255)                                                                                   cTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTtct                                                                  >  4:361695/1‑149 (MQ=255)                                                                                      gACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCa                                                               >  3:397319/1‑149 (MQ=255)                                                                                          aTTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGa                                                           >  1:508922/1‑149 (MQ=255)                                                                                           ttAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGaa                                                          <  2:138101/149‑1 (MQ=255)                                                                                             aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAAc                                                        >  4:430974/1‑149 (MQ=255)                                                                                             aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAAc                                                        <  4:232438/149‑1 (MQ=255)                                                                                             aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAAc                                                        <  3:217027/149‑1 (MQ=255)                                                                                             aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAAc                                                        <  3:561854/149‑1 (MQ=255)                                                                                                 ggTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGa                                                    <  2:376435/149‑1 (MQ=255)                                                                                                 ggTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACAAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGa                                                    <  3:430986/149‑1 (MQ=255)                                                                                                  gTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTg                                                                                                  <  3:389617/102‑1 (MQ=255)                                                                                                  gTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTg                                                                                                  >  4:389617/1‑102 (MQ=255)                                                                                                   aagcgtatcgcgagctggaccaggcgctggtagtgcgtgggcgataaaAAACGCCACCTTTTCGGTGGCGttt                                                                                                                              >  3:290057‑M1/49‑73 (MQ=255)                                                                                                   aagcgtatcgcgagctggaccaggcgctggtagtgcgtgggcgataaaAAACGCCACCTTTTCGGTGGCGttt                                                                                                                              <  4:290057‑M1/25‑1 (MQ=255)                                                                                                       ccGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTaaa                                              >  4:316070/1‑149 (MQ=255)                                                                                                              aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTgg                                       >  1:533154/1‑149 (MQ=255)                                                                                                              aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTgg                                       <  2:46177/149‑1 (MQ=255)                                                                                                              aGGCGTTACGCCGCACCCGGCACTACGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTgg                                       <  3:370673/149‑1 (MQ=255)                                                                                                                   ttACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGAttt                                  <  4:397319/149‑1 (MQ=255)                                                                                                                   ttACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGAttt                                  <  1:144260/149‑1 (MQ=255)                                                                                                                     aCGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTgg                                <  4:371075/149‑1 (MQ=255)                                                                                                                       gccgcACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGgca                              <  4:105873/149‑1 (MQ=255)                                                                                                                            aCCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCAcc                         <  2:99398/149‑1 (MQ=255)                                                                                                                                ggCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACACTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCaaa                     >  1:416894/1‑149 (MQ=255)                                                                                                                                     tAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGc                <  4:392046/149‑1 (MQ=255)                                                                                                                                       ggCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGt              <  2:555439/149‑1 (MQ=255)                                                                                                                                            tctGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAAc         >  4:296565/1‑149 (MQ=255)                                                                                                                                             gataaaAAACGCCACCGAAAAGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAAcc        <  3:394024‑M1/143‑1 (MQ=255)                                                                                                                                                 acaAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATg    >  2:471916/1‑149 (MQ=255)                                                                                                                                                   AAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGgc                               >  2:385504‑M1/1‑120 (MQ=255)                                                                                                                                                   AAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGgc                               <  1:385504‑M1/120‑1 (MQ=255)                                                                                                                                                   AAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGt  <  3:195356‑M1/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CCGCGAACGTATCCGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGT  >  NC_000913/3214795‑3215089

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