New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1207790 =	22 (0.430)	12 (0.250)	10/452	1.9	36.6%	coding (290/630 nt)	ycfK	protein StfP
	?	NC_000913	1209619 =	21 (0.440)					intergenic (‑500/‑66)	tfaE/pinE	putative tail fiber assembly protein TfaE/site‑specific DNA recombinase

AACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGC  >  NC_000913/1207555‑1208033                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    aaCGGTTTTTCCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGt                                                                                                                                                                                                                                                 <  2:89659/240‑1 (MQ=255) aaCGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGt                                                                                                                                                                                                                                                 >  1:89659/1‑240 (MQ=255)            cagTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTAATTTGCCATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTCCTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAACCGGCTCTCGGAAACCTTGGTTTGGGCGACGGCCCAAAACTc                                                                                                                                                                                                                            >  2:157426/3‑250 (MQ=17)                  gtcctgagcagttttttcagtcccggcctttcggactttgtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                                                                                                                                                                                                       <  1:372463‑M1/31‑1 (MQ=255)                   aaCCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                                                                                                                                                                                     <  1:186845/250‑1 (MQ=255)                    aCCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagca                                                                                                                                                                                                                   <  1:523426/251‑1 (MQ=255)                               ggCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAtt                                                                                                                                                                                                        >  1:486115/1‑251 (MQ=255)                               ggCACGCCAGCCACCGATCTGCATGATGAATACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCCTCCGGTGCCAGCATGGAGAAGGGGCGGCCCGGCCAGTTACTTCACGCACTTCGCGCGCTGATGTTAAGCCGCAAGAATCCGTTTGGAGCTATCAACTCGGATGGCACTGTGCAAACGGCTCTCGAACACCTTGGTTTGGGAGAAGGAGCACAACTCCATGAAGCAACGGCTACAtt                                                                                                                                                                                                        >  2:516536/1‑251 (MQ=255)                                  ccgCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtg                                                                                                                                                                                                     <  1:11403‑M1/250‑34 (MQ=255)                                              gattgtgaagatgaCTACTTTGACAAGTTCCAGGAGGAACGTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                                                                                                                                                                         <  2:548737/242‑1 (MQ=255)                                                 cggacttggtgccgcttatcttgtcttcaaaaatcagctcgcatcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCAt                                                                                                                                                                                      >  1:497383‑M1/188‑251 (MQ=255)                                                       tatgacgttgtgtacaggTTGCAGGAGGAACTTTGCAGCGGGGTGGAGGCCTCCGATGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTTCAAACGGCCCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATTCAGCAACGGCTACATTAGGACTCACCTGTTTCATAGCta                                                                                                                                                                                 <  2:376925/233‑1 (MQ=21)                                                       gatgaCTACTTGGCCAGGTGGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCATTTACTTACCGCCCTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCCCCGGTTTCATAGCta                                                                                                                                                                                 <  2:86162/250‑1 (MQ=21)                                                        atgaCTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCta                                                                                                                                                                                 >  1:376925/1‑249 (MQ=255)                                                                    tttgtcttcaaaaagcagctcgcatcctgcacagttcggcgaattacgttgtaaaactgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagccggactgatgccatttgcttgttatttcttcattatcataaacgTTGGTTTGGGAGAGGGTGCACAACTCAATGCAGCAAAGGCTACATTAGGACGCACCGGTTtgatagctatacaggtactgatt                                                                                                                                                                   <  2:268445‑M1/83‑12 (MQ=255)                                                                            cAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCCGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCAGCAAGAATCCGTTTGGCGATATCCAATAGGATGGCACTGTGCAAACGTCTCTCGAAAACCTTGGTTTGGGTGAAGGAGCAATACTCAATGCAGCAACGGCTAAATTAGGACGCACCGGTTTCATAGCTAAACGGGGTATTATGGGTGGTAt                                                                                                                                                           >  2:499284/1‑251 (MQ=255)                                                                            cAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCCAGAATCCGTTTGGCGATATCACATCGGATGGCACTGTGCAAACGGCTCTCGAAAACATTGGTTTGGGAGAAGGctctgcattacctgtaggtgtcccggttccgtggccttcggccactccgccaacaggctggctggaatgcaatgg                                                                                                                                                           >  2:253795‑M1/1‑176 (MQ=255)                                                                             aggaggAACTTGGTAGCGGTGTGGGGGCCTCCGGTGCCAGCCTTGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAAGCCGTTTGGCGATATCAAATCGGATGGCTCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcacggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggt                                                                                                                                                          <  2:26861‑M1/251‑77 (MQ=255)                                                                                aggaACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGCAGGctctgcatgacctgttggtgtcccggttccgcagccttcagcccctccgccaaacggctggctgaaatgcacaggtgcc                                                                                                                                                       >  2:5529‑M1/1‑172 (MQ=255)                                                                                       ttGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGCGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgcttttt                                                                                                                                                >  1:128393‑M1/1‑165 (MQ=255)                                                                                       ttGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATc                                                                                                                                                  <  1:367394/249‑1 (MQ=255)                                                                                        cgcatcctgcacaattcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGa                                                                                                                                                                     >  1:395474‑M1/149‑229 (MQ=255)                                                                                        cgcatcctgcacaattcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAATCAATGCAGCAACGGCTACATTAGGACGCACCGGtttcatagctataccggttatgaaactcggaagagcgtcgtggcg                                                                                                                                               <  2:395474‑M1/103‑23 (MQ=255)                                                                                                         aTCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGATGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctggagaataaccgg                                                                                                                              >  1:227354‑M1/1‑147 (MQ=255)                                                                                                            gcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaattatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGa                                                                                                                           >  1:564329‑M1/129‑251 (MQ=255)                                                                                                                          agaAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttat                                                                                                             <  1:222429‑M1/251‑122 (MQ=255)                                                                                                                           gaAGGGACGGCCCGACCAGTTCCTTACCGCAACTCGCGCGCTGTTGTTAAGCCGCAAGAATCCGGTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCAt                                                                                                            <  2:140525/251‑1 (MQ=255)                                                                                                                             aGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccg                                                                                                          <  1:353609‑M1/251‑125 (MQ=255)                                                                                                                              gtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCt                                                                                                          >  1:138573‑M1/111‑250 (MQ=255)                                                                                                                                      aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataa                                                                                                  >  1:98083‑M1/1‑118 (MQ=255)                                                                                                                                          ccAGTTTCTTACCTCACTTCGCGCGCTGCTGTTAAGCCGCAAGAAGCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGATGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCTCAAAAGCATCTGCCTCTggggggg                                                                                             <  2:275965/251‑1 (MQ=255)                                                                                                                                          ccAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgc                                                                                             <  1:5529‑M1/251‑138 (MQ=255)                                                                                                                                                   acgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAAt                                                                                     <  1:435013‑M1/161‑1 (MQ=255)                                                                                                                                                          cTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccgtaactggcaaaggattatccg                                                                                                          >  2:238548‑M1/1‑98 (MQ=255)                                                                                                                                                          cTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccg                                                                                                          <  1:238548‑M1/222‑125 (MQ=255)                                                                                                                                                                 cggtgccgttaagccGCAAGAATCCGTGTGGCGATATCAAATCGGATGGCACTGTGCATACGGCGCTCGAAAACCTTGGTTTGGGAGAAGTctcggcattaccttttggtgtcccggttccgtggccgtcagccactccgccaacaggctggctgaaatgcaatggtttcgctttttctgctgaagaatgcccggaactggcaaaggtttatccgacaaagaaattgccggatttacgtggtgagtttagt                                                                      <  2:101836‑M1/244‑161 (MQ=255)                                                                                                                                                                  gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCgg                                                                     >  1:299905/1‑251 (MQ=255)                                                                                                                                                                     gctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtg                                                                  >  1:138603‑M1/1‑87 (MQ=255)                                                                                                                                                                      gtgtTAAGCCGCAAGAATCCGTTTGGCGATATCAAAGCGGATGGCACTGTGCAAACGGCTCTCGAAAACCGTGGTTTGGGAGAAGGctctgcattacctgttggtggcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtgg                                                                 <  1:312801‑M1/250‑166 (MQ=255)                                                                                                                                                                       tgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCg                                                                >  1:66758/1‑251 (MQ=255)                                                                                                                                                                             gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATTGCACTGTGCAAACGGTTCTCGAAAGCCTTGGTTTGGGAGAAGGctctgcattacctgttggtttcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaatggtttagccgacaaataaattgcctgatttacgttgtgagtttattcttggctgggat                                                          <  2:523302‑M1/251‑173 (MQ=255)                                                                                                                                                                                          gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCtgt                                                 <  2:490813/247‑1 (MQ=255)                                                                                                                                                                                          gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCtgt                                                 >  1:490813/1‑247 (MQ=255)                                                                                                                                                                                                aggaatgatgccatttgcttgttatttcttcattttcagaaacgTTGGGTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGGGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTAGTCCCGGTTGCTTTTAATAATGCCTGTGTTGCCgttg                                       >  2:318421‑M1/45‑251 (MQ=255)                                                                                                                                                                                                 gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgattc                                      >  1:412301‑M1/1‑59 (MQ=255)                                                                                                                                                                                                    atgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCATGTGTTGCCGTTGTTGc                                   <  1:253739‑M1/211‑1 (MQ=255)                                                                                                                                                                                                      cAAAACGGAGGGCACTGTGCAAACGGCTCTCGAAAATCTTGGTTTGGGAGAAGGctctgcattaccttttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtc                                 <  2:327094‑M1/251‑198 (MQ=255)                                                                                                                                                                                                         aTCGGCTGGCACTGTGCAAACGGCTCTCGAAACCATTGGTTTGGGAGAAGGAGCAAACCTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGATATGATTGGTGGTATTCAGCAATCAGTAATAATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCATCTGGGGGGGAGGGAAATACAGTATTATACCCGGTTGGATTTAATATTACCGGTGTTGCCgtttgtgcaatgt                              >  2:414156/1‑241 (MQ=255)                                                                                                                                                                                                              aTGGCACGGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGCGAAGGctctgcattaactgttgttgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggattatacgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgaagggcgggggatgggtacaggtggctctttt                         >  2:328728‑M1/1‑46 (MQ=255)                                                                                                                                                                                                                ggtaCTGTGCAAACGGATCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCACCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACa                        <  1:264795/247‑1 (MQ=255)                                                                                                                                                                                                                      gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattttaagta                  <  1:253795‑M1/250‑213 (MQ=255)                                                                                                                                                                                                                        gCAAACTGCTCTCGAAAACCTGGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccggtgccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattttaagtatt                <  2:233402‑M1/250‑215 (MQ=255)                                                                                                                                                                                                                             cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatttacgtggtgagtttattcgtggctgggatgacgggcgggggattgatacaggtcgctctattttaagtattcagggg          <  2:98083‑M1/251‑221 (MQ=255)                                                                                                                                                                                                                              cattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAACGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGCTGGAAATACAGTTGAATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACATTGTCAACGCacca          <  2:465359‑M1/236‑2 (MQ=255)                                                                                                                                                                                                                                 tttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATATTGCCTGTGTTGCCGttgttgcaatttatgacatttttagctaacctatc       <  2:257352‑M1/239‑18 (MQ=255)                                                                                                                                                                                                                                       aaaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACAGCTACATTAGGACGCACCGGTTTCATAACTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGTAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGATGCCGTTGTTGCAAATTATGACAATGTCAGCGCaacattcaatgc  >  2:443424/1‑238 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    AACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGC  >  NC_000913/1207555‑1208033

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