New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4168568121 (1.050)45 (0.420) 17/184 9.9 28.6% intergenic (+121/‑73) gltT/rrlB tRNA‑Glu/23S ribosomal RNA of rrnB operon
?NC_000913 = 4209974 NA (NA)intergenic (+125/‑69) gltV/rrlE tRNA‑Glu/23S ribosomal RNA of rrnE operon
Rejected: Coverage evenness skew score above cutoff.

TTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGC  >  NC_000913/4209921‑4210073
                                                     |                                                                                                   
tttGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGt                                                       >  1:221333/1‑100 (MQ=11)
tttGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGt                                                       <  1:2049240/100‑1 (MQ=11)
 ttGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTa                                                      >  2:2081373/1‑100 (MQ=31)
 ntGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTa                                                      >  2:2784354/2‑100 (MQ=14)
       tttAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAAc                                                <  2:1090240/100‑1 (MQ=33)
        ttAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACa                                               <  2:1079558/100‑1 (MQ=33)
        ttAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACa                                               <  2:2711325/100‑1 (MQ=33)
        ttAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACa                                               <  1:87297/100‑1 (MQ=33)
           aaaaTCTGGAACAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTATCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACAt                                              >  2:1170641/1‑98 (MQ=12)
              aTCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCg                                         <  1:185845/100‑1 (MQ=33)
              aTCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCg                                         >  2:108320/1‑100 (MQ=33)
               tCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAGATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCgg                                        <  2:1586119/100‑1 (MQ=18)
               tCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCgg                                        <  2:39375/100‑1 (MQ=33)
               tCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCgg                                        <  2:1651186/100‑1 (MQ=33)
                cTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCggg                                       <  1:2651675/100‑1 (MQ=33)
                 tGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGt                                      <  2:904081/100‑1 (MQ=33)
                 tGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGt                                      <  2:214527/100‑1 (MQ=33)
                  ggATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGtt                                     >  2:2139131/1‑100 (MQ=33)
                   gATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTtg                                    <  1:2417927/100‑1 (MQ=33)
                   gATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTtg                                    <  2:944458/100‑1 (MQ=33)
                    aTCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTtgt                                   >  1:275252/1‑100 (MQ=33)
                    aTCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTtgt                                   <  2:1517352/100‑1 (MQ=33)
                    aTCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTtgt                                   >  1:1839948/1‑100 (MQ=33)
                       aaGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAg                                >  1:71350/1‑100 (MQ=33)
                       aaGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAg                                <  1:1761754/100‑1 (MQ=33)
                       aaGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAg                                >  2:610889/1‑100 (MQ=32)
                       aaGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAg                                <  2:871299/100‑1 (MQ=33)
                        aGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAgg                               <  1:314042/100‑1 (MQ=33)
                         gCTGAAAATTGCAACACTGACCAACGAAAGTTGTTCGTGAGTCTCTCAACTTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGt                              <  2:958646/100‑1 (MQ=255)
                         gCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGt                              >  2:780985/1‑100 (MQ=33)
                          cTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGtt                             <  1:815960/100‑1 (MQ=33)
                           tGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTa                            >  2:88557/1‑100 (MQ=33)
                           tGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTa                            >  2:866728/1‑100 (MQ=33)
                           tGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCGCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTGGGGTTGTGAGGTTa                            >  1:918251/1‑100 (MQ=16)
                            gAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTaa                           >  1:896710/1‑100 (MQ=33)
                            gAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTaa                           >  1:2225466/1‑100 (MQ=33)
                             aaaaTTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAg                          <  1:1648622/100‑1 (MQ=33)
                             aaaaTTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAg                          <  2:2761479/100‑1 (MQ=33)
                              aaaTTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGc                         <  1:2118440/100‑1 (MQ=33)
                               aaTTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCg                        >  1:155931/1‑100 (MQ=33)
                                aTTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGa                       >  1:379044/1‑100 (MQ=33)
                                  tGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACt                     >  1:1307326/1‑100 (MQ=33)
                                  tGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACt                     <  1:88557/100‑1 (MQ=33)
                                   gAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTa                    <  1:1908430/100‑1 (MQ=33)
                                   gAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTa                    >  2:271706/1‑100 (MQ=33)
                                     aaCACTGAACCACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGACTCGTAAGAATCATCTTCGGGTTGTGAGGTTAAGCGACTAAg                  <  2:1416952/100‑1 (MQ=255)
                                      ctcacgAACAACGAAAGTTGTTCGTGAGTCTCTCGAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGc                 >  1:2451140/6‑100 (MQ=255)
                                       cacTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                >  2:692339/1‑100 (MQ=33)
                                        acTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGt               >  2:2832603/1‑100 (MQ=33)
                                        acTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGt               <  1:548378/100‑1 (MQ=33)
                                         cTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg                >  2:1754989/1‑98 (MQ=32)
                                          tGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTac             <  1:2392392/100‑1 (MQ=33)
                                          tGAACAACGAAAGATGTTCGTGAGTATCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTac             >  2:1786760/1‑100 (MQ=12)
                                              caacGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg         >  1:327307/1‑100 (MQ=33)
                                               aacGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGt        >  2:245033/1‑100 (MQ=32)
                                               aacGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGt        <  1:868914/100‑1 (MQ=33)
                                               aacGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGt        >  1:1754527/1‑100 (MQ=33)
                                                acGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg       >  1:2583808/1‑100 (MQ=32)
                                                acGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg       >  2:1149269/1‑100 (MQ=33)
                                                acGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg       <  1:879526/100‑1 (MQ=33)
                                                acGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg       >  2:2759676/1‑100 (MQ=33)
                                                 cGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTgg      <  2:243891/100‑1 (MQ=33)
                                                 cGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTgg      <  2:243890/100‑1 (MQ=33)
                                                 cGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTgg      <  2:2152713/100‑1 (MQ=32)
                                                 cGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTgg      <  2:2152493/100‑1 (MQ=33)
                                                 cGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTgg      <  1:953718/100‑1 (MQ=32)
                                                     aGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGc  >  1:1562649/1‑100 (MQ=33)
                                                     |                                                                                                   
TTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGTAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGC  >  NC_000913/4209921‑4210073

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