New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 148496737 (0.930)5 (0.130) 5/246 4.3 13.2% coding (1907/3903 nt) hrpA putative ATP‑dependent helicase
?NC_000913 = 1484972 31 (0.840)coding (1912/3903 nt) hrpA putative ATP‑dependent helicase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GCGCGTGCGCGAATGGCAGGATATCTACACCCAGTTGCGTCAGGTGGTGAAAGAACTTGGCATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGG  >  NC_000913/1484827‑1485100
                                                                                                                                                 |                                                                                                                                
gcgcgtgcgcgAATGGCAGGATATCTACACCCAGTTGCGTCAGGTGGTGAAAGAACTTGGCATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGAt                                                                                                                                <  4:91097/148‑1 (MQ=255)
 cgcgtgcgcgAATGGCAGGATATCTACACCCAGTTGCGTCAGGTGGTGAAAGAACTTGGCATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAg                                                                                                                                  <  6:123353/145‑1 (MQ=255)
 cgcgtgcgcgAATGGCAGGATATCTACACCCAGTTGCGTCAGGTGGTGAAAGAACTTGGCATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGAt                                                                                                                                <  6:154284/147‑1 (MQ=255)
 cgcgtgcgcgAATGGCAGGATATCTACACCCAGTTGCGTCAGGTGGTGAAAGAACTTGGCATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGc                                                                                                                              <  6:177056/149‑1 (MQ=255)
               gCAGGATATCTACACCCAGTTGCGTCAGGTGGTGAAAGAACTTGGCATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGaa                                                                                                                 >  4:23014/1‑148 (MQ=255)
                cAGGATATCTACACCCAGTTGCGTCAGGTGGTGAAAGAACTTGGCATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAata                                                                                                               <  3:23014/149‑1 (MQ=255)
                                              gtgAAAGAACTTGGCATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAAcgc                                                                                             <  8:76704/137‑1 (MQ=255)
                                              gtgAAAGAACTTGGCATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAAcgc                                                                                             >  7:76704/1‑137 (MQ=255)
                                                gAAAGAACTTGGCATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCg                                                                                                          >  6:34883/1‑122 (MQ=255)
                                                gAAAGAACTTGGCATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCg                                                                                                          <  5:34883/122‑1 (MQ=255)
                                                            cATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTg                                                                   <  6:147596/149‑1 (MQ=255)
                                                             aTTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTgg                                                                  >  5:65826/1‑149 (MQ=255)
                                                             aTTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTATCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAAc                                                                                               >  6:28496/1‑120 (MQ=255)
                                                                  ggTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTAt                                                             >  8:51178/1‑149 (MQ=255)
                                                                       aCAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGTATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCaa                                                         >  4:100993/1‑148 (MQ=255)
                                                                          gCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCaaaaaa                                                     >  8:87793/1‑149 (MQ=255)
                                                                                  gcggAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTa                                              >  2:181062/1‑148 (MQ=255)
                                                                                  gcggAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTa                                              >  4:32642/1‑148 (MQ=255)
                                                                                   cggAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCaaaaaa                                                     >  4:17020/1‑140 (MQ=255)
                                                                                   cggAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCaaaaaa                                                     <  3:17020/140‑1 (MQ=255)
                                                                                    ggAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAAt                                           >  5:88985/1‑149 (MQ=255)
                                                                                       gTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATggg                                        <  1:181062/149‑1 (MQ=255)
                                                                                            gcgAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAAtg                                   <  7:87793/149‑1 (MQ=255)
                                                                                                 aTTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGTATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGc                              <  3:100993/149‑1 (MQ=255)
                                                                                                        ttGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTg                       <  6:65826/149‑1 (MQ=255)
                                                                                                        ttGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTg                       <  6:88985/149‑1 (MQ=255)
                                                                                                           cGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGt                     <  7:51178/148‑1 (MQ=255)
                                                                                                                cTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCt                                               >  4:114580/1‑117 (MQ=255)
                                                                                                                cTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCt                                               <  3:114580/117‑1 (MQ=255)
                                                                                                                 tGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTa                    >  2:143441/1‑143 (MQ=255)
                                                                                                                 tGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTa                    <  1:143441/143‑1 (MQ=255)
                                                                                                                           aCTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCtgt    >  3:64556/1‑149 (MQ=255)
                                                                                                                             tttCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAAt                                           >  2:101115/1‑108 (MQ=255)
                                                                                                                             tttCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAAt                                           <  1:101115/108‑1 (MQ=255)
                                                                                                                             tttCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTgg  >  1:165959/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                
GCGCGTGCGCGAATGGCAGGATATCTACACCCAGTTGCGTCAGGTGGTGAAAGAACTTGGCATTCCGGTTAACAGCGAACCGGCGGAGTATCGCGAAATTCACATTGCGTTGCTGACCGGTTTACTTTCCCATATCGGCATGAAAGATGCCGATAAACAAGAATATACCGGCGCACGTAACGCGCGTTTCTCCATCTTCCCCGGTTCTGGTTTATTCAAAAAACCGCCTAAATGGGTAATGGTGGCGGAACTGGTAGAAACCAGCCGCCTGTGG  >  NC_000913/1484827‑1485100

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