New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	15 (0.370)	14 (0.370)	13/448	1.2	59.6%	coding (305/630 nt)	ycfK	protein StfP
	?	NC_000913	= 1209602	5 (0.130)					intergenic (‑483/‑83)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

GGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATT‑CCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGG  >  NC_000913/1207558‑1208047                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    ggTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                                                                                                                                                                                                                  >  2:11662/1‑251 (MQ=255)       aCCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcat                                                                                                                                                                                                                                             <  2:314008‑M1/250‑9 (MQ=255)              gTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                                                                                                                                                                                                                      >  1:152150/1‑250 (MQ=255)                 aCCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagca                                                                                                                                                                                                                                  <  2:3371/251‑1 (MQ=255)                   ccccAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgt                                                                                                                                                                                                                                >  1:179539‑M1/1‑230 (MQ=255)                      cAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGGTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                                                                                                                                                                                                                             <  1:190651‑M1/251‑25 (MQ=255)                       aGACCGGCACTCCTGCCACCGATCTGGATGATGACTACTTTGTCATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAAGCCGTTTGTCGATATCAAATCGGATGGCGCTGTGCAAACGTCTCTCGAAAACCGTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccg                                                                                                                                                                                                                            <  1:235753‑M1/251‑26 (MQ=255)                            ttttcagtcccggcctttcggactttgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAt                                                                                                                                                                                                                        >  1:91493‑M1/206‑250 (MQ=255)                             tttcagtcacggcctttcggactttgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaccgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAtt                                                                                                                                                                                                                       <  2:193416‑M1/46‑1 (MQ=255)                                 agtcccggcctttcggactttgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAc                                                                                                                                                                                                                  <  1:98984‑M1/51‑1 (MQ=255)                                                 gatgatGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcc                                                                                                                                                                                                  >  1:417450‑M1/1‑200 (MQ=255)                                                   ctttgtgccgcttatcttgtcttcaaaaatcagatcgcatcctgcccagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccattgattgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt                                                                                                                                                                                                 <  1:355981‑M1/68‑1 (MQ=255)                                                    gatgaCTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtat                                                                                                                                                                                               <  1:367511/251‑1 (MQ=255)                                                      gtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtata                                                                                                                                                                                              >  2:98984‑M1/180‑250 (MQ=255)                                                      gtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttatggtcatttgttgatacgcgtacatagccaataagcatgttaaatccacctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGATTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                                                                                                                                                                             <  2:291131‑M1/72‑1 (MQ=255)                                                          tactTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTTGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtttgtgtcccggttccgtggccttcagccactccgccaacaggctg                                                                                                                                                                                         >  1:151654‑M1/1‑191 (MQ=255)                                                           actTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt                                                                                                                                                                                         <  1:56495/250‑1 (MQ=255)                                                              ttGACATGTTGCCGGAGGAACTTTGCCGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGGAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCCAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctg                                                                                                                                                                                     >  2:249741‑M1/1‑187 (MQ=255)                                                                tccatGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggcgga                                                                                                                                                                                    <  2:90556‑M1/248‑66 (MQ=255)                                                                   tgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                                                                                                                                                                                <  2:341763‑M1/85‑1 (MQ=255)                                                                        gCAGGAGGATCTTTGCAGCGGTGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATAGCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaat                                                                                                                                                                            <  1:409110‑M1/250‑74 (MQ=255)                                                                          aggaggAACTTTGCAGCGTGGTGGAGTCATCCGGTGCCCGCCTGGAGAAGGGGCGGCACGACCAGTTTCTTACCGCCATTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTTGCGATATTAAATCGGATGGCACTTTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                                                                                                                                         <  2:94378/251‑1 (MQ=255)                                                                                                       tCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                                                                                                                                                                                                                     <  2:143799‑M1/178‑33 (MQ=255)                                                                                                       tCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGAGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattccctgttggtgtcccggttccgt                                                                                                                                                                                                                     >  1:143799‑M1/1‑146 (MQ=255)                                                                                                       gcgcattacgttgtagatctgtgttatggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTg                                                                                                                                            >  1:291131‑M1/131‑251 (MQ=255)                                                                                                           gTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAAt                                                                                                                                        <  1:11662/251‑1 (MQ=255)                                                                                                                    gggAGAAGGTGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATAAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGGAATTAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                                                                                                                               <  1:262416/250‑1 (MQ=255)                                                                                                                          aGGGGCGGCACGACCAGTTTCTTACCGCTCTTCGCGCGCTGCTGGTAAGCCGCAAGAATCCGTTTGTCGATATCAAAGCGGATGGCACTGTGCAAACGGCGCTCGAAAACCTTGGTTTGGGATAAGGAGCAAAACTCAATGCAGCATCGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAGCAGTAATCATTCAGTGGGGTTGGAATGCCGCATAAGCATCt                                                                                                                         <  2:333463/251‑1 (MQ=255)                                                                                                                                         aGTTACTTACAGCACTTCGCGCGCTGCTGTTACGCCGCAAGAATCCGTTTGGCGATCTCAAATCGGATGCCACTGTGCAAACGGCTCTCGACAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATACGGAAGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGAt                                                                                                          >  2:128107/1‑251 (MQ=255)                                                                                                                                                   cGCACTTCGCGCGCTGCGGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCCCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaaggcaatggtgccgctgtttctgctgaagaatacccggaactggcaaaggctgatccgacaaataaattgcctgatttacgtg                                                                                                <  2:133762‑M1/251‑150 (MQ=255)                                                                                                                                                                       tAAGCCGCAAGAATCCCCTTCGCGATATCAAGTCGGAAGGCACTGTGCAAAAGGCTCTCGATACCCTTGGATTGGGAGACGGAGCCAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGCTTCATAGCTATACAGGTTATGATTGGTGACATTGAGCAATCCGTAATCATTCAGTGGGAGTGGAATGCCGCAAACGCATCGGCCTCTGGGGGGGCTGGAAATACAGTTGTCTT‑CCCGTATGCGttt                                                                           >  2:203959/1‑251 (MQ=255)                                                                                                                                                                             gCAAGAATACGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCACGAAAACCGTGGTTTGGGGGAAGGAGCAAAACTCAAGGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAGTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATCCATTTGTATT‑CCCGGTTGCGTTTaataat                                                                     <  2:286484/251‑1 (MQ=255)                                                                                                                                                                                          tGGCGATAGCAAATCGGATGGCACTGTGCAAACGGCGCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgcgttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgttgtgagtttatgcgtggctgggatgacggtcgggggatt                                                          <  2:179539‑M1/250‑188 (MQ=255)                                                                                                                                                                                           ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCGCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCCGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGGTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGAGGGAAATAAGGTTG‑ATTACACTGTTGCGTTTCAAAATGCTTGTGTTGCCgg                                                       >  2:235081/1‑250 (MQ=255)                                                                                                                                                                                             aggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATT‑CCCGGTTGCGTTTAATAATGCCTGTGTTGCCgttg                                                     <  2:320699‑M1/207‑1 (MQ=255)                                                                                                                                                                                                               cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATT‑CCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATg                                   <  1:293817/251‑1 (MQ=255)                                                                                                                                                                                                                    ggcAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatgtacgtggtgagtttattcgtggctgggatgactggagggggattgatacaggtcgctctattttaagtatt                               <  1:249999‑M1/250‑215 (MQ=255)                                                                                                                                                                                                                                  ataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGtt                                                                                              <  2:21765‑M1/167‑1 (MQ=255)                                                                                                                                                                                                                                  ataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGtt                                                                                              >  1:21765‑M1/8‑174 (MQ=255)                                                                                                                                                                                                                                        gTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAAGCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTTCCTCTGGGGGGGATGGAAATACAGTTGTATT‑CCCGGTTGCGTTTAATAATGCCTGTGTGGCCGTTGTTGCAAATTATGACAATGGCAGCGCACCTATCAATGCAGTGGc          >  2:246954‑M1/2‑251 (MQ=255)                                                                                                                                                                                                                                                gggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATT‑CCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggggg  >  1:267681/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    GGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATT‑CCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGG  >  NC_000913/1207558‑1208047

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.