New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4168165 =52 (0.610)12 (0.150) 5/168 NT 17.5% noncoding (1507/1542 nt) rrsB 16S ribosomal RNA of rrnB operon
?NC_000913 4168187 = 64 (0.790)noncoding (1529/1542 nt) rrsB 16S ribosomal RNA of rrnB operon

ACACCATGGGAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCG  >  NC_000913/4168066‑4168264
                                                                                                   |                                                                                                   
gcctcactcccaacgcgtaaacgccttgcttttcactttctatcagacaatctgtgtgagcactacaaagtacgcttctttaaggtaaggaggtgatcca                                                                                                     <  1:1625261‑M1/1‑1 (MQ=255)
gcctcactcccaacgcgtaaacgccttgcttttcactttctatcagacaatctgtgtgagcactacaaagtacgcttctttaaggtaaggaggaaatcca                                                                                                     <  1:2278780‑M1/1‑1 (MQ=255)
 actcactcccaacgcgtaaacgccttgcttttcactttctatcagacaatctgtgtgagcactacaaagtacgcttctttaaggtaaggaggtgatccaa                                                                                                    <  1:1644057‑M1/2‑1 (MQ=255)
  ctcactcccaacgcgtaaacgccttgcttttcactttctatcagacaatctgtgtgagcactacaaagtacgcttctttaaggtaaggaggtgatccAAc                                                                                                   <  1:1814849‑M1/3‑1 (MQ=255)
  cgcactcccaacgcgtaaacgccttgcttttcactttctatcagacaatctgtgtgagcactacaaagtacgcttctttaaggtaaggaggtgatccAAc                                                                                                   <  2:1035876‑M1/3‑1 (MQ=255)
    cactcccaacgcgtaaacgccttgcttttcactttctatcagacaatctgtgtgagcactacaaagtacgcttctttaaggtaaggaggtgatccAACCg                                                                                                 <  1:870246‑M1/5‑1 (MQ=255)
                      gccttgcttttcactttctatcagacaatctgtgtgagcactacaaagtacgcttctttaaggtaaggaggtgatccAACCGTAGGGGAACCTGCGGTTg                                                                               >  1:2097463‑M1/78‑100 (MQ=255)
                        ttttactgctcattttcatcagacaatctgtgtgagcactgcaaagtacgcttctttaaggtaaggaggtgatccAACCGTAGGGGAACCTGCGGTTGGa                                                                             >  1:1504898‑M1/76‑100 (MQ=255)
                        cttgcttttcactttctatcagacaatctgtgtgagcactacaaagtacgcttctttaaggtaaggaggtgatccAACCGTAGGGGAACCTGCGGTTGGa                                                                             >  1:125694‑M1/76‑100 (MQ=255)
                         ttgcttttcactttctatcagacaatctgtgtgagcactgcaaagaacgcttctttaaggtaaggaggtgatccAACCGTAGGGGAACCTGCGGTTGGAt                                                                            >  1:36325‑M1/75‑100 (MQ=255)
                          ttactgctcattttcatcagacaatctgtgtgagcactgcaaagtacgcttctttaaggtaaggaggtgatccAACCGTAGGGGAACCTGCGGTTGGATc                                                                           >  2:68080‑M1/74‑100 (MQ=255)
                          tgcttttcactttctatcagacaatctgtgtgagcactacaaagtacgcttctttaaggtaaggaggtgatccAACCGTAGGGGAACCTGCGGTTGGATc                                                                           >  1:138234‑M1/74‑100 (MQ=255)
                           gcttttcactttctatcagacaatctgtgtgagcactacaaagtacgcttctttaaggtaaggaggtgatccAACCGTAGGGGAACCTGCGGTTGGATCa                                                                          >  1:2279999‑M1/73‑100 (MQ=255)
                            actgctcattttcatcagacaatctgtgtgagcactgcaaaggacgcttctttaaggtaaggaggtgatccAACCGTAGGGGAACCTGCGGTTGGATc                                                                           >  2:651261‑M1/72‑98 (MQ=255)
                             ttttcactttctatcagacaatctgtgtgagcacttcaaagaacgcttctttaaggtaagtaggtgatccAACCGTAGGGAAACCTGCGGTTGGATCa                                                                          >  2:1047400‑M1/71‑98 (MQ=255)
                                ctcattttcatcagacaatctgtgtgagcactgcaaagtacgcggctttaaggtaaggaggtgatccAACCGTAGGGGAACCTGCGGTTGGATCAcctcc                                                                     >  2:453182‑M1/68‑100 (MQ=255)
                                 cactttctatcagacaatctgtgtgagcactacaaagtacgcttctttaaggtaaggaggtgatccAACCGTAGGGGAACCTGCGGTTGGATCAcctcct                                                                    >  1:1982763‑M1/67‑100 (MQ=255)
                                                        tACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGTTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGt                                             <  2:99993/100‑1 (MQ=11)
                                                        tACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGt                                             <  2:1212062/100‑1 (MQ=11)
                                                        tACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGt                                             >  2:1580333/1‑100 (MQ=11)
                                                        tACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGt                                             <  2:1601495/100‑1 (MQ=11)
                                                          ccACTTGGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAg                                           <  1:35686/100‑1 (MQ=17)
                                                             cTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGc                                        >  2:1156408/1‑100 (MQ=33)
                                                             cTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGc                                        >  1:309460/1‑100 (MQ=32)
                                                             cTTTGTGATTCACGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGc                                        >  1:1570504/1‑100 (MQ=17)
                                                               ttGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTc                                      <  1:274326/100‑1 (MQ=33)
                                                               ttGTGATTCATGACTGGAGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTAACTTAAAGAAGCGTACTTTGTAATGCTc                                      >  2:1493438/1‑100 (MQ=255)
                                                                tgtgATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTca                                     >  1:479549/1‑100 (MQ=33)
                                                                tgtgATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTca                                     <  2:989727/100‑1 (MQ=33)
                                                                tgtgATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTca                                     >  2:661511/1‑100 (MQ=33)
                                                                tgtgATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTca                                     <  1:1184640/100‑1 (MQ=33)
                                                                tgtgATTCATGACTGGGGTGAAGTAGGAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTc                                      <  1:1559803/99‑1 (MQ=12)
                                                                 gtgATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcac                                    >  1:1782253/1‑100 (MQ=33)
                                                                 gtgATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcac                                    <  1:1143401/100‑1 (MQ=33)
                                                                 caaagaacgcttctttaaggtaaggaggtgatccAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGAAGTGCTcac                                    >  1:368774‑M1/35‑100 (MQ=255)
                                                                  tgATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcaca                                   <  1:1475467/100‑1 (MQ=33)
                                                                   gATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcacac                                  >  2:693751/1‑100 (MQ=33)
                                                                      tCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGa                               <  1:588141/100‑1 (MQ=33)
                                                                       cATGACTGGGGTGAAGTCGTAACAAGGTACCCGTAGGGGAACCTGCGGTTGGATCACCTACTTACCTTAAAGAAGCGTACTTTGTAGTGCTcacac                                  >  2:1727579/1‑96 (MQ=255)
                                                                       cATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAt                              >  2:458241/1‑100 (MQ=33)
                                                                       cATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACATTAAAGAAGCGTACTTTGTAGTGCTCACACAGAt                              >  2:869085/1‑100 (MQ=17)
                                                                           aCTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTc                          <  2:1985800/100‑1 (MQ=33)
                                                                           aCTGGGGTGAAGTCGTAACAAGGTAACCGTAGGACAAACTGCGGTTGGATCACCTACTTACCTTAAAGAAGCGTACTTTGTAGTGCTtacaca                                 <  2:1419499/93‑1 (MQ=11)
                                                                            cTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                         >  1:944949/1‑100 (MQ=33)
                                                                            cTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                         <  1:119013/100‑1 (MQ=33)
                                                                            cTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                         >  2:225208/1‑100 (MQ=33)
                                                                            cTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                         <  2:784584/100‑1 (MQ=33)
                                                                            cTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                         <  1:1059803/100‑1 (MQ=33)
                                                                              ggggTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGa                       >  2:1162213/1‑100 (MQ=33)
                                                                              ggggTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGa                       >  2:1997789/1‑100 (MQ=32)
                                                                                ggTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATa                     >  2:1041765/1‑100 (MQ=35)
                                                                                 gTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAg                    <  2:1875897/100‑1 (MQ=35)
                                                                                   gAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGaa                  >  1:2304616/1‑100 (MQ=38)
                                                                                     aGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAg                >  1:1182501/1‑100 (MQ=38)
                                                                                      gTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGt               >  2:533730/1‑100 (MQ=38)
                                                                                      gTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGt               >  1:636430/1‑100 (MQ=38)
                                                                                      gTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGt               >  1:1113394/1‑100 (MQ=38)
                                                                                      gTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGt               >  1:1236959/1‑100 (MQ=38)
                                                                                      gTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGt               >  1:1929312/1‑100 (MQ=38)
                                                                                        cGTTACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGa             >  2:325846/1‑100 (MQ=38)
                                                                                        cGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGa             >  2:547573/1‑100 (MQ=38)
                                                                                        cGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGa             >  1:2044358/1‑100 (MQ=38)
                                                                                        cGTAACAAGGTAACCGCAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTTATAGAAAGTGa             >  2:1310347/1‑100 (MQ=255)
                                                                                         gTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGaa            >  2:1217072/1‑100 (MQ=38)
                                                                                         gTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGaa            >  1:691466/1‑100 (MQ=38)
                                                                                         gTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCCTACTTAGTAGCGCTCACACAGATTGTCTGAAAGAAAGTGaa            >  1:1399586/1‑100 (MQ=11)
                                                                                          tAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGaaa           <  1:2121627/100‑1 (MQ=35)
                                                                                           aaCAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGaaaa          <  2:2334618/100‑1 (MQ=38)
                                                                                           aaCAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGaaaa          <  1:1860622/100‑1 (MQ=38)
                                                                                            aCAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAg         <  2:1704348/100‑1 (MQ=38)
                                                                                               aGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCaa      >  1:793265/1‑100 (MQ=38)
                                                                                                ggTAACCGTAGGGGAACCTGCGGTTGGATCACCGCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAg     <  1:375521/100‑1 (MQ=38)
                                                                                                 gTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGCTTGTCTGATAGAAAGTGAAAAGCAAgg    >  2:174510/1‑100 (MQ=38)
                                                                                                 gTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAgg    >  1:1511197/1‑100 (MQ=38)
                                                                                                   AACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCg  <  1:388638‑M1/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ACACCATGGGAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCG  >  NC_000913/4168066‑4168264

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.