New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1397238	NA (NA)	11 (0.340)	9/276	1.7	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	= 1978494	0 (0.000)					intergenic (‑297/+24)	flhD/insB‑5	DNA‑binding transcriptional dual regulator FlhD/IS1 protein InsB

GCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGA  >  NC_000913/1397107‑1397377                                                                                                                                    |                                                                                                                                            gCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTTAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCctaactaa                                                                                                                                     >  2:1663397/1‑140 (MQ=17)     tAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCaa                                                                                                                                 >  1:1308186/1‑140 (MQ=34)         tGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCg                                                                                                                             <  2:1640101/140‑1 (MQ=39)         tGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCg                                                                                                                             <  2:1279041/140‑1 (MQ=39)          ggCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAACCATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGt                                                                                                                            <  2:1264159/140‑1 (MQ=255)                 aaaCGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGt                                                                                                                     >  2:1769019/1‑140 (MQ=255)                      gTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAGGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATc                                                                                                                <  2:1511692/140‑1 (MQ=255)                                 cTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAgagg                                                                                                     >  1:1670369/1‑140 (MQ=255)                                     ttCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAgaggtgag                                                                                                 <  1:1347220/140‑1 (MQ=255)                                     ttCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAgaggtgag                                                                                                 >  2:1428684/1‑140 (MQ=255)                                          ggCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGtttt                                                                                            <  2:1565613/140‑1 (MQ=255)                                                   gggAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGa                                                                                   <  2:1686763/140‑1 (MQ=255)                                                     gagaAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACa                                                                                 <  1:1783043/140‑1 (MQ=255)                                                      agaAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACaa                                                                                <  1:1769019/140‑1 (MQ=255)                                                         aaaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACgg                                                                             <  1:1601559/140‑1 (MQ=255)                                                             aTCGGCTTAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGaa                                                                         <  1:1663397/140‑1 (MQ=255)                                                                ggCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATg                                                                      <  2:1615694/140‑1 (MQ=255)                                                                   tCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCa                                                                   <  1:1675287/140‑1 (MQ=255)                                                                      aaCATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATgg                                                                <  2:1272355/140‑1 (MQ=255)                                                                            aagaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGt                                                          <  1:1428684/140‑1 (MQ=255)                                                                                            gAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACaa                                          >  2:1568236/1‑140 (MQ=255)                                                                                               tCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGaa                                       >  1:1673329/1‑140 (MQ=255)                                                                                                     gagaAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGt                                 <  2:1670369/140‑1 (MQ=255)                                                                                                       gaagaaTTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCt                               <  1:1578214/140‑1 (MQ=255)                                                                                                                 cccGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGAt                     <  2:1348223/140‑1 (MQ=255)                                                                                                                     cTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTAca                 >  1:1700984/1‑140 (MQ=255)                                                                                                                                   cccTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCg   <  1:1403996/140‑1 (MQ=255)                                                                                                                                    ccTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGa  >  1:1334085/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                            GCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGA  >  NC_000913/1397107‑1397377

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