New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3653230	NA (NA)	46 (0.830)	37/276	0.4	92.0%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	= 3766089	4 (0.070)					coding (3907/4134 nt)	rhsA	rhs element protein RhsA

GTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCG  >  NC_000913/3653107‑3653368                                                                                                                            |                                                                                                                                           gTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAttt                                                                                                                            >  2:757392/1‑140 (MQ=37)  tCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAtttt                                                                                                                           >  2:943010/1‑140 (MQ=38)         gcgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACAt                                                                                                                    <  1:863153/140‑1 (MQ=255)         gcgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACAt                                                                                                                    <  1:847027/140‑1 (MQ=255)          cgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATa                                                                                                                   <  1:42438/140‑1 (MQ=255)              tACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGa                                                                                                               <  2:879059/140‑1 (MQ=255)                ctgctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAac                                                                                                             >  1:453050/1‑140 (MQ=255)                      aaTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTa                                                                                                       >  2:909409/1‑140 (MQ=255)                      aaTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTa                                                                                                       <  1:757392/140‑1 (MQ=255)                           aCTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAccccc                                                                                                  >  2:290457/1‑140 (MQ=255)                           aCTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAccccc                                                                                                  <  1:1369/140‑1 (MQ=255)                              gTCGGCGAAGGTAAGTTGACGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaaa                                                                                               <  2:635668/140‑1 (MQ=255)                                  gCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATtat                                                                                           >  1:547040/1‑140 (MQ=255)                                        gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAg                                                                                     <  1:943010/140‑1 (MQ=255)                                        gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAg                                                                                     <  1:796151/140‑1 (MQ=255)                                                  gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAata                                                                           >  1:253029/1‑140 (MQ=255)                                                       atgatgAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTg                                                                      >  2:764841/1‑140 (MQ=255)                                                            gaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACCCCCTACCCCCAAATTATATAAAGATGGACAATATCTTGacaaa                                                                 <  1:810393/140‑1 (MQ=255)                                                            gaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGacaaa                                                                 <  1:764841/140‑1 (MQ=255)                                                                   gTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCg                                                          >  1:719895/1‑140 (MQ=255)                                                                    ttCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGa                                                         >  2:250697/1‑140 (MQ=255)                                                                    ttCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGa                                                         >  2:264235/1‑140 (MQ=255)                                                                         tGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGt                                                     <  1:715673/139‑1 (MQ=255)                                                                         tGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGt                                                     >  2:715673/1‑139 (MQ=255)                                                                              ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTAttt                                               >  2:641609/1‑140 (MQ=255)                                                                                  aTGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCaa                                           <  1:641609/140‑1 (MQ=255)                                                                                  aTGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCaa                                           <  2:453050/140‑1 (MQ=255)                                                                                    gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACa                                         <  1:909409/140‑1 (MQ=255)                                                                                           aTGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTaa                                  <  1:250697/140‑1 (MQ=255)                                                                                                     atCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGa                        >  2:557406/1‑140 (MQ=255)                                                                                                          ggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAAcc                   >  1:323989/1‑140 (MQ=255)                                                                                                                 ttCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACg            >  1:772892/1‑140 (MQ=255)                                                                                                                       ccttccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCaaa      <  1:369866/140‑1 (MQ=255)                                                                                                                           ccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCg  >  2:878399/1‑140 (MQ=255)                                                                                                                            |                                                                                                                                           GTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGCG  >  NC_000913/3653107‑3653368

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