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

AAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCA  >  NC_000913/3214823‑3215114                                                                                                                                                |                                                                                                                                                     ggttcgacgcagatttgggtgctgccaaatcccagcggtttaagtcgcgtttcactggagaaactggttgaagcgtatcgcgagctggaccaggcgctggtagtgcgtgggcgataaaaaacgccaccttttcggtggcgtttTTTATc                                                                                                                                                 <  3:394024‑M2/6‑1 (MQ=255)   ggCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGc                                                                                                                                               <  3:148948/149‑1 (MQ=255)              aaCTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAc                                                                                                                                    >  2:106345/1‑149 (MQ=255)              aaCTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAc                                                                                                                                    <  1:34062/149‑1 (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                                                                                                   <  1:315497/149‑1 (MQ=255)                                               gTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAAc                                                                                                   >  2:576738/1‑149 (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:232438/149‑1 (MQ=255)                                                                 aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAAc                                                                                 <  3:217027/149‑1 (MQ=255)                                                                 aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAAc                                                                                 <  3:561854/149‑1 (MQ=255)                                                                 aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAAc                                                                                 >  4:430974/1‑149 (MQ=255)                                                                     ggTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGa                                                                             <  2:376435/149‑1 (MQ=255)                                                                     ggTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACAAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGa                                                                             <  3:430986/149‑1 (MQ=255)                                                                      gTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTg                                                                                                                           >  4:389617/1‑102 (MQ=255)                                                                      gTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTg                                                                                                                           <  3:389617/102‑1 (MQ=255)                                                                           ccGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTaaa                                                                       >  4:316070/1‑149 (MQ=255)                                                                                  aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTgg                                                                <  2:46177/149‑1 (MQ=255)                                                                                  aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTgg                                                                >  1:533154/1‑149 (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)                                                                                                                       aaacgccaccgaaaaggtggcgtttTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCtt                                                                                                       >  3:290057‑M2/26‑73 (MQ=255)                                                                                                                       aaacgccaccgaaaaggtggcgtttTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCtt                                                                                                       <  4:290057‑M2/48‑1 (MQ=255)                                                                                                                       aaacgccaccgaaaaggtggcgtttTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGt                           >  1:182690‑M2/26‑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)                                                                                                                             CACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGt                     <  4:32140‑M1/149‑1 (MQ=255)                                                                                                                                 ttttCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCgtttgt                 <  1:366691/149‑1 (MQ=255)                                                                                                                                           cgtttTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACt       <  4:394024‑M2/144‑1 (MQ=255)                                                                                                                                                tttATTGCCCACGGCCCACCCACGCCTGGGCCAGCTCGCGGTTCGGTTCCACCCGTTTTTCCAATGGAACCCCACTTTAACCGCTGGGGATTTGCAACACCCAAATCTTCGGTCAACCAATTGTGAGCGTTTGTTTCCCCCACTGTGCa  <  2:533154/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     AAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTGTTTCCCCCACTGTGCA  >  NC_000913/3214823‑3215114

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.