New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2288919 =	NA (NA)	32 (0.940)	31/484	0.2	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	= 4577768	0 (0.000)					coding (191/1047 nt)	mcrC	Type IV methyl‑directed restriction enzyme EcoKMcrBC subunit

TTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGG  >  NC_000913/4577523‑4578004                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             ttAACTGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGCGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATTGCAGGAATGATCTCGGTGTTAGgga                                                                                                                                                                                                                                           <  2:278509‑M2/249‑4 (MQ=255)                     cccGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaat                                                                                                                                                                                                                     <  1:344800‑M2/251‑26 (MQ=255)                           aaTTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGAGTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttct                                                                                                                                                                                                               >  2:119526‑M2/1‑220 (MQ=255)                                     aaaGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTGTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgact                                                                                                                                                                                                      >  2:211922‑M2/1‑210 (MQ=255)                                         tGAACGCGCTTCATATCTGATAGTTGAATTTAACTTTTCATGTTTAATTAATATGGCTCATGTGCTTATTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGAATATCAAAAGTACTGACGGTTTTCCCATGATTCAGCTGGAAGACGAGAATTCTTTTAGCACACTCTATTCGCCCTTTGATGCCAGGAATGGTCTCGGGATTAGggaaggcgcgaataagcggggaaattcttatcggctgactcagt                                                                                                                                                                                                  <  2:320146‑M2/250‑45 (MQ=255)                                              gAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttc                                                                                                                                                                                            <  1:303277‑M2/251‑51 (MQ=255)                                                                 tGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgcctcagtcatttcatttcttcatgtttgggcc                                                                                                                                                                         <  2:123312‑M2/251‑70 (MQ=255)                                                                      ttAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgatt                                                                                                                                                                     >  2:310005‑M2/1‑177 (MQ=255)                                                                                         taatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGCTTAAGATGTAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgcctt                                                                                                                                                 <  1:195238‑M2/251‑94 (MQ=255)                                                                                           atatGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttg                                                                                                                                                <  2:187549‑M2/250‑95 (MQ=255)                                                                                            tatGGCTAATGTGCTTTTTATAATTCGGTTAGCCCGCGTGTCTTAATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTATATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcatcatgtttgagccgattttttctcccgtaaatgccttga                                                                                                                                               <  1:322176‑M2/250‑96 (MQ=255)                                                                                                       tgCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctattt                                                                                                                                   <  1:178975‑M2/251‑108 (MQ=255)                                                                                                          tttttATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctattt                                                                                                                                   <  1:157324‑M2/248‑108 (MQ=255)                                                                                                             ttATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGGACTGACGGTTTTCCCATGATTCAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccg                                                                                                                             <  1:299073‑M2/251‑114 (MQ=255)                                                                                                                  aTTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAATAACTGACGGTTTTCCCATGATTAAGATGGAAGACGCGTATTGTTTTAGAAAACTCTATTCGCCCTTTGATGCCAGGACTGATCTCGGTGTTAGggaaggtgcgaataagcggggacattcttatcggctgactcagtcatttcatttcttcatgtttgagccgattttttctccggtaaatgacttgcatcagcgtatttagcacgtttct                                                                                                                        <  2:310594‑M2/251‑119 (MQ=255)                                                                                                                      ggTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgc                                                                                                                    <  1:264702‑M2/251‑123 (MQ=255)                                                                                                                         tAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccat                                                                                                                 <  1:158241‑M2/251‑126 (MQ=255)                                                                                                                             aagCGTGGCTTCATTAAGCATATCAAAAGTACTGACTGTGTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaa                                                                                                             >  2:400‑M2/2‑122 (MQ=255)                                                                                                                                cGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggc                                                                                                          <  2:278043‑M2/251‑133 (MQ=255)                                                                                                                                cGTGTCTTCATTAAGCATATCAAAAGTAATGCCCGTTTTCCCATGCTTAAGATCGATGCATAGCATTGTTTTCGCAAACTCTATACGCCCTATGATGCCAGGAATGATCTCGGTGTTAGggaagatgcgaataagcggggcaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaactgccttgaatcagcctatttagaccgtttcttcgccatttacggc                                                                                                          <  2:20320‑M2/251‑133 (MQ=255)                                                                                                                                        cATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgcttttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatccc                                                                                                  <  2:338335‑M2/251‑141 (MQ=255)                                                                                                                                                            cTGACGGTTTTCCCATGATTAAGATGGCAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcatcctatttagaccgtttcttcgccatttaaggcgttatccccagtttttagtgagatctct                                                                              <  2:259170‑M2/251‑161 (MQ=255)                                                                                                                                                                      tCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatccccagtttttagtgagatctctcccactgacg                                                                    <  2:131297‑M2/251‑171 (MQ=255)                                                                                                                                                                          aTGATTCAGATGGAAGACGAGTCTTGTTTTAGCAAACTCTATTCGCCCTTTGATACCAGGAATGATATCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactccgacatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatccccagtttttagtgagatctctcccactgacgtatc                                                                <  2:314321‑M2/251‑175 (MQ=255)                                                                                                                                                                             aTTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatccccagtttttagtgagatctctcccactgacgtatcatt                                                             <  1:252697‑M2/251‑178 (MQ=255)                                                                                                                                                                                aaGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttattcgccatttaaggcgttatccccagtttttagtgagatctctcccactgacgtatcatttg                                                           >  2:195238‑M2/1‑71 (MQ=255)                                                                                                                                                                                           cgcgTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatccccagtttttagtgagatctctcccactgacgtatcatttggtccgcccgaaa                                               >  2:345219‑M2/1‑60 (MQ=255)                                                                                                                                                                                                        aGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatccccagtttttagtgagatctctcccactgacgtatcatttggtccgcccgatacaggttggccag                                   >  2:297397‑M2/1‑47 (MQ=255)                                                                                                                                                                                                                                   aGGAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatctccagtttttagtgagatctctcccactgacgtatcatttggtccgcccgaaacaggttggccagcgtgaataacatcgccagttggttatcg       >  1:245241‑M2/1‑20 (MQ=255)                                                                                                                                                                                                                                     gAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatccccagtttttagtgagatctctcccactgacgtatcatttggtccgcccgaaacaggttggccagcgtga                              >  1:52581‑M2/1‑18 (MQ=255)                                                                                                                                                                                                                                     gAATGATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttctcggctgactcagtcatttcatttcttcatgtttgagccgattttttctcccgtaaatgccttgaatcagcctatttagaccgtttcttcgccatttaaggcgttatccccagtttttagtgagatctctcccactgacgtatcatttggtccgcccgaaacaggttggccagcgtga                              <  2:52581‑M2/226‑209 (MQ=255)                                                                                                                                                                                                                                         gATCTCGGTGTTAGggaaggtgcgaataagcggggaaattcttttcggcggactcagtcatttcagtgcttcatgtttgagccgattttttctcccgtaaatgccttgaaccagcctatttagaccgtttcttcgccatttaatgcgttatccccagtttttagtgtgatctctcccactgacgtaccatgtggtccgctctaaataggttggcgagcgtgtataacatcgccagttggtttccgttttt  >  2:132552‑M2/1‑14 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             TTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATATCAAGAAGATTGTTACCGGGTATGGCTTCAAGGTTTGCCTGCTTAATTTCCTGTAAATAACCCCATGCATAGGTAAGCATGTAATAGATATTACGGACAGGTATCACGGGCTGTTCCACTATGAGTCCCCTAATAATTTGTTGGTCCATTTCTGTTGTTTATAGG  >  NC_000913/4577523‑4578004

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