New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_pgi	2514277 =	NA (NA)	65 (0.870)	56/290	0.3	91.7%	intergenic (+33/‑54)	nupC/insL1	nucleoside (except guanosine) transporter/IS186 transposase
	?	NC_000913_3_hsa_pgi	4543633 =	6 (0.080)					coding (491/549 nt)	fimA	major type 1 subunit fimbrin (pilin)

AAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGC  >  NC_000913_3_hsa_pgi/2514131‑2514419                                                                                                                                                   |                                                                                                                                               aaaGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATa                                                                                                                                              >  2:477293/1‑149 (MQ=255) aaaGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAACACCATACATAAAAAAGCCGGGGATAATTCCCATa                                                                                                                                              >  1:324854/1‑149 (MQ=255) aaaGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCCGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATa                                                                                                                                              >  1:144991/1‑149 (MQ=255)        tgaatgaaGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGAGAGTGTTCAGTCAGCGTCAATCGCAGCACTGGTGCTGGAAGCCCTAACATAAAAAAGCCGGGGATAATTACCAGAAGAGCTa                                                                                                                                       >  2:737856/1‑149 (MQ=255)           atgaaGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACt                                                                                                                                    >  1:930953/1‑149 (MQ=255)                   cAAGGTAACGTGGTTTCTCGCTTCGTGCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTTCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGtt                                                                                                                            <  2:139367/149‑1 (MQ=255)                   cAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGtt                                                                                                                            <  1:477293/149‑1 (MQ=255)                     aGGTAACGTGGTTTCTCGCTTCGGTCTGAATCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAATGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTtgt                                                                                                                          <  1:826167/149‑1 (MQ=255)                     aGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTtgt                                                                                                                          >  1:896507/1‑149 (MQ=255)                        tAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGt                                                                                                                       >  1:635181/1‑149 (MQ=255)                         aaCGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGTCTATAGCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGATGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                                                                                                      <  1:1051937/149‑1 (MQ=255)                          aCGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCGCATAAGCGCTAACTTAAGGCATGTGGTAt                                                                                                                     >  1:614216/1‑149 (MQ=255)                           cGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                                                                                                    <  1:519732/149‑1 (MQ=255)                               gTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                                                                                >  2:1141099/1‑149 (MQ=255)                               gTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                                                                                                <  2:1017951/149‑1 (MQ=255)                                             tCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                                                                                                  >  1:371840/1‑149 (MQ=255)                                               tGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCACATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                                                                >  1:70461/1‑149 (MQ=255)                                                   gCTGGTTTACGGCTCTACCCTGGTGTGTGTGCTTTCTGCGTCAATAGCAGAAATGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAAATCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc                                                                                            <  1:581119/149‑1 (MQ=37)                                                     tGGTTTACGGCGCTACCCTGGTGAGTGTGCTGGCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa                                                                                          <  2:1171922/149‑1 (MQ=255)                                                      ggTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                         >  1:1077092/1‑149 (MQ=255)                                                       gTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg                                                                                        <  2:930953/149‑1 (MQ=255)                                                         ttACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa                                                                                      <  1:1141099/149‑1 (MQ=255)                                                            cGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                                                                   <  2:1112154/149‑1 (MQ=255)                                                              gCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                                                 <  2:201749/149‑1 (MQ=255)                                                                 ctACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc                                                                              >  1:889769/1‑149 (MQ=255)                                                                 ctACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc                                                                              <  2:896507/149‑1 (MQ=255)                                                                   aCCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc                                                                            <  2:1027128/149‑1 (MQ=255)                                                                                  cTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCa                                                             >  1:1236851/1‑149 (MQ=255)                                                                                   tGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCaa                                                            <  2:371840/149‑1 (MQ=255)                                                                                    gTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAt                                                           >  1:543320/1‑149 (MQ=255)                                                                                      cTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGGTGTGGTATTACGCGTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTc                                                         >  2:471963/1‑149 (MQ=255)                                                                                        gCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTg                                                       >  2:1164159/1‑149 (MQ=255)                                                                                          ggcAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTATCGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                                     <  1:95950/147‑1 (MQ=255)                                                                                          gTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                                     <  2:635181/149‑1 (MQ=255)                                                                                          gTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTACCATAAGAGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgtc                                                     >  1:1120483/1‑147 (MQ=255)                                                                                                   agcaCTGGTGCTGTTAGAGCTAACATAAAAATGCCGGGGATAATTCCCATAAGCGCAAACTTAATTGTTGTGGTATTACGCCTTTTATGTTTTAAGGTGCAGATGAATTACTCTCACGATCACTGGTCAGCGATTCTGGCCCATATTgg                                            <  1:255656/149‑1 (MQ=255)                                                                                                        ttgTGATGTAAGACAATATATAAAAAAGTGGGGGATAATTCCCATAAGCTATAACTTAAGGGTTGTGGTATTACGATTTATATGATTTAACTTGCCGTTGAATTAGTCTCACTATAACTGTTCAGCAATTCTGTCCCATATTGGTAAgg                                       <  2:614216/147‑2 (MQ=14)                                                                                                         agTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGcc                                      >  1:1136171/2‑149 (MQ=255)                                                                                                          gTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTTAGGGTTGTGGTATTACGCCTGATATTATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCACATATTGGTAAGcca                                     >  2:1166985/1‑148 (MQ=255)                                                                                                             cTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCAGATGATTTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTATGCCCgaa                                  >  2:890064/1‑149 (MQ=255)                                                                                                              tGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaag                                 >  2:1198308/1‑149 (MQ=255)                                                                                                               gTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaaga                                >  1:1056805/1‑149 (MQ=255)                                                                                                                  aGACCATACAAAAAAAAGCCGGCGATAATTCAGATAAGCGCTAACTTAAGGGTTGTGGAATTACGCCTGATATGATTTAACGTGCCGATGAATTACTATCACGATAACAGGTCAGCAATTCTTGCCCATATTGGTAAGCCCAAAGAAca                             >  1:771231/1‑148 (MQ=14)                                                                                                                    aCCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTgg                           <  1:369432/149‑1 (MQ=255)                                                                                                                      catacataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCTCGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGAt                         >  2:279590/1‑149 (MQ=255)                                                                                                                        tacataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGCACTGGATAc                       >  1:266701/1‑149 (MQ=255)                                                                                                                           ataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTc                    >  2:1092962/1‑149 (MQ=255)                                                                                                                             aaaaaaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCgg                  >  2:477319/1‑149 (MQ=40)                                                                                                                               aaaaTCCGGGGATAATTCTCATTATGGCTAACTTACGGTTTGTGGTATTACGCCTGATATGATTTAACTTGCCTATGAATTACTCTATCGTTAACTGGTCAGCAATTCTGGCCCATATTGGTACGCCCGAAGAACTGGATACTTCGGCa                <  1:737856/149‑1 (MQ=11)                                                                                                                                 aaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACg              >  1:613183/1‑149 (MQ=38)                                                                                                                                    ccGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTaa           >  2:366044/1‑149 (MQ=35)                                                                                                                                     cGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAAt          >  1:543696/1‑149 (MQ=34)                                                                                                                                       gggATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGc        >  1:1209062/1‑149 (MQ=34)                                                                                                                                          aTAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATTCATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCAGAATAACTGGATACTTCGGCACGTAATGCCgg     >  2:782320/1‑149 (MQ=255)                                                                                                                                             aTTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGc  >  1:881339/1‑149 (MQ=31)                                                                                                                                                   |                                                                                                                                               AAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGC  >  NC_000913_3_hsa_pgi/2514131‑2514419

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