New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3944977NA (NA)11 (0.090) 7/284 16.2 17.8% noncoding (1275/2904 nt) rrlC 23S ribosomal RNA of rrnC operon
?NC_000913 = 4038810 53 (0.430)noncoding (1293/2905 nt) rrlA 23S ribosomal RNA of rrnA operon
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAGGCGA  >  NC_000913/4038664‑4038875
                                                                                                                                                   |                                                                 
aTGCACCGAAGCTGCGGCAGCGGCACTATGTGTTGTTGGGTCAGGGCGCGTGCTGTACGCCT‑GTGAGGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGAGAATGCTGACATCAGTAACGATAAAGCGGGTGAAAAGCGCGCTCg                                                                 >  1:1500252/1‑149 (MQ=14)
 tGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGc                                                                >  2:826605/1‑149 (MQ=255)
 tGCACCGAAGCTGCGGCAGCGACACTATGTGTTGATGGTTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGc                                                                >  2:1964023/1‑149 (MQ=255)
 tGCACCGAAGCCGCGGCAGCGACACGATGTGTTGGTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGc                                                                <  1:1577661/149‑1 (MQ=255)
  gCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGcc                                                               <  1:1312313/149‑1 (MQ=255)
  gCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGcc                                                               <  1:203166/149‑1 (MQ=40)
  gCACCGAAGCTGCGGCAGCGACACTATGCGGTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTCACATTGGCAACGGTAAAGCGGGTGAAAAGCCCGCTCGcc                                                               <  2:1967168/149‑1 (MQ=11)
    aCCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGAGAAGGCTGACATAAGTAACGGTAAAACGGGCGAAAAGcccgcccgtggg                                                             >  2:1400944/1‑145 (MQ=14)
    aCCGAAGCTGCGGCAGCGACACTATGTGGTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAAGGCTGACATAAGTAAAGATAAAGCGGGTGAAAAGCCCGCTCGCCgt                                                             >  2:928661/1‑148 (MQ=255)
      ccaaGCGGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                           <  1:546089/147‑1 (MQ=255)
      cGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGaa                                                           <  2:107634/149‑1 (MQ=255)
       gAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAg                                                          >  1:738244/1‑149 (MQ=255)
       gAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAAGGGAGCGTTCTGTAAGCCT‑GGGAAGGTGTGCTGTGAGGCATGCTGGAGGAATCAGAAGTGCGAATCCTGACAGAAGTAACGATAAACCGGGGGAAAAGGCCCCACGCAGGAAg                                                          >  2:1235390/1‑149 (MQ=11)
        aggCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCTTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGACgcagta                                                         <  2:1773960/147‑6 (MQ=255)
        aaGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGCAGGTGTGCTGTGAGGCATGCTGGAGGGATCAGAAGTGCGAAGGCTGACATAAGTAACGAAAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGa                                                         >  2:127073/1‑149 (MQ=255)
              cGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAgg                                                   >  2:1664572/1‑149 (MQ=255)
              cGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAgg                                                   <  2:693467/149‑1 (MQ=40)
              cGGCAGCGACACTATGTGTTGTTGGGCAGGGGAGCGTTCTGTAAGCCT‑GTGACGGTGTGATGTGAGGCATGCTGGAGGTAGCAGAAGTGCGAAAGCTGACACAAGTCAAGATAAAGCGGGGAAAACGCCCGCGCGGCGGCAGACCACgg                                                   >  2:1083916/1‑149 (MQ=14)
               ggCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGACGCCCaaccg                                                  >  2:1913063/1‑146 (MQ=255)
                  aGCGACACTGTGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTc                                               <  2:1281427/149‑1 (MQ=255)
                  aGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGATGTGCGAATGCTGACATAAGTAACGATAAAGTGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTc                                               <  1:1602870/149‑1 (MQ=255)
                  aGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTc                                               <  2:1855525/149‑1 (MQ=40)
                  aGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTc                                               <  1:234065/149‑1 (MQ=40)
                     gACACTCTGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTGACGATAAAGCGGGAGAAAAGCCCGCTCGCCGGAAGGCCAGGGGTTCCTg                                            >  2:1123553/1‑149 (MQ=11)
                      acacTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGt                                           >  2:1994727/1‑149 (MQ=38)
                       cacTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTc                                          <  2:546692/149‑1 (MQ=38)
                        acTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTcc                                         >  2:1717327/1‑149 (MQ=35)
                          tATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCaa                                       >  1:710514/1‑149 (MQ=35)
                          tATGTGTTGTTGGGTAGGGGAGCGTTCTGGAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCaa                                       <  1:1658896/149‑1 (MQ=25)
                          tATGTGTTGTTGGGTAGGGGAGCGTGCTGTAAGCCT‑GTGGAGGTGTGCTGTGAGGCATGCTGGGGGTATCGGGAGTGCGAGTGCTGACGTAGGTAACGATAAAGCGGGTGGAAAGCCCGGTCGCCGGAAGCCCAAGGGTTGCTGTCCag                                       >  2:571794/1‑148 (MQ=11)
                            ggtgtTGGTGGGGATGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCCGGGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTAGCCGGAAGACCATGGGTTCCTGTCCAACg                                     <  1:1283671/148‑1 (MQ=11)
                              tgttgtTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGtt                                   >  1:1966299/1‑149 (MQ=33)
                               gttgttGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGACCGCTCGCCGGAAGACCAAGGGTTCCCGTCCCACGTTa                                  >  2:395430/1‑149 (MQ=12)
                                   tggggTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCg                              <  2:1898733/147‑1 (MQ=17)
                                     gggTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCggg                            >  1:1643643/1‑149 (MQ=32)
                                      ggTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCGGCTCGCCGGAAGCCCAAGGGTTCCTGCCCAAAGTTAATCgggc                           >  2:505100/1‑148 (MQ=255)
                                          agggAGAGTTCTGGCATCCTGGTGAA‑GTGTGCTGCGCGGCATGCTTGAGGTATCAGAAGTGCGAATGCGGGCATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAATACCAAGGGTTCCTGTCCAACGTTAATTGGGGCAgg                       <  2:653996/148‑1 (MQ=11)
                                           gggAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAggg                      >  2:133753/1‑149 (MQ=32)
                                            ggAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGt                     >  2:703412/1‑149 (MQ=32)
                                              aGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTTACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGa                   <  1:1031556/149‑1 (MQ=17)
                                                  ttCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCg               >  2:1364628/1‑149 (MQ=32)
                                                     tGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGTGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAcc            >  2:267351/1‑149 (MQ=17)
                                                     tGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAcc            >  1:1155601/1‑149 (MQ=32)
                                                     tGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAcc            >  1:735961/1‑149 (MQ=32)
                                                      gtgaGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAccc           <  2:1190716/146‑1 (MQ=17)
                                                      gTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCCGAAGCGCGAACGCTGACATAAGTAACGATAAAGCGGGGGGAGAGCCCGCTCGCGGGAAGACCCAGGCCTGATGTCCAACGTTAATCGGGGACGGATGCGTCGAccc           >  2:784940/1‑149 (MQ=11)
                                                       tAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGAcccc          >  2:748290/1‑149 (MQ=32)
                                                          gCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTaa       >  1:1567818/1‑149 (MQ=32)
                                                          gCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTaa       <  1:1053368/149‑1 (MQ=32)
                                                          gCCT‑GTGAAGGTGGGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTaa       <  2:87218/149‑1 (MQ=17)
                                                            cT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGAGGCAGGGTGAGGCGACCCCTAAgg     >  1:840391/1‑149 (MQ=17)
                                                               gtgAAGGTGTGCTGTGAGGAATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGCTAAAGCGGGTGAAAAGCCCGCACGCGGGAAGACCAGGGCTTCCTGTCCAACGTTAATCGGGGCAGGGCGAGTCGACCCCGAAGGca   >  2:789894/1‑148 (MQ=11)
                                                                tgAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAGGCGa  >  1:600145/1‑149 (MQ=11)
                                                                                                                                                   |                                                                 
ATGCACCGAAGCTGCGGCAGCGACACTATGTGTTGTTGGGTAGGGGAGCGTTCTGTAAGCCT‑GTGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACATAAGTAACGATAAAGCGGGTGAAAAGCCCGCTCGCCGGAAGACCAAGGGTTCCTGTCCAACGTTAATCGGGGCAGGGTGAGTCGACCCCTAAGGCGA  >  NC_000913/4038664‑4038875

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