New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3426836 =NA (NA)12 (0.200) 10/184 6.8 NA intergenic (‑53/+122) rrlD/alaU 23S ribosomal RNA of rrnD operon/tRNA‑Ala
?NC_000913 3426854 = NA (NA)intergenic (‑71/+104) rrlD/alaU 23S ribosomal RNA of rrnD operon/tRNA‑Ala
Rejected: Coverage evenness skew score above cutoff.

CAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAAC  >  NC_000913/3426749‑3426910
                                                                                       |                                                                          
cAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCgt                                                                >  1:1000888/1‑100 (MQ=31)
  gggCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCgttg                                                              >  2:246646/1‑100 (MQ=33)
  gggCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCgttg                                                              <  2:1349834/100‑1 (MQ=32)
  gagCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCgttg                                                              >  1:1008242/3‑100 (MQ=17)
     cATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTc                                                           <  2:1110913/100‑1 (MQ=33)
        ccACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAgt                                                        >  2:1585833/1‑100 (MQ=32)
         cACCGTGTACGCTTAGTCGCTTAACCTCCCAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAgtg                                                       <  2:1331831/100‑1 (MQ=18)
         cACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAgtg                                                       <  1:1104153/100‑1 (MQ=33)
         cACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAgtg                                                       <  2:539219/100‑1 (MQ=33)
           ccGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGtt                                                     >  2:56657/1‑100 (MQ=32)
            cGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGttt                                                    >  1:618907/1‑100 (MQ=33)
             gtgtACGCTTAGTCGCTTAACCTCACAACCCGACGATGTGTCACTTCGGTGTTGCGAACATTTGAGAGACTCACGACCATCTTTCGTTGTTCAGTGttt                                                    >  1:760222/1‑99 (MQ=11)
             gtgtACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTc                                                   >  1:402227/1‑100 (MQ=32)
                 aCGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAACATTTGAGAGACTCACTACCAACTTTCGTTGTTCAGCGCTTCAAtt                                               >  1:674446/1‑100 (MQ=11)
                  cGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAAttt                                              >  2:581895/1‑100 (MQ=33)
                    cTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGTACAACTTTCGTTGTTCACTGTTTCAATTTTc                                            <  2:1075031/100‑1 (MQ=16)
                    cTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTc                                            <  2:341659/100‑1 (MQ=33)
                      tAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAg                                          >  1:753211/1‑100 (MQ=32)
                       aGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGc                                         >  2:1501407/1‑100 (MQ=33)
                        gTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCt                                        >  2:274006/1‑100 (MQ=33)
                          cGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTg                                      >  2:462300/1‑100 (MQ=33)
                          cGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTg                                      >  2:1073142/1‑100 (MQ=33)
                           gCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGa                                     >  2:1555550/1‑100 (MQ=33)
                            cTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGAt                                    >  1:1522498/1‑100 (MQ=33)
                             ttAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATc                                   >  1:75525/1‑100 (MQ=32)
                                aCCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAg                                <  2:1008252/100‑1 (MQ=33)
                                aCCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAg                                >  2:454991/1‑100 (MQ=33)
                                  cTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGAt                              >  2:387241/1‑100 (MQ=32)
                                   tCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGAtt                             >  1:70506/1‑100 (MQ=33)
                                     acaACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGAtttt                           <  1:1169510/100‑1 (MQ=33)
                                       gatctgAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTa                         <  2:772854/95‑1 (MQ=18)
                                         cccGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTaaa                       <  1:697895/100‑1 (MQ=32)
                                          ccGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAag                      <  1:1525578/100‑1 (MQ=33)
                                            gAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTaagatc                    >  1:772847/1‑96 (MQ=18)
                                            gAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAagag                    <  2:860355/100‑1 (MQ=33)
                                            gAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAagag                    <  1:28563/100‑1 (MQ=33)
                                             aaGATTTTTCACTTCTGAGTTGGGAAAATTTGTTTTTCTCACGAACCACTTCCGTTGTTCAGTGTTTCACTTTTCACCTTGATCCCGATTTTTAAaga                     <  1:273766/98‑1 (MQ=2)
                                                  gTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAAc              <  2:1261800/100‑1 (MQ=35)
                                                    ttCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACtt            >  1:786425/1‑100 (MQ=38)
                                                      cACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCg          >  1:811798/1‑100 (MQ=40)
                                                              agTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAAc  <  1:1006240/100‑1 (MQ=11)
                                                              agTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAAc  <  1:1006190/100‑1 (MQ=11)
                                                                                       |                                                                          
CAGGGCATCCACCGTGTACGCTTAGTCGCTTAACCTCACAACCCGAAGATGTTTCACTTCAGAGTTGCGAAAATTTGAGAGACTCACGAACAACTTTCGTTGTTCAGTGTTTCAATTTTCAGCTTGATCCAGATTTTTAAAGAGCAAAACTTCGCAGTGAAC  >  NC_000913/3426749‑3426910

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