New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3995709 =	49 (0.870)	3 (0.060)	3/250	NT	6.4%	coding (127/708 nt)	yigA	conserved protein YigA
	?	NC_000913	3995732 =	42 (0.790)					coding (150/708 nt)	yigA	conserved protein YigA

GGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATC  >  NC_000913/3995569‑3995844                                                                                                                                             |                                                                                                                                         ttgctctgggctgcaggtagaggaggcgataaaacaggccttcgttggcgatagcctgttccatcaacggcgccatgttctcttccagaacatgaatatgattacgtgcgcgggccatgtgccactcgaccaacgaaacgGTg                                                                                                                                        >  2:614171‑M1/141‑143 (MQ=255)   agactgcgctgcaggtagagtaggcgataaaacaggccttcgttggcgatagcctgttccatcaacagcgccatgttctcttccagaacatgcacatgagtacctgcgagggctacggggaacgcgcacaccgaaacgGTGc‑c                                                                                                                                     <  2:1065929‑M1/5‑1 (MQ=255)    tttATTCATCTATGAAGCAACCAGGGGAATAACTGCAGGAAACACTCACGGAGCTCGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑Cg                                                                                                                                    <  1:200880/143‑1 (MQ=255)     ttATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGc                                                                                                                                   >  1:344243/1‑143 (MQ=255)     gcggctctgcaggcagagtgggcgataaaaccgggcttcgttggcgatagcctgttccatcaacagcgccatgttctcttccagaacatgaatatgattacgtgcgcgggccatgtgccactcgaccaacgaaacgGTGC‑CGc                                                                                                                                   >  2:977978‑M1/137‑143 (MQ=255)          cATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCg                                                                                                                              <  1:5635/143‑1 (MQ=255)            tCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGt                                                                                                                            >  1:776186/1‑143 (MQ=255)               aTGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑GGCAGCCGGGAcg                                                                                                                         >  2:800577/1‑143 (MQ=255)                 gAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTAcgcg                                                                                                                       >  2:152070/1‑143 (MQ=255)                      aaCCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCAcc                                                                                                                  >  1:217805/1‑143 (MQ=255)                       aCCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCCGTAGACGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCg                                                                                                                 >  2:401178/1‑143 (MQ=255)                             ggAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAACATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGt                                                                                                           >  1:328556/1‑143 (MQ=255)                               aagaaCTGCAGGAACCACTCACGGACCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTg                                                                                                         <  1:951085/143‑1 (MQ=255)                               aaacaggccttcgttggcgatagcctgttccatcaacagcgccatgttctcttccagaacatgaatatgattacgtgcgcgggccatgtgccactcgaccaacgatacgtGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTg                                                                                                         >  2:503223‑M1/111‑143 (MQ=255)                                 gaaCTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGt                                                                                                       >  1:945520/1‑143 (MQ=255)                                 caaCTGGAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGCGTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGt                                                                                                       <  1:527409/142‑1 (MQ=255)                                  gcaggccctcgttgggtacagcctgttcactcaaccgcgcgatgttctcttccagaacatgaatatgattacgtgcgcgggccatgtgccactcgaccaacgaaacgGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTc                                                                                                      >  1:201970‑M1/108‑143 (MQ=255)                                  aaCTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTc                                                                                                      >  1:53025/1‑143 (MQ=255)                                     tGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACGGTTTCGTTGGGCGAg                                                                                                   >  2:199601/1‑143 (MQ=255)                                     tGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTGGAg                                                                                                   >  1:780366/1‑143 (MQ=255)                                        aGGACACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTgg                                                                                                <  2:1036215/143‑1 (MQ=255)                                         ttcgttggcgatagcctgttccatcaacagcgccatgttctcttccagaacatgaatatgattacgtgcgcgggcgatgtgccactcgaccaacgaaacgGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGc                                                                                               >  1:719435‑M1/101‑143 (MQ=255)                                            atcacTCCCCGCGCTTGACGACCGGGCGGTTGTCGATCATCTGATTAAAAATCCTGAGTTTTTTATCCGTAGTGCGCGCGTAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGcaca                                                                                            <  1:1078761/141‑1 (MQ=255)                                               acTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTGGAGTGGCACGTgg                                                                                         >  1:90488/1‑143 (MQ=255)                                                 tCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTTGTTGGTCGAGTGGCACATGGcc                                                                                       >  2:1087690/1‑143 (MQ=255)                                                       aGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGAAGCAACGAGACGTGTGC‑AGAAACAGGCACGCGGCACCGTTTCGTTGAGCGCGTGGCACATGGCCCgcgcc                                                                                 >  2:233302/1‑142 (MQ=255)                                                          ttGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGCACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGACCCGCGCACgc                                                                              >  2:145905/1‑142 (MQ=255)                                                            gatgaCCGGGCTGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTaa                                                                            <  1:622107/143‑1 (MQ=255)                                                            gaggcccGGGGGGCTGCCGATTCTGTGATTAAAAATCCTGAGTTTTTTCTCCGTAATGCGCGCGCAGTAGAAGCGAGACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTaa                                                                            <  1:1048207/138‑1 (MQ=255)                                                                 ccGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCata                                                                       >  1:706938/1‑143 (MQ=255)                                                                         ttGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCCTATTCATGtt                                                               >  1:1094327/1‑143 (MQ=255)                                                                              gATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTgta                                                          >  2:3069/1‑141 (MQ=255)                                                                               aTTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGaa                                                         <  1:593756/143‑1 (MQ=255)                                                                                 tATCGGATGAGAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAGGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAaga                                                       <  2:112389/143‑1 (MQ=255)                                                                                    cTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATATGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCCCGTAATCATATTCCTGTTCTGGCAgagac                                                    >  2:898822/1‑142 (MQ=255)                                                                                      gATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCGTATTCATGTTCTGGAAGAGAACa                                                  >  1:892319/1‑143 (MQ=255)                                                                                          aaaaaTCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTTTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGgc                                              >  1:1051360/1‑143 (MQ=255)                                                                                           aaaaTCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGCAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCCTATTCATCTCCTGGAAGGGCACATGgcg                                             >  2:741391/1‑143 (MQ=255)                                                                                            aaaTCCTGAGTTTTTTATCCGTAAGGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTCCTTGAAGAGAACATGgcgc                                            <  2:969805/143‑1 (MQ=255)                                                                                            aaaTCCTGAGTGTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGGGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGgcgc                                            <  1:485358/143‑1 (MQ=255)                                                                                               tCCTGGGGTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGt                                         <  1:734309/143‑1 (MQ=255)                                                                                                  tGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCCCATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAAAGAACATGGGGCTGTTGa                                      >  2:380077/1‑143 (MQ=255)                                                                                                  tGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGa                                      <  1:472381/143‑1 (MQ=255)                                                                                                   gtcttgtttATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGCGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGAt                                     <  1:196473/137‑1 (MQ=255)                                                                                                   gggTTTTTTATCGGTACGGCGCGCGCAGTAGAGGCGAGACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGAt                                     <  2:584472/141‑1 (MQ=255)                                                                                                   gAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGAt                                     <  1:428283/143‑1 (MQ=255)                                                                                                    aGTTTTTTCGCGGGAATGCGAGGGCAGTGGAAGGGATACGGGTGC‑GGCAGCCGGTACGCGGCACCGTTTCGTTGGTCGGGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATg                                    <  1:227683/143‑1 (MQ=255)                                                                                                    aGTTTCTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATg                                    <  1:542413/143‑1 (MQ=255)                                                                                                      ttttttATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACACCAGCATATTCATGTTCTTGAAGAGAACATGGCGCCGTTGATGGa                                  >  2:272715/1‑143 (MQ=255)                                                                                                      gtttttATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTATCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAGGAGAACATGGCGCTGTTGATGGa                                  >  1:540140/2‑143 (MQ=255)                                                                                                           tATCCGTAGTGCGCGCGCAGTAGAGGCTATACGTGGGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAgg                             <  2:800177/143‑1 (MQ=255)                                                                                                           tATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAgg                             >  2:85137/1‑143 (MQ=255)                                                                                                               ggTAATGCGCGCGCAGTCGAGGCTGTACTTCAGCACGTTTTC‑GTACTCGTCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGGACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTAt                         <  2:1071021/142‑1 (MQ=255)                                                                                                                                 gAAGCCATACGCGGGC‑CGCATCCGGTACGCGGCCCCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGtt       <  2:911483/143‑1 (MQ=255)                                                                                                                                      gATACGTTTGC‑CGCATCCGGTACGTGGTACCGTTTCGTGTGTCGAGTGGCACATGGCCCGCGTATGTTATCATATTCATGTTCTGGAAGAGAACATGGCGCTGGTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATc  <  1:122178/143‑1 (MQ=255)                                                                                                                                             |                                                                                                                                         GGATTTATTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGC‑CGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATC  >  NC_000913/3995569‑3995844

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