New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AE015451	= 605066	28 (0.690)	3 (0.080)	3/322	NT	9.5%	coding (379/969 nt)	thiL	thiamin monophosphate kinase
	?	AE015451	= 605092	30 (0.770)					coding (405/969 nt)	thiL	thiamin monophosphate kinase

GATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCG  >  AE015451/604865‑605327                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            gATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCggccgcg                                                                                                                                                                                                                                      <  1:655111‑M2/235‑8 (MQ=255)    ccGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGc                                                                                                                                                                                                                                   <  1:500611/235‑1 (MQ=255)           tgctTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGc                                                                                                                                                                                                                           >  1:359002/1‑236 (MQ=255)                     gCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCGTgc                                                                                                                                                                                                                 >  1:543429/1‑236 (MQ=255)                                         ccAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCg                                                                                                                                                                                             <  1:268177/236‑1 (MQ=255)                                                                     tGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCggtggt                                                                                                                                                                     <  2:655110/232‑1 (MQ=255)                                                                                              aCCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTgcg                                                                                                                                                                                                                <  2:109581/164‑1 (MQ=255)                                                                                              aCCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGc                                                                                                                                         <  1:459486/235‑1 (MQ=255)                                                                                               cccTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGc                                                                                                                                              >  2:459485/1‑229 (MQ=255)                                                                                                cctgcctgACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTgcgc                                                                                                                                                                                                               <  1:492110/163‑1 (MQ=255)                                                                                                      tgACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGGGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGc                                                                                                                                      >  1:698174/1‑230 (MQ=255)                                                                                                                 cTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAgcgc                                                                                                                     >  1:476972/1‑236 (MQ=255)                                                                                                                       ggctggcGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGc                                                                                                                                                                                               >  2:669313/1‑156 (MQ=255)                                                                                                                       ggctggcGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGc                                                                                                                                                                                               <  1:669314/156‑1 (MQ=255)                                                                                                                                   gtGCCGGCGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCaa                                                                                                                                                          <  2:227289/180‑1 (MQ=255)                                                                                                                                     cccgACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACCCCACGCGCGGTCCGTTGAGCATCACCGTGACCGTGTTCg                                                                                                                                                                                                                                            >  2:596134‑M2/2‑96 (MQ=255)                                                                                                                                                 ggTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCgg                                                                                                                                                                                            >  2:239682/1‑133 (MQ=255)                                                                                                                                                            ccACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTAc                                                                           <  1:646779/235‑1 (MQ=255)                                                                                                                                                                  ttGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGctgctg                                                                                    <  2:315650/220‑1 (MQ=255)                                                                                                                                                                       gCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGtt                                                                                                                                                                    >  2:152292/1‑135 (MQ=255)                                                                                                                                                                        cATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTg                                                                         >  1:505568/1‑225 (MQ=255)                                                                                                                                                                        cATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGc                                                              >  1:647137/1‑236 (MQ=255)                                                                                                                                                                                     ggtggtGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCg                                                                                                                                                                          <  2:749925/115‑1 (MQ=255)                                                                                                                                                                                     ggtggtGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCg                                                                                                                                                                          >  1:749926/1‑115 (MQ=255)                                                                                                                                                                                     ggtggtGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCAccc                                                  <  1:332244/235‑1 (MQ=255)                                                                                                                                                                                               ccacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGGGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGGGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGctgctg                                                                                    >  1:151213/1‑191 (MQ=255)                                                                                                                                                                                                 acGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTg                                                                                                                                    <  2:108985/141‑1 (MQ=255)                                                                                                                                                                                                     gcgGCCCGTTGAGCATCACCGTGACCGTGTTCggccgcgcgtggtgtccccaccgatcaggctcatgcggcaacggtgggccatgcgaccaaggccgtcggcataggccgcc                                                                                                                                                             >  1:67426‑M2/1‑32 (MQ=255)                                                                                                                                                                                                     gcgGCCCGTTGAGCATCACCGTGACCGTGTTCggccgcgcgtggtgtcaccaccgatcaggctcatgcggcaacggtgggccatgcgaccaaggccgtcggcataggccgcc                                                                                                                                                             <  2:67426‑M2/112‑81 (MQ=255)                                                                                                                                                                                                       ggCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCgg                                              >  1:20401/1‑221 (MQ=255)                                                                                                                                                                                                           cGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCg                                                                                               >  1:496005/1‑168 (MQ=255)                                                                                                                                                                                                           cGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCg                                                                                               <  2:496004/168‑1 (MQ=255)                                                                                                                                                                                                                aGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATg                                                               <  2:359001/195‑1 (MQ=255)                                                                                                                                                                                                                     cgccgtgaccgtgTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGcc                                                                                                                                        <  2:318043/115‑1 (MQ=255)                                                                                                                                                                                                                         gtgaccgtgTTCggccgcgcgtggtgtcaccaccgatcaggctcatgcggcaacggtgggccatgcgaccaaggccgtcggcataggccgccagccagtcagggccaacgtcaggcagggtcagtgcgagggtgaagccgattgcagtcgcgcccatggctgccaggtcactggccgccacagccaacgaacgctggccaagcagcaacggatcgcaaacagcaggaaagtgcacc             <  1:536616‑M2/236‑225 (MQ=255)                                                                                                                                                                                                                          tgaccgtgTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGCTCTGTTGTGTGTCGGTGGTTCGCTGGGCaa                                                                                                                                                          >  2:240514/1‑94 (MQ=255)                                                                                                                                                                                                                                    cGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCt                                                                                                                               <  2:414286/111‑1 (MQ=255)                                                                                                                                                                                                                                    cGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCg  >  1:289395/1‑236 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            GATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGCGTGGCAAGGCCACTGCGGCCCTGGACATTTCCG  >  AE015451/604865‑605327

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