New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	28 (0.540)	23 (0.470)	22/456	1.0	51.6%	coding (305/630 nt)	ycfK	protein StfP
	?	NC_000913	= 1209602	17 (0.350)					intergenic (‑483/‑83)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

CCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGT  >  NC_000913/1209357‑1209849                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        ccATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTACAATCGAGCGACCTGTATCAATCCCCCGCCCGTCATCCAAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTtc                                                                                                                                                                                                                                                    >  1:66607/1‑251 (MQ=255)   aTGAGCATGATACTAAGTTGCATACCACTGACTACTTAAAATCGCGAGACTTGTATCAATACCCAGCCCGTCATCCAAGACACGAATAAACTCACCACGTAAATCAGGCAATTTATTAATCGCATAAGCCTTTGACAGTTCCGGGTATTCTTCAGCAGAAAACGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGAGGTGGCTGAAGGCCACGGAACAGGGACACCAACCGGTAATGCCGAGCCTtctc                                                                                                                                                                                                                                                  >  2:343031/1‑251 (MQ=255)     gAGCATGCACCTCCGTTGCATACCCCTGAATACTTAAAATAGAGAGGACTGTATCAATCCCCCGCCCGTCATCCCAACCACGAATAAACTCACCACGTAAATCAGGCCCTTTATTTGTAGGCTAAGACTTTGCCAGATCAGGGTATTCTTCAGCAGCACAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCCCTGAACAGGGACACCAACAGGTAATGCTGAGACTTCTccc                                                                                                                                                                                                                                                >  2:455609/1‑251 (MQ=255)         aTGATCCTCCGTTGCATACCCCTGAATTCTTAAAATCGAGCGCCCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAACAAACTCACCACGTCAATCAGGCAATTTATTTGTAGGATACGCCTTTGCCCGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCCTTTCAGCCAGCCTGTTGGCGGAGTGGCTGCAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAAc                                                                                                                                                                                                                                            >  1:467000/1‑251 (MQ=255)         aTGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCATCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAAc                                                                                                                                                                                                                                            <  1:204589/251‑1 (MQ=255)                  cGTTGCATACCCCTGAATACTTAAAATATAGCGACCTGTATCAATCCCCCGCCCTTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtttt                                                                                                                                                                                                                                   >  1:422026‑M2/1‑229 (MQ=255)                     cgcctacccCTGAAGACTGACAATAGAGAGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACTAATAAACTCCCCCCGTAAATCAGTTAATTTATTTGCCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGACAGTCTTTTGGCGGAGTGGCTGAAGGCCACGGAACCTGGACACCAACAGGTAATGCAGAGccttccaccaaaccaaggttttcg                                                                                                                                                                                                                                 >  2:157‑M2/5‑226 (MQ=255)                        aTACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGACGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccacaccaaggttttcgag                                                                                                                                                                                                                               <  1:243179‑M2/249‑27 (MQ=255)                                 aaTACGTAAAATAGGACGTCCTGGATCAATCCCCCCCCCGGCCTCCCATCCACGAATAAACTCACCACGTAAATCAGGCAATTTGTTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAATGCGGCTCCATGGCAATTCACCCAGCCTGTTCGCTAAGTGGCTGAAGTCCACGGATCCGGGACATCAACAGGTAAGGCAGAGCCTTCTCCCATACCAACGTTTATGaagaggtagaaat                                                                                                                                                                                                                    <  2:467000/251‑12 (MQ=255)                                           aTAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccttc                                                                                                                                                                                                          >  2:180056‑M2/1‑204 (MQ=255)                                           aTAGAGCGACCTGTATCAATCCCCCCCCCGTCATCCCAGCCACGAATAACCTCACCACGTAAATCAGTCAATTTATTTGTCGGATAAGTCTTTGCCAGTTCCGGCTATGCTTCAGCAGAAAAAGCGGCTCCATTGCATTTCAGCCAGCCTTTTGGGGGAGTGTCTGAATGCCACGTAACCTGGACTCCAACAGGTAATTCAGGGgcttctccgaaaccaaggttttcgagagccgtttgcacagtgccatc                                                                                                                                                                                                          >  1:343031‑M2/1‑204 (MQ=255)                                              gagCGACCTGTATAAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgt                                                                                                                                                                                                       <  2:263718‑M2/251‑51 (MQ=255)                                                           cAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCATGCt                                                                                                                                                                                          >  2:85253/1‑251 (MQ=255)                                                               cccccGCCCGTCATACCAGCCACGACTAAACTCCCCACGTAATTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAACGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGCGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggacttccagagccgtttgcacagtgcaatccgttttgatatcgccaaacg                                                                                                                                                                                      <  2:216283‑M2/251‑68 (MQ=255)                                                                             cccAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAGAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                                                                                                        <  1:341731/251‑1 (MQ=255)                                                                                       aaTAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCCGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGAAACTGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt                                                                                                                                                              >  1:466097/1‑251 (MQ=255)                                                                                       aaTAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAGGAAAATGAAGAAATAACAAGCAAATGGCATCAGTCCTGCTTTTACCAGGGGGATTGAACATGCTTAtt                                                                                                                                                              <  2:445308/251‑1 (MQ=255)                                                                                                 ccacGTAAATCAGGCAATTTATTTGTTGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttcccccaaaccaaggttttcgagagccgtttgcacagagccatccgatttgatatcgccaaacggattcttgcggcttcacagcagcgcgcgcagt                                                                                                                                                      <  2:179928‑M2/249‑100 (MQ=255)                                                                                                     gTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggta                                                                                                                                                >  1:442099‑M2/1‑146 (MQ=255)                                                                                                             ggCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATAGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGACAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaaca                                                                                                                                         >  1:215221/1‑250 (MQ=255)                                                                                                              gCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactgg                                                                                                                                       <  1:333811‑M2/251‑115 (MQ=255)                                                                                                               cAATTTAGTTGTCTGATAAGCCTGTGCCAGTTCCGTGTATTCTTCAGCAGAAAAAGCGCCACCAGTGCATTTCAGCCAGCCTCTTGGCGGAGTGGCTGAAGGCCAAGGAACCGGGACACCAACAGGTAATGCAGAGCCGTCTCCCAAACCAACGTGTATGAAAATGAAGAAATAACAAGCAAATGGCATCATGCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                                                                                                                      <  2:466097/251‑1 (MQ=255)                                                                                                                       ttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGa                                                                                                                               >  1:535204/1‑250 (MQ=255)                                                                                                                         gTCGGAAAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGGAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAac                                                                                                                             <  2:139847/250‑1 (MQ=255)                                                                                                                           cGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcccct                                                                                                                          >  1:213864‑M2/1‑124 (MQ=255)                                                                                                                            ggATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGTATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                                                                                                                         <  1:426168/251‑1 (MQ=255)                                                                                                                            ggATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTGCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACTGAACCGGGACACCAACAGGTAATTCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTGATTGGCTATTTACGCGTATCAACAAATGACCAGAACACAg                                                                                                                         <  1:315570/251‑1 (MQ=255)                                                                                                                                               ccGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggat                                                                                                      >  2:333811‑M2/1‑104 (MQ=255)                                                                                                                                                            agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAgg                                                                                          >  1:245929/1‑250 (MQ=255)                                                                                                                                                                gAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatctccaaacggattcttgcggct                                                                                                                                                                         >  2:468231‑M2/1‑87 (MQ=255)                                                                                                                                                                gAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggct                                                                                                                                                                         <  1:468231‑M2/167‑81 (MQ=255)                                                                                                                                                                 aaaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGa                                                                                    <  1:56051/251‑1 (MQ=255)                                                                                                                                                                       cGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcc                                                                              <  1:448447‑M2/251‑172 (MQ=255)                                                                                                                                                                          cACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgccaagttcctcc                                                                           <  2:402880‑M2/251‑175 (MQ=255)                                                                                                                                                                                 gCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaac                                                                     >  1:402880‑M2/1‑70 (MQ=255)                                                                                                                                                                                       cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCTGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattggatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtcaa                                                              >  2:243179‑M2/1‑64 (MQ=255)                                                                                                                                                                                                          gTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGaaa                                           <  2:245929/251‑1 (MQ=255)                                                                                                                                                                                                                     ccACGGAACCGAGACACCAACAGGTAATGCAGAGccttctcccaatccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttccgcctgcaacatgtcaaagtagtcatcatccagatcggtggcaggcg                                >  1:263718‑M2/1‑34 (MQ=255)                                                                                                                                                                                                                         ggAACCGGGACACCAACAGTTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgttccgccccttctccagtctgtcaccggatgcctccacctcgctgcaaagttcctcctgcaacatttcaaagtagtcatcatccagatcggtggcaggcgtgcc                            >  1:455609‑M2/1‑30 (MQ=255)                                                                                                                                                                                                                           aaCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGaaaaa                           >  1:458889/1‑250 (MQ=255)                                                                                                                                                                                                                                            gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATc         <  2:458889/251‑1 (MQ=255)                                                                                                                                                                                                                                                  cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtcaaagtagtcatcatccagatcggtggcaggcgtgccggtctgggggttaccacgggtaaa    <  1:303452‑M2/250‑246 (MQ=255)                                                                                                                                                                                                                                                   agagCCTTCTACCCAACCAACCTTTATGAAAATGAAGAAATAAAAAGCAAATGGCATAATTCCTGCTTTTACCAGGGGGATCTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACAAGAACAAAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGAGGCACAAATTCCGAAAGGCAGTGACTGAAAACAATGCTCAGGCCATTAtaggaaggg  >  2:363310/1‑243 (MQ=255)                                                                                                                                                                                                                                                    gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtcaaagtagtcatcatccagatcggtggcaggcgtgccg                           >  2:445721‑M2/1‑3 (MQ=255)                                                                                                                                                                                                                                                    gagccttctcccaaaaaaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggatgcctccaccacgctgcaaagttcctcctgcaacatgtcaaagtagtcatcatccagatcggtggcaggcgtgccg                           <  1:445721‑M2/225‑223 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        CCATGAGCATGATCCTCCGTTGCATACCCCTGAATACTTAAAATAGAGCGACCTGTATCAATCCCCCGCCCGTCATCCCAGCCACGAATAAACTCACCACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGT  >  NC_000913/1209357‑1209849

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