New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3653230	NA (NA)	54 (1.110)	38/276	0.2	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	= 3766089	0 (0.000)					coding (3907/4134 nt)	rhsA	rhs element protein RhsA

TTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAG  >  NC_000913/3653102‑3653366                                                                                                                                 |                                                                                                                                          ttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGacac                                                                                                                                <  2:319890/140‑1 (MQ=33)  tCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACt                                                                                                                               >  1:734360/1‑140 (MQ=33)     ggTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAtt                                                                                                                            >  1:669429/1‑140 (MQ=37)     ggTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAtt                                                                                                                            <  1:206579/140‑1 (MQ=37)            cGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTAc                                                                                                                     <  1:224982/140‑1 (MQ=255)            cGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTAc                                                                                                                     <  2:716140/140‑1 (MQ=255)                gcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATaa                                                                                                                 >  1:300190/1‑140 (MQ=255)                   tACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGa                                                                                                              <  2:352177/140‑1 (MQ=255)                   tACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGa                                                                                                              <  2:695086/140‑1 (MQ=255)                    aCTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGaa                                                                                                             <  2:219580/140‑1 (MQ=255)                       gctgAATTCACTGTCGGCGAAGGAAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAacac                                                                                                          >  2:788984/1‑140 (MQ=255)                         tgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACAccc                                                                                                        <  2:762110/140‑1 (MQ=255)                           aaTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTa                                                                                                      >  1:113793/1‑140 (MQ=255)                                aCTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAccccc                                                                                                 <  1:160113/140‑1 (MQ=255)                                   gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaaa                                                                                              <  1:299/140‑1 (MQ=255)                                   gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaaa                                                                                              <  2:635190/140‑1 (MQ=255)                                     cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaaa                                                                                              <  1:402954/138‑1 (MQ=255)                                     cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaaa                                                                                              >  2:402954/1‑138 (MQ=255)                                        cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATtata                                                                                         <  1:389850/140‑1 (MQ=255)                                            ggTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATaaa                                                                                     >  1:439529/1‑140 (MQ=255)                                             gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAg                                                                                    <  1:42857/140‑1 (MQ=255)                                              tAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGa                                                                                   <  2:23781/140‑1 (MQ=255)                                              tAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACATTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaaa                                                                                              <  1:711442/129‑1 (MQ=255)                                              tAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACATTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaaa                                                                                              >  2:711442/1‑129 (MQ=255)                                                  ttGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGa                                                                              <  2:397998/140‑1 (MQ=255)                                                   tgctgaCTCATGATCAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGAc                                                                             <  2:509875/137‑1 (MQ=255)                                                    gatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACa                                                                            <  2:113793/140‑1 (MQ=255)                                                      tgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGAc                                                                             <  1:80050/137‑1 (MQ=255)                                                      tgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGAc                                                                             >  2:80050/1‑137 (MQ=255)                                                       gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAata                                                                         <  2:439529/140‑1 (MQ=255)                                                           cATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCtt                                                                     >  1:281081/1‑140 (MQ=255)                                                             tgatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGa                                                                   <  1:472212/140‑1 (MQ=255)                                                                    cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGacaaacaa                                                            >  1:89366/1‑140 (MQ=255)                                                                              tGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGtt                                                  >  2:134787/1‑140 (MQ=255)                                                                                 cTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTAt                                               >  2:26091/1‑140 (MQ=255)                                                                                  tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTAtt                                              <  2:89366/140‑1 (MQ=255)                                                                                   ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGt                                                   <  1:653126/134‑1 (MQ=255)                                                                                   ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGt                                                   >  2:653126/1‑134 (MQ=255)                                                                                   ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTAttt                                             >  2:436806/1‑140 (MQ=255)                                                                                       aTGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCaa                                         <  1:26091/140‑1 (MQ=255)                                                                                            aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCa                                          <  2:506187/134‑1 (MQ=255)                                                                                            aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCa                                          >  1:506187/1‑134 (MQ=255)                                                                                                          atCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTCCTAAAAGAACACCC‑ACCCCCAAATTATTAAAAGGTGGGACATTAT‑TTGAACAACAAATCCAAAGGTTAATTTTAAAAATATTAACCATTAAttaa                     >  2:430459/1‑137 (MQ=255)                                                                                                          atCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGa                      >  2:341818/1‑140 (MQ=255)                                                                                                          atCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGa                      >  1:682801/1‑140 (MQ=255)                                                                                                           tCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGAc                     <  2:452067/140‑1 (MQ=255)                                                                                                              gggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCAt                            <  2:117268/130‑1 (MQ=255)                                                                                                              gggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCAt                            >  1:117268/1‑130 (MQ=255)                                                                                                              gggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAAc                  <  2:300190/140‑1 (MQ=255)                                                                                                               ggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAAcc                 >  2:201730/1‑140 (MQ=255)                                                                                                               ggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAAcc                 >  2:797904/1‑140 (MQ=255)                                                                                                                  cTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCttt              >  2:608194/1‑140 (MQ=255)                                                                                                                  cTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCttt              >  1:117657/1‑140 (MQ=255)                                                                                                                        cGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAg        <  2:117657/140‑1 (MQ=255)                                                                                                                              ttccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAg  <  1:134787/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                          TTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAT‑GGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAG  >  NC_000913/3653102‑3653366

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.