New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	1235874 =	49 (0.860)	3 (0.060)	3/250	NT	6.1%	coding (315/852 nt)	glpF	aquaglyceroporin
	?	NC_002947	1235912 =	46 (0.840)					coding (277/852 nt)	glpF	aquaglyceroporin

GGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAGATGATACTG  >  NC_002947/1235765‑1236048                                                                                                                                                    |                                                                                                                                         ggTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCACACCTGGGCCAGCATg                                                                                                                                         <  7:339859/149‑1 (MQ=255) atggccatctacctcaccgccggtatttccggtgcgcacctgaacccagccgtgagtatcgccctcacactgttcgcagggttcgacaagcgcaagctgcccttctacaTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATg                                                                                                                                         >  1:109484‑M1/110‑149 (MQ=255)                ccGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTg                                                                                                                          >  7:37897/1‑149 (MQ=255)                 cGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGt                                                                                                                                        >  8:263819/1‑134 (MQ=255)                 cGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGt                                                                                                                                        <  7:263819/134‑1 (MQ=255)                       cAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGt                                                                                                                   >  7:91042/1‑149 (MQ=255)                             cAGGCTGCCTTCGCTACCTCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAAc                                                                                                              <  3:314840/148‑1 (MQ=255)                             cAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAAcc                                                                                                             <  2:324590/149‑1 (MQ=255)                                gCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTg                                                                                                          >  3:80036/1‑149 (MQ=255)                                      ttCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGaa                                                                                                      >  2:191534/1‑147 (MQ=255)                                      ttCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGaa                                                                                                      >  1:48212/1‑147 (MQ=255)                                         gCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGGTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAgtg                                                                                                 <  8:259193/149‑1 (MQ=255)                                             cccagccgtgagtatcgccctcacactgttcgcagggttcgacaagcgcaagctgcccttctacaTGCTGCGCCACAGAATGCGCCGAAAACCTGGGCCAGCATGTAg                                                                                                                                      <  4:313269‑M1/43‑1 (MQ=255)                                             cccagccgtgagtatcgccctcacactcttcgcagggttcgacaagcgcaagctgcccttctacaTGCTGCGCCACAGAATGCGCCGCATACCTCGGCCAGCATGTAg                                                                                                                                      >  3:313269‑M1/66‑108 (MQ=255)                                              cgcgCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAgg                                                                                            >  4:349931/1‑149 (MQ=255)                                               gcgcAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAggg                                                                                           >  5:106800/1‑149 (MQ=255)                                               gcgcAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAggg                                                                                           <  4:80036/149‑1 (MQ=255)                                                         gCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCAc                                                                                 <  6:154502/149‑1 (MQ=255)                                                           gttggCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACgg                                                                               <  8:335980/146‑1 (MQ=255)                                                           gTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACgg                                                                               >  3:311599/1‑149 (MQ=255)                                                               gCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTgg                                                                           <  8:91042/149‑1 (MQ=255)                                                                 tttttCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGAGCCACAGAATGCGCCGCATACCTTGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGt                                                                         >  6:367562/4‑149 (MQ=255)                                                                 ttGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGt                                                                         >  2:203848/1‑149 (MQ=255)                                                                    ttCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTc                                                                       <  4:311599/148‑1 (MQ=255)                                                                       gAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGt                                                                   >  5:219726/1‑149 (MQ=255)                                                                         aaTCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTg                                                                            <  1:69789/138‑1 (MQ=255)                                                                         aaTCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTg                                                                            >  2:69789/1‑138 (MQ=255)                                                                                 aaCAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCgg                                                          >  7:223570/1‑148 (MQ=255)                                                                                   cAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGaaa                                                       >  4:28242/1‑149 (MQ=255)                                                                                    aGATTGCTGTACAGCGTGTAGACCAGTGCTCCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGaaa                                                       >  1:30739/1‑148 (MQ=255)                                                                                     gATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGAAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACTGTGTGAGGGCGATACTCACGGCTCGGTTCAGGTGCGCACCGGAAATa                                                     >  8:285989/1‑149 (MQ=255)                                                                                      aTTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATAc                                                    >  6:288216/1‑149 (MQ=255)                                                                                           tGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACcggcg                                               >  4:17256/1‑149 (MQ=255)                                                                                           tGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACcggcg                                               >  1:32971/1‑149 (MQ=255)                                                                                              aCAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGt                                             >  7:99482/1‑148 (MQ=255)                                                                                               cAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGa                                           <  6:106800/149‑1 (MQ=255)                                                                                               cAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGa                                           <  3:349931/149‑1 (MQ=255)                                                                                               cAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGa                                           <  3:28242/149‑1 (MQ=255)                                                                                                  cGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGt                                        <  3:17256/149‑1 (MQ=255)                                                                                                         ctacatgctggccgaggtatgcggcgcattctgtggcgcagcaTGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAAGAGTGTGAGGGCGATACTCACGGCTggg                                                                          >  3:313269‑M2/44‑108 (MQ=255)                                                                                                         ctacatgctggcccaggttttcggcgcattctgtggcgcagcaTGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTggg                                                                          <  4:313269‑M2/65‑1 (MQ=255)                                                                                                         aCCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGc                                 <  5:288216/149‑1 (MQ=255)                                                                                                            catgctggcccaggtatgcggcgcattctgtggcgcagcaTGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCAt                              >  1:109484‑M2/41‑149 (MQ=255)                                                                                                             gTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATg                             <  2:32971/149‑1 (MQ=255)                                                                                                                    gcCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacgc                      >  3:270134/1‑149 (MQ=255)                                                                                                                      cacaGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGCGGGCGATACTCACGGCGGGGTTCAGGTGCGCACCGGAAATACcggcgg                                              >  2:277572/1‑123 (MQ=255)                                                                                                                      cacaGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacgcc                     >  2:278156/1‑148 (MQ=255)                                                                                                                        caGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacgcccac                  <  2:30739/149‑1 (MQ=255)                                                                                                                           aaTGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacgcccacgcc               >  8:234444/1‑149 (MQ=255)                                                                                                                           aaTGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGacgcccacgcc               <  6:219726/149‑1 (MQ=255)                                                                                                                                gccgcATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAga          <  2:48212/149‑1 (MQ=255)                                                                                                                                    cATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAAAAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACAGGAAAGACCGGCTGTGAGGTAGATTGCACTGATGCCGCCCACGCCCCAgatga       >  6:98391/1‑148 (MQ=255)                                                                                                                                        ccTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACATTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAGATGATACt   <  5:367562/148‑1 (MQ=255)                                                                                                                                        ccTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATAGTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAGATGATACTg  <  7:234444/149‑1 (MQ=255)                                                                                                                                        ccTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAGATGATACt   <  1:203848/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                         GGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTACAGCGTGTAGACCAGTGCTGCGCCACAGAATGCGCCGCATACCTGGGCCAGCATGTAGAAGGGCAGCTTGCGCTTGTCGAACCCTGCGAACAGTGTGAGGGCGATACTCACGGCTGGGTTCAGGTGCGCACCGGAAATACCGGCGGTGAGGTAGATGGCCATGCTGACGCCCACGCCCCAGATGATACTG  >  NC_002947/1235765‑1236048

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