New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4576863 =	53 (1.130)	3 (0.070)	3/244	NT	5.8%	intergenic (+8/+49)	yjiV/mcrC	putative uncharacterized protein YjiV/5‑methylcytosine‑specific restriction enzyme subunit McrC
	?	NC_000913	4576890 =	48 (1.110)					intergenic (+35/+22)	yjiV/mcrC	putative uncharacterized protein YjiV/5‑methylcytosine‑specific restriction enzyme subunit McrC

GCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAA  >  NC_000913/4576723‑4577004                                                                                                                                             |                                                                                                                                              gaaacatcgttataaaattaatggcttcgatattggcttgtgtaccgtcaatttaggtcaggaatggccgtgtatacatcaagaattactcgacattttcgatgaatatctcaaataagcggtaagttttatagaagctaCCg                                                                                                                                             >  2:373213‑M1/141‑143 (MQ=255)  cATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCGGCAAGGAAGGGAATGTCGCCAACAGCGAAGGGAGTTGGGCAACGGATGTGCTGGGGGAGGGGGTCGCCTCCAGATGATGAGCGGc                                                                                                                                            >  1:816090/1‑143 (MQ=255)   aTTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGCTGAGCCGCt                                                                                                                                           >  2:703531/1‑143 (MQ=255)      atcgttataaaattaatggcttcgatattggcttgtgtaccgtcaatttaggtcaggaatggccgtgtatacatcaagaattactcgacattttcgatgaatatctcaaataagcggtaagttttatagaaaataCCGCTccc                                                                                                                                        <  2:761218‑M1/8‑1 (MQ=255)       cgttGACTGCCGCGATGACGTTTCTTTCCGGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCg                                                                                                                                       <  2:767689/141‑1 (MQ=255)        cgttataaacgtaatggcttcgatattggcttgtgtaccgtcaatttaggtcaggaatggccgtgtatacatcaagaattactcgacattttcgatgaatatctcaaataagcggtaagttttatagaaaataCCGCTCCCGa                                                                                                                                      >  2:551598‑M1/134‑143 (MQ=255)                   cgATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCgg                                                                                                                           <  2:248445/143‑1 (MQ=255)                     aTGACGTTCCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGa                                                                                                                         <  1:627948/143‑1 (MQ=255)                         cGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCgg                                                                                                                     >  2:156854/1‑143 (MQ=255)                            ttctttcCAGTTAGTCACAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTAt                                                                                                                  <  1:798759/143‑1 (MQ=255)                            ttctttcCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTAt                                                                                                                  <  1:158441/143‑1 (MQ=255)                            ttcgttgCAGTTAGTCAATAGGGCTGAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTAt                                                                                                                  <  1:320714/143‑1 (MQ=255)                               tttcCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGCGCGGTCtttt                                                                                                               >  1:137453/1‑143 (MQ=255)                                 tcCAGTTGGTCAAAATGGCTCAGGCCCGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCt                                                                                                             <  1:393011/143‑1 (MQ=255)                                 tcCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCt                                                                                                             >  2:547915/1‑143 (MQ=255)                                    aGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGACAACGGATGTGCTGGTGGAGGTGATCGCCTCCGTATGATGAGCCGCTCCCGATGTGGTGTCCGGGGCCGCAtttgcgatc                                                                                                          >  1:11304/1‑137 (MQ=255)                                       tAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATaaaa                                                                                                       <  1:240981/143‑1 (MQ=255)                                                 ggcggtggCCGGCTAGGAGGGGAGTGGCGGTAACAGGGGAGAGAGTGGGGCGAGGGATGTGCTGGTGGAGGTGGTCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCtt                                                                                             <  2:783486/137‑1 (MQ=255)                                                    tCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCCCCTCCTGATGCTGAGCCGCCCCCGATGTGGTGTCGGGAGCGGTATTTTCTATACAACCTCCCGCTTAtt                                                                                          >  2:628776/1‑143 (MQ=255)                                                       ggCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTga                                                                                       <  1:453558/143‑1 (MQ=255)                                                         ccTGCAAGGAAGCGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTgaga                                                                                     <  1:184869/143‑1 (MQ=255)                                                           tGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGCGata                                                                                   >  1:156029/1‑143 (MQ=255)                                                           tGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGata                                                                                   >  1:12106/1‑143 (MQ=255)                                                            gCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGatat                                                                                  <  1:10355/143‑1 (MQ=255)                                                                 gaaggGGATGGCGGCAAGAGCGAAGAGGGGTGGGCAACGGAGGTGCTGGTGGAGGTGATCGCCTCCGGATGAGGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATc                                                                             <  1:420098/143‑1 (MQ=255)                                                                 gaaggGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATc                                                                             <  1:881783/143‑1 (MQ=255)                                                                   gccgtgtatacatcaagaattactcgacattttcgatgaatatctcaaataagcggtaagttttatagaaaataCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGa                                                                           >  1:264606‑M1/75‑143 (MQ=255)                                                                    gggAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGCTGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAGCACTTACCGCTTATTGGGGATAGTCATCgga                                                                          >  2:841097/1‑141 (MQ=255)                                                                       gttGTCGCGAGCGGCGTAGATACGTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTGCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGGGCGGGATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAAt                                                                       <  2:748610/141‑1 (MQ=255)                                                                          gTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTc                                                                    <  1:684266/143‑1 (MQ=255)                                                                           tCGCCAACAGCGAAGAGAGTTGGGCCACTGCTGTGCTGGTGGAGGTGAGCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCg                                                                   <  2:8422/143‑1 (MQ=255)                                                                            cGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGa                                                                  <  1:761110/143‑1 (MQ=255)                                                                            cGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGa                                                                  <  1:331066/143‑1 (MQ=255)                                                                             gCCAACCGCGAAGAGAGCTGGGCAACGGGTGTGGTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAg                                                                 <  1:50215/143‑1 (MQ=255)                                                                              ccAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGCGTCGGGAGCGGTATTTTCTATAAAACTTACCGCGTATTTGAGATATTCATCGAAAATGGCGAGt                                                                >  2:576248/1‑143 (MQ=255)                                                                                agccgcgcAGATAGGTGGGCAACGGATTCGCTGGTTGACGTGAGCGCCTCCTGATGATGTGCCGCTCCCGATTCGGTAACAAGAGCTGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTaa                                                              <  2:604165/139‑1 (MQ=255)                                                                                aaCAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGCGTaa                                                              >  2:516363/1‑143 (MQ=255)                                                                                aaCAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTaa                                                              >  1:72781/1‑143 (MQ=255)                                                                                   aGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGCGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGCCGGGCGCGGTATTTTCTATAAACCTTACAGCTTAGTTGAGAAATACATAGAAAATGGCGacacatcc                                                           >  2:890894/1‑136 (MQ=255)                                                                                       aaGGGAGTTGGGCGACGGATGGGCTGGTGGAGGTGATCGGCTGCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGa                                                       <  2:687170/143‑1 (MQ=255)                                                                                       aaGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGa                                                       <  2:801087/143‑1 (MQ=255)                                                                                                  gCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACgg                                            >  1:387254/1‑143 (MQ=255)                                                                                                   cAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGc                                           >  1:453551/1‑143 (MQ=255)                                                                                                         aatatctcaaataagcggtaagttttatagaaaataCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTcc                                     <  1:570191‑M1/107‑1 (MQ=255)                                                                                                          tgtgCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCt                                    <  1:756296/143‑1 (MQ=255)                                                                                                            tgCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTACTTTCTATCAAACTTACCGCTTATTTGAGCTATTCATCGAAAATGTCGAGTAATTCTTGCTGTATACAGCGCCATTCCtta                                  >  2:42500/1‑141 (MQ=255)                                                                                                             gCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGAc                                 <  1:601145/143‑1 (MQ=255)                                                                                                              cTGGTGGAGGTGATCGCCTCCTGATGATGCGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTGATTTGAGATAGTCATCGAAAAACTCCGCTAAGACTCGATGTATACACGGCCATTCCTGAcc                                >  2:189559/1‑143 (MQ=255)                                                                                                              cTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCGAGATGTGGTGGCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATAGAAAATGTCGAGTACTTCTTGATGGATACACGGCCATTCCTGccc                                >  2:749643/1‑143 (MQ=255)                                                                                                                  tggAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAGCTTACCGCTTATTTGAGATATTCATCGAAAATGTAGAGTAATTCTTGATGTATACACGGCCCTTCCTGACCTaga                            >  1:300665/1‑141 (MQ=255)                                                                                                                      ggTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGa                        <  2:672797/143‑1 (MQ=255)                                                                                                                       gcggtaagttttatagaaaataCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGAc                       >  2:570191‑M1/23‑143 (MQ=255)                                                                                                                       gTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGAc                       >  1:492539/1‑143 (MQ=255)                                                                                                                         gATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACgg                     >  1:108738/1‑143 (MQ=255)                                                                                                                          aTCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGt                    <  1:424357/143‑1 (MQ=255)                                                                                                                            cGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGGTGTATACCCGGCCATTCCTGAGCTAGATTTACGGTac                  >  2:684466/1‑143 (MQ=255)                                                                                                                             gCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACCGCCATTCCTGCCCTAAATTGACGGTaca                 >  2:184462/1‑143 (MQ=255)                                                                                                                               ctcctGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTAcaca               <  1:852849/143‑1 (MQ=255)                                                                                                                               ctcctGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTAcaca               <  1:638055/143‑1 (MQ=255)                                                                                                                                 gcgcagcatggGCCGCTCCCGATGTGGTGTCTGGTCTGTTATTTTCTATAAAACTTACCGCTTATTTGAGATACTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAg             <  2:405614/136‑1 (MQ=255)                                                                                                                                  ctgctgATGCGCCGCTGCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGATATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGc            <  2:71321/143‑1 (MQ=255)                                                                                                                                    gatgatgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCa          >  2:142625/1‑143 (MQ=255)                                                                                                                                        atgaGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAatat      >  2:198207/1‑143 (MQ=255)                                                                                                                                            gCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCGTCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGaa  <  1:406884/143‑1 (MQ=255)                                                                                                                                             |                                                                                                                                              GCATTACTTTGATTTCCGCGATGACGTTTCTTTCCAGTTAGTCAAAATGGCTCAGGCCTGCAAGGAAGGGAATGTCGCCAACAGCGAAGAGAGTTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGAGCCGCTCCCGATGTGGTGTCGGGAGCGGTATTTTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGTCGAGTAATTCTTGATGTATACACGGCCATTCCTGACCTAAATTGACGGTACACAAGCCAATATCGAA  >  NC_000913/4576723‑4577004

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