New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	9 (0.230)	15 (0.410)	15/388	1.3	52.3%	coding (290/630 nt)	ycfK	protein StfP
	?	NC_000913	1209619 =	19 (0.520)					intergenic (‑500/‑66)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

TAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATA  >  NC_000913/1207511‑1208053                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                        cccccccacccccccccccccgttacagcttccataacaccaggcggtcactttcccttaatgtccttatcagctttttcagtcccggccgttcgtacgttcagccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                                                                                                                                                                                                                                     >  2:457787‑M1/280‑284 (MQ=255)          aaaCCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGACCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGA                                                                                                                                                                                                                                                             >  1:436810‑M1/1‑283 (MQ=255)            gccccagccgatccagcttccagacaaccagagtgtcacctgccgataatgtcctgagcagttttttcagtcccggcctttcggactttgtgccgattatcttgtcttcaaaaatcagctcgcatcatgcacagttgagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagcgaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTgggaaa                                                                                                                                                                                                                                                             <  1:573988‑M1/13‑3 (MQ=255)                             ttccagacaaccagagtgtcacctgccgataatgtcctgagcagttttttcagtcccggcctatcggactttgtgcagcttatcttgtcttcaaaaatcagcgcgcatcctgcacagatcagcgaattacgttgtagatctgtgttctggtcatttgttgatatgcgtacatagctaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTGTGGGAGAAGGAGCAAAACtcaataccg                                                                                                                                                                                                                                         <  1:380021‑M1/33‑5 (MQ=255)                                                                 ctgagcagttttttcagtcccggcctttcggactttgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctttgttctggtcgtttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatattcttgttatttcttcattttcataaacaTTGGTTTGGGAGCAGTAGCATAACTCAATgc                                                                                                                                                                                                                                           <  1:270313‑M1/31‑1 (MQ=255)                                                                         ccGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacct                                                                                                                                                                                                                                             >  1:283782‑M1/1‑223 (MQ=255)                                                                                                      gaCTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgacatgcaatggtgccgctgtttctgc                                                                                                                                                                >  1:317769‑M1/1‑194 (MQ=255)                                                                                                      cacannttagacaggtngcAGGAGGACGTTTGCGGCGTGGTGGCGGCATACGGGGCCGGCCGGGGGAAGAGGCGGCACGACCAGTTACTTACCGCACTCCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgc                                                                                                                                                                <  2:436810‑M1/275‑91 (MQ=255)                                                                                                           ctTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                                                                                                                                                                                                       >  1:651135/1‑240 (MQ=255)                                                                                                             tatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttaatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgatatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAg                                                                                                                                                         >  1:491598‑M1/172‑284 (MQ=255)                                                                                                                    gTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAG                                                                                                                                                                                                                                                              <  1:29435‑M1/176‑1 (MQ=255)                                                                                                                    gTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAG                                                                                                                                                                                                                                                              >  2:29435‑M1/1‑176 (MQ=255)                                                                                                                     cttcaaaaatcagctcgcagcctgcacagttcagcgcattacgttgtagatctgtgttctgttcatttgttgttacgcgtacatagccaataagcaagttaaatccccctggtaaaagccggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGgag                                                                                                                                                                                                                                                              <  1:622844‑M1/12‑1 (MQ=255)                                                                                                                      tGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGACAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggct                                                                                                                                                                                               >  1:124213‑M1/1‑178 (MQ=255)                                                                                                                                aCTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgacatgcaatggtgccgctttttctgctgaagaatacccgcaactggaaaagg                                                                                                                                      >  1:363216‑M1/1‑168 (MQ=255)                                                                                                                                     gcatcctgaacagttcagcgcattacgttttagttctgttttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                                                                                                                                 >  1:288553‑M1/148‑284 (MQ=255)                                                                                                                                                        gcattacgttgtagctctgtggtctggtcatttgttgatactcgtacatagccaataagcatgttaaatccccctgttaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGGTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCCTAGCTATACCGGTTATGAGTGGTGGTAt                                                                                                                                                                               <  1:621602‑M1/91‑1 (MQ=255)                                                                                                                                                             acgttgtagatctgtgttctggtcatttgttgttacgcgtacatagccattaagcaggttaaatcctcctggtaaaagcaggaataatgcaatttgcttgttgttttttcattttcataaacgTTGGTTTGGGAGAAGGAGCAACACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTGTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAGTGCGGCAAAAGCATCTGCCTCGGGGGGGGATGGTAAt                                                                                                         <  1:395394‑M1/161‑1 (MQ=255)                                                                                                                                                                  cTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaagaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtgg                                                                                                    <  1:202817‑M1/284‑151 (MQ=255)                                                                                                                                                                      agaAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGTCACTGTGCAAACGGCTCTCGAAGACCTTGGTCTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaataccaggaaatggccaagggttatacgacaaaaaaattgcgagatttacgggcgngn                                                                                                >  2:202817‑M1/1‑130 (MQ=255)                                                                                                                                                                         aGGCGCGGCACGAACAGTTACTTACCGCACTTCGCGTGCTGCTTTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCATACGGCTCTCGAAACCCTTGGTGTGGGAGAAGGAGGAAAACCCAATGCAGCAACGGCTACATTAGGAGGCACCGGTTTCATAGCTATACCGGTTATGAGtg                                                                                                                                                                                      <  1:611575/195‑1 (MQ=25)                                                                                                                                                                           gggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                                                                                                                                                                                                                                    <  1:369690‑M1/131‑7 (MQ=255)                                                                                                                                                                           gggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                                                                                                                                                                                                                                    >  2:369690‑M1/1‑125 (MQ=255)                                                                                                                                                                            gttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                                                                                         >  1:88100‑M1/109‑269 (MQ=255)                                                                                                                                                                            ggcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgaggttatt                                                                                          >  1:5794‑M1/1‑124 (MQ=255)                                                                                                                                                                                 ggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTa                                                                                                                                                                                <  2:292374‑M1/90‑1 (MQ=255)                                                                                                                                                                                 ggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTa                                                                                                                                                                                >  1:292374‑M1/104‑193 (MQ=255)                                                                                                                                                                                   tcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCCGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCGTTGGGGGTGGAATGCAGAAAATGCATCTgcccgtgggggggagggagagacatgtg                                                                                                  <  2:88100‑M1/168‑11 (MQ=255)                                                                                                                                                                                                acCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATGAAAGCGGATGGCACTGTTCAAACGGCTCTCGAAAACCTTGGATTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggagaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcg                                                                      <  1:479465‑M1/284‑181 (MQ=255)                                                                                                                                                                                                   gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATTTCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGGTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAAt                                                                                                                                                                >  1:183157/1‑191 (MQ=255)                                                                                                                                                                                                          gccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGATTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCgtt                                                            <  1:17434‑M1/206‑1 (MQ=255)                                                                                                                                                                                                                          ccGCAATAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgacatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattt                                            <  1:55532‑M1/284‑207 (MQ=255)                                                                                                                                                                                                                             cAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgataaaggtcgctctattttaa                                         >  1:9325‑M1/1‑75 (MQ=255)                                                                                                                                                                                                                             cAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAAt                                          >  1:179559/1‑283 (MQ=255)                                                                                                                                                                                                                                               aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                                                                                                                                                                                                                    >  1:392311‑M1/42‑79 (MQ=255)                                                                                                                                                                                                                                               aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                                                                                                                                                                                                                    <  2:392311‑M1/38‑1 (MQ=255)                                                                                                                                                                                                                                                    tgccatttgcttgttatttcttcattttcataaacgTTGTTTTGGGAGAAGGAGCAAAACTCAATgcagc                                                                                                                                                                                                                                        >  1:351574‑M1/37‑70 (MQ=255)                                                                                                                                                                                                                                                    tgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                                                                                                                                                                                                                        <  2:351574‑M1/34‑1 (MQ=255)                                                                                                                                                                                                                                                      tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                                                                                                                                                                                   <  1:380102/73‑1 (MQ=255)                                                                                                                                                                                                                                                      tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                                                                                                                                                                                   >  2:380102/1‑73 (MQ=255)                                                                                                                                                                                                                                                             gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                                                                                         <  1:464854/188‑1 (MQ=255)                                                                                                                                                                                                                                                              cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattttaagtattcaggggtatgcaacggaggatcatgct         <  1:629904‑M1/283‑242 (MQ=255)                                                                                                                                                                                                                                                                    gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcccacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaatagcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggcttgagaacggttgctcttttttaagttttcaggggtatgcaacggaggagctgttgcttatac  >  1:316981‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                                                                    gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                                                                                         <  2:358891/181‑1 (MQ=255)                                                                                                                                                                                                                                                                    gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                                                                                         >  1:358891/1‑181 (MQ=255)                                                                                                                                                                                                                                                                         cGGCTCTCGAAAACCTTGGCTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccacgccgccaacaggctgactgaaatgcaatggtgccgctttttctactgaagaatacccgaaactggcaaagactaatccgacaaataaattgcctgttttacgtggtgagtttattcgtggctgggatgacgggcgggggagt                                                               >  1:324054‑M1/1‑31 (MQ=255)                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                        TAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATA  >  NC_000913/1207511‑1208053

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