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

GGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGC  >  AE015451/604843‑605294                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     ggTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCgaacacggtcac                                                                                                                                                                                                                          <  1:536616‑M1/236‑13 (MQ=255)          ctctgctgTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGa                                                                                                                                                                                                                                <  1:717377/219‑1 (MQ=255)                       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)                                                    ggCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACCCCACGCGCGGCCCGt                                                                                                                                                                                                                                   >  1:108985/1‑176 (MQ=255)                                                               ccAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCg                                                                                                                                                            <  1:268177/236‑1 (MQ=255)                                                                      ccTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGa                                                                                                                                                                                                                                >  2:381255/1‑161 (MQ=255)                                                                                          cTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCaccgtgaccg                                                                                                                                                                                                                 >  1:214076/1‑156 (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)                                                                                                                                                 ggcggcCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCgaacacggtcacggtgatgctcaacgggccgc                                                                                                                                                                                                      <  2:67426‑M1/112‑33 (MQ=255)                                                                                                                                                 ggcggcCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGGGACACCACGCGCGGCCgaacacggtcacggtgatgctcaacgggccgc                                                                                                                                                                                                      >  1:67426‑M1/1‑80 (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)                                                                                                                                                                                                                          cgcgGCCgaacacggtcacggtgatgctcaacgggccgcgcgtggtgtcaccaccgatcaggctcatgcggcaacggtgggccatgcgaccaaggccgtcggcataggccgccagccagtcagggccaacgtcaggcagggtcagtgcgagggtgaagccgattgcagtcgcgcccatggctgccaggtcactggccgccacagccaacgaacgctggccaagcagcaacggatc  <  1:655111‑M1/235‑229 (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)                                                                                                                                                                                                                                cgaacacggtcacggtgatgctcaacggaccgcgcgtggggtcaccaccgatcaggctcatgcggcaacggtgggccatgcgaccaaggccgtcggg                                                                                                                                      >  2:596134‑M1/1‑1 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     GGTGCACTTTCCTGCTGTTTGCGATCCGTTGCTGCTTGGCCAGCGTTCGTTGGCTGTGGCGGCCAGTGACCTGGCAGCCATGGGCGCGACTGCAATCGGCTTCACCCTCGCACTGACCCTGCCTGACGTTGGCCCTGACTGGCTGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATGCCGCAGCTCACCCTCGGGGCTCTGCTGC  >  AE015451/604843‑605294

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