New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3670807 =	69 (0.670)	3 (0.030)	3/252	NT	3.8%	coding (1216/1650 nt)	treF	cytoplasmic trehalase
	?	NC_000913	3670847 =	86 (0.870)					coding (1256/1650 nt)	treF	cytoplasmic trehalase
Rejected: Frequency below/above cutoff threshold.

AACTGGAGAGCGCCATCGCCAACATCTCGGCGCTGAAAGGCGAGAAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACG  >  NC_000913/3670661‑3670954                                                                                                                                                   |                                                                                                                                                    tgcccagccgttgggtttatcccactgtccaccggtttcgtactcgcttgccagaatcccgccaggtgtcagtaaccgactgcgcacggcgtttgccagacgatcggcctgttcatggttcgccataccgacatagagtggcacaaTgg                                                                                                                                                   >  2:606029‑M1/147‑149 (MQ=255) aaCTGGAGAGCGCCATCGCCAACATCTCGGCGCTGAAAGGCGAGAAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTgg                                                                                                                                                   <  1:654273/149‑1 (MQ=255)  gcccagccgttgggtttatcccactgttcaccggtttcgtactcgcttgccagaatcccgccaggtgtcagtaaccgactgcgcacggcgtttgccagacgatcggcctgttcatggttcgccataccgacatagagtggcacaaTGGc                                                                                                                                                  >  2:195696‑M1/146‑149 (MQ=255)  gcccagccgttgggtttatcccactgttcaccggtctcgtactcgcttgccagaatcccgccaggtgtcagtaaccgactgcgcacggcgtttgccagacgatcggcctgttcatggttcgccataccgacatagagtggcacaaTGGc                                                                                                                                                  >  3:836851‑M1/146‑149 (MQ=255)      gagaGCGCCATCGCCAACATCTCGGCGCTGAAAGGCGAGAAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGAcg                                                                                                                                              >  4:286705/1‑149 (MQ=255)                    cAACATCTCGGCGCTGAAAGGCGAGAAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGgcgc                                                                                                                                <  3:771265/149‑1 (MQ=255)                      aCATCTCGGCGCTGAAAGGCGAGAAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTg                                                                                                                              <  4:188018/149‑1 (MQ=255)                            cGGCGCTGAAAGGCGAGAAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTcc                                                                                                                        <  4:889085/149‑1 (MQ=255)                                   gAAGGGCGAGAAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCGgccg                                                                                                                 >  2:153185/1‑149 (MQ=255)                                   gAAAGGCGAGAAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTgccg                                                                                                                 >  2:255932/1‑149 (MQ=255)                                   gAAAGGCGAGAAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTgccg                                                                                                                 >  1:776344/1‑149 (MQ=255)                                      aGGCGAGAAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCa                                                                                                              >  4:681249/1‑149 (MQ=255)                                      aGGCGAGAAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCa                                                                                                              <  3:286705/149‑1 (MQ=255)                                          gagaAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATtgt                                                                                                          >  2:781837/1‑149 (MQ=255)                                          actcgcttgccagaatcccgccaggtgtcagtaaccgactgcgcacggcgtttgccagacgatcggcctgttcatggttcgccataccgacatagagtggcaccaTGGCGACGCGAACAACTGGCGCTGTTTTCCTCGGCCGCCATtgt                                                                                                          >  4:665261‑M1/106‑149 (MQ=255)                                           agaAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATtgtg                                                                                                         <  4:445951/149‑1 (MQ=255)                                            gaAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGc                                                                                                        >  1:641296/1‑149 (MQ=255)                                            gaAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGc                                                                                                        <  2:415549/149‑1 (MQ=255)                                               agagACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCAc                                                                                                     <  3:917144/149‑1 (MQ=255)                                                  gaCAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCActct                                                                                                  <  1:601231/149‑1 (MQ=255)                                                   aCAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTa                                                                                                 >  2:510009/1‑149 (MQ=255)                                                    cagaatcccgccaggtgtcagtaaccgactgcgcacggcgtttgccagacgatcggcctgttcatggttcgccataccgacatagagtggcacaat                                                                                                                                                     <  4:364212‑M1/1‑1 (MQ=255)                                                    cagaatcccgccaggtgtcagtaaccgactgcgcacggcgtttgccagacgatcggcctgttcatggttcgccataccgacatagagtggcacaat                                                                                                                                                     >  3:364212‑M1/96‑96 (MQ=255)                                                      gAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGt                                                                                              >  2:780337/1‑149 (MQ=255)                                                      gAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGt                                                                                              >  2:891888/1‑149 (MQ=255)                                                        aGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCg                                                                                            <  1:781837/149‑1 (MQ=255)                                                          cACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGt                                                                                          <  1:638667/149‑1 (MQ=255)                                                           aCTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTa                                                                                         <  1:891888/149‑1 (MQ=255)                                                            cTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGTGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTAt                                                                                        >  4:688933/1‑149 (MQ=255)                                                              gTTGCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATgg                                                                                      <  4:315166/149‑1 (MQ=255)                                                                   gCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAAc                                                                                 >  2:907367/1‑149 (MQ=255)                                                                    ccAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAAcc                                                                                >  2:182159/1‑149 (MQ=255)                                                                      aGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCAt                                                                              >  1:12487/1‑149 (MQ=255)                                                                      aGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCAt                                                                              >  4:358748/1‑149 (MQ=255)                                                                        aaaGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCActct                                                                                                  >  2:308198/1‑127 (MQ=255)                                                                        aaaGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCActct                                                                                                  <  1:308198/127‑1 (MQ=255)                                                                         aaGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGaa                                                                           <  1:510009/149‑1 (MQ=255)                                                                          aGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAAc                                                                          >  2:512661/1‑149 (MQ=255)                                                                           gCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACa                                                                         >  1:330902/1‑149 (MQ=255)                                                                           gCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACa                                                                         >  2:445858/1‑149 (MQ=255)                                                                            ccAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAg                                                                        >  4:454291/1‑149 (MQ=255)                                                                             cAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAgg                                                                       >  1:206271/1‑149 (MQ=255)                                                                              aGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGc                                                                      <  4:32494/149‑1 (MQ=255)                                                                              aGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGc                                                                      >  2:252057/1‑149 (MQ=255)                                                                              aGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGc                                                                      <  4:274538/149‑1 (MQ=255)                                                                                tGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCg                                                                    >  3:555561/1‑149 (MQ=255)                                                                                   ccGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATc                                                                 <  1:892216/149‑1 (MQ=255)                                                                                   ccGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATc                                                                 <  3:733421/149‑1 (MQ=255)                                                                                      tCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACTGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTc                                                              <  4:496702/149‑1 (MQ=255)                                                                                       cgcgATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCt                                                             <  3:454291/149‑1 (MQ=255)                                                                                        gcgATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGGCTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTg                                                            <  1:153185/149‑1 (MQ=255)                                                                                        gcgATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTg                                                            <  3:91350/149‑1 (MQ=255)                                                                                          gATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGc                                                          >  4:445219/1‑149 (MQ=255)                                                                                              cGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAAc                                                      >  4:384902/1‑149 (MQ=255)                                                                                               ggTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACg                                                     >  2:787919/1‑149 (MQ=255)                                                                                                gTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGc                                                    >  1:488107/1‑149 (MQ=255)                                                                                                   aaCCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTgg                                                           <  4:606612/139‑1 (MQ=255)                                                                                                   aaCCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTgg                                                           >  3:606612/1‑139 (MQ=255)                                                                                                    aCCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTg                                                <  3:578471/149‑1 (MQ=255)                                                                                                      cGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTgcg                                              >  3:353539/1‑149 (MQ=255)                                                                                                           ccTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTTGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTc                                         <  2:597437/149‑1 (MQ=255)                                                                                                            ctctGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCg                                        <  1:68811/149‑1 (MQ=255)                                                                                                             tctGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCgg                                       >  1:207430/1‑149 (MQ=255)                                                                                                               tGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGtt                                     >  2:158289/1‑149 (MQ=255)                                                                                                                ttcatggttcgccataccgacatagagtggcacaaTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCAt                                                                                                             >  2:213707‑M1/36‑76 (MQ=255)                                                                                                                ttcatggttcgccataccgacatagagtggcacaaTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCAt                                                                                                             <  1:213707‑M1/41‑1 (MQ=255)                                                                                                                    tgatgaAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGa                                <  3:793336/149‑1 (MQ=255)                                                                                                                       tgaAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTTGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGacac                             >  3:528047/1‑149 (MQ=255)                                                                                                                       tgaAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGacac                             <  1:842407/149‑1 (MQ=255)                                                                                                                        gaAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACAcc                            >  4:715616/1‑149 (MQ=255)                                                                                                                              ggCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCgg                      <  3:817311/149‑1 (MQ=255)                                                                                                                                cATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGa                    <  2:702565/149‑1 (MQ=255)                                                                                                                                    tACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCt                <  4:555561/149‑1 (MQ=255)                                                                                                                                     aCCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTg               <  4:574332/149‑1 (MQ=255)                                                                                                                                      ccGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTgg              >  1:457383/1‑149 (MQ=255)                                                                                                                                       cgcgATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGc             <  3:732353/149‑1 (MQ=255)                                                                                                                                             tACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGa       <  2:335422/149‑1 (MQ=255)                                                                                                                                             tACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGa       <  1:24428/149‑1 (MQ=255)                                                                                                                                              aCGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCAt                                                                              <  4:712633/77‑1 (MQ=255)                                                                                                                                              aCGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCAt                                                                              >  3:712633/1‑77 (MQ=255)                                                                                                                                                 aCTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTAc   <  1:445858/149‑1 (MQ=255)                                                                                                                                                  cTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACg  >  3:822521/1‑149 (MQ=255)                                                                                                                                                  cTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACg  >  2:557574/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    AACTGGAGAGCGCCATCGCCAACATCTCGGCGCTGAAAGGCGAGAAAGAGACAGAAGCACTGTTCCGCCAGAAAGCCAGTGCCCGTCGCGATGCGGTAAACCGTTACCTCTGGGATGATGAAAACGGCATCTACCGCGATTACGACTGGCGACGCGAACAACTGGCGCTGTTTTCCGCTGCCGCCATTGTGCCACTCTATGTCGGTATGGCGAACCATGAACAGGCCGATCGTCTGGCAAACGCCGTGCGCAGTCGGTTACTGACACCTGGCGGGATTCTGGCAAGCGAGTACG  >  NC_000913/3670661‑3670954

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