New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3694553 =85 (0.810)4 (0.040)
+C
4/258 NT 4.5% coding (684/2619 nt) bcsA cellulose synthase, catalytic subunit
?NC_000913 3844657 = 86 (0.820)coding (694/1767 nt) adeD cryptic adenine deaminase
Rejected: Frequency below/above cutoff threshold.

TTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGC  >  NC_000913/3844511‑3844804
                                                                                                                                                  |                                                                                                                                                   
ttGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGaa                                                                                                                                                   >  3:714664/1‑149 (MQ=255)
ttGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGa                                                                                                                                                    >  1:866334/1‑148 (MQ=255)
       aaTGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGcc                                                                                                                                            >  4:164461/1‑149 (MQ=255)
             ggACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGaaa                                                                                                                                      <  3:837252/149‑1 (MQ=255)
              gACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAg                                                                                                                                     <  3:387794/149‑1 (MQ=255)
                cTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCAGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGtt                                                                                                                                   <  4:96300/149‑1 (MQ=255)
                      tGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGc                                                                                                                             >  4:829509/1‑149 (MQ=255)
                       ggCGTAATTAGCGGTCAGAATGCGCTGCTGGATAAACTGGATGCATTTCGCCACCTGACCCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGc                                                                                                                             <  1:696554/148‑1 (MQ=255)
                       ggCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCt                                                                                                                            <  4:689652/149‑1 (MQ=255)
                           tAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGga                                                                                                                         <  1:159955/148‑1 (MQ=255)
                             aTTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGgaaga                                                                                                                      >  3:727948/1‑149 (MQ=255)
                              ttAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGgaaga                                                                                                                      <  2:258834/148‑1 (MQ=255)
                               tAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGaagaag                                                                                                                    <  3:164461/149‑1 (MQ=255)
                               tAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGaagaag                                                                                                                    <  3:60558/149‑1 (MQ=255)
                                  cGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGAc                                                                                                                 <  3:722668/149‑1 (MQ=255)
                                    ggCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACg                                                                                                                <  2:229767/148‑1 (MQ=255)
                                    ggCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGc                                                                                                               >  3:937957/1‑149 (MQ=255)
                                        gAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGa                                                                                                           >  4:532180/1‑149 (MQ=255)
                                         aaTGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGa                                                                                                           <  1:950160/148‑1 (MQ=255)
                                             cgcTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAAtt                                                                                                       <  2:833178/148‑1 (MQ=255)
                                                   tCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTc                                                                                                >  4:406082/1‑149 (MQ=255)
                                                    cGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAAtt                                                                                                       >  1:263247/1‑141 (MQ=255)
                                                    cGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAAtt                                                                                                       <  2:263247/141‑1 (MQ=255)
                                                      aTAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCgg                                                                                              >  2:840197/1‑148 (MQ=255)
                                                      aTAAACGGGATGCATTTCGCCACCTGACGCTTGACGGTAACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGGAAACTGCCACGAAAGTTATCAGCTGTAAGAAGGACGCCGGAAATTACAACTCAGc                                                                                             >  4:735036/1‑149 (MQ=255)
                                                          aCTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATg                                                                                          <  1:700041/148‑1 (MQ=255)
                                                           cTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTc                                                                                        >  4:596214/1‑149 (MQ=255)
                                                              gATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGt                                                                                      <  2:392484/148‑1 (MQ=255)
                                                                tGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtga                                                                                   <  3:38606/149‑1 (MQ=255)
                                                                tGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtg                                                                                    >  2:728322/1‑148 (MQ=255)
                                                                 gCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtga                                                                                   <  1:521085/148‑1 (MQ=255)
                                                                  cATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtgatg                                                                                 <  4:577888/149‑1 (MQ=255)
                                                                    tttCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtgatga                                                                                >  2:462733/1‑148 (MQ=255)
                                                                           aCCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGaa                                                                        >  3:256273/1‑149 (MQ=255)
                                                                           aCCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGaa                                                                        <  4:714664/149‑1 (MQ=255)
                                                                           aCCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGa                                                                         <  2:46803/148‑1 (MQ=255)
                                                                           aCCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCAGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGaa                                                                        <  3:523318/149‑1 (MQ=255)
                                                                             cTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAgg                                                                      >  4:707981/1‑149 (MQ=255)
                                                                             cTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAgg                                                                      >  1:257838/1‑148 (MQ=255)
                                                                              tGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAggg                                                                     >  4:434179/1‑149 (MQ=255)
                                                                               gACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATTAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGt                                                                    >  4:259083/1‑149 (MQ=255)
                                                                               gACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACg                                                                                                                <  2:615590/105‑1 (MQ=255)
                                                                               gACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACg                                                                                                                >  1:615590/1‑105 (MQ=255)
                                                                               gACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGt                                                                    <  3:652276/149‑1 (MQ=255)
                                                                                aCGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGt                                                                    >  2:796936/1‑148 (MQ=255)
                                                                                       aCGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGc                                                             >  2:566866/1‑148 (MQ=255)
                                                                                       aCGGTCACTGACCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGcc                                                            <  4:248670/149‑1 (MQ=255)
                                                                                          gTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGaaa                                                                                                         <  4:258039/101‑1 (MQ=255)
                                                                                          gTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGaaa                                                                                                         >  3:258039/1‑101 (MQ=255)
                                                                                           tCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGc                                                         >  2:907671/1‑148 (MQ=255)
                                                                                            cACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCa                                                        >  1:410789/1‑148 (MQ=255)
                                                                                             aCTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAAt                                                      >  3:569151/1‑149 (MQ=255)
                                                                                              cTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAAt                                                      <  2:707389/148‑1 (MQ=255)
                                                                                              cTGCCCGGGATTGCGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGGGAAGGGTCCGCTGCCCGCAAt                                                      >  2:497360/1‑148 (MQ=255)
                                                                                                gCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTAGTGAAATCTGCCACGTAAGTTATCAGCTGGAAGTAGGACGCCGGAAATTACAACTCGGCATGTCGGTGATGATCCGCGAAGGGTTCCCTGCCCGCAAtctc                                                   <  4:929430/149‑1 (MQ=255)
                                                                                                gCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAAtct                                                    >  2:806642/1‑148 (MQ=255)
                                                                                                 cccGGGTTTGGGTGGTAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCa                                                  <  2:257838/148‑1 (MQ=255)
                                                                                                    atctggtggcgttacacctctacgctgaactgggacgatccggtccGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGccc                                                           >  3:849525‑M2/47‑137 (MQ=255)
                                                                                                    atctggtggcgttacacctctacgctgaactgggacgatccggtccGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGccc                                                           <  4:849525‑M2/91‑1 (MQ=255)
                                                                                                     tctggtggcgttacacctctacgctgaactgggacgatccggtccGAAAACTGCCACGAAAGTTGTCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGcc                                                            >  1:245079‑M2/46‑135 (MQ=255)
                                                                                                     tctggtggcgttacacctctacgctgaactgggacgatccggtccGAAAACTGCCACGAAAGTTGTCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGcc                                                            <  2:245079‑M2/90‑1 (MQ=255)
                                                                                                        ttGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACgcgc                                            <  1:892116/148‑1 (MQ=255)
                                                                                                         tGGGTGGTAAAGAACTTCACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTg                                          <  4:536333/149‑1 (MQ=255)
                                                                                                         tGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCt                                           <  1:533725/148‑1 (MQ=255)
                                                                                                                  aaGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTg                                 >  3:35178/1‑149 (MQ=255)
                                                                                                                   aGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtga                                                                                   >  4:139736/1‑98 (MQ=255)
                                                                                                                   aGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtga                                                                                   <  3:139736/98‑1 (MQ=255)
                                                                                                                   aGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTg                                 >  1:518407/1‑148 (MQ=255)
                                                                                                                    gAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGa                                <  1:299482/148‑1 (MQ=255)
                                                                                                                    gAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGAt                               >  3:200400/1‑149 (MQ=255)
                                                                                                                     aaCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGAt                               <  1:441158/148‑1 (MQ=255)
                                                                                                                     aaCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGAt                               <  2:732619/148‑1 (MQ=255)
                                                                                                                       cTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCa                             <  1:120436/148‑1 (MQ=255)
                                                                                                                       cTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCa                             <  1:711321/148‑1 (MQ=255)
                                                                                                                        ttAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGtt                                  >  4:275438/1‑142 (MQ=255)
                                                                                                                        ttAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGtt                                  <  3:275438/142‑1 (MQ=255)
                                                                                                                        ttAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCaa                            <  1:907671/148‑1 (MQ=255)
                                                                                                                        ttAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAAc                           <  4:569151/149‑1 (MQ=255)
                                                                                                                         tAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAAc                           >  1:652315/1‑148 (MQ=255)
                                                                                                                            cGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGaa                        <  1:444684/148‑1 (MQ=255)
                                                                                                                            cGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAAt                       <  3:759354/149‑1 (MQ=255)
                                                                                                                             gCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGccc                                                           >  3:824804/1‑112 (MQ=255)
                                                                                                                             gCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGccc                                                           <  4:824804/112‑1 (MQ=255)
                                                                                                                             gCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAAc                           <  3:662730/144‑1 (MQ=255)
                                                                                                                             gCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAAtt                      >  4:171905/1‑149 (MQ=255)
                                                                                                                             gCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAAGTACAACTCGGCATGTCGTTGATGATCCGCAAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAAc                           >  4:662730/1‑144 (MQ=255)
                                                                                                                              ccTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAAtt                      >  2:182893/1‑148 (MQ=255)
                                                                                                                               cTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTa                    >  3:83555/1‑149 (MQ=255)
                                                                                                                                  tatTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACa                 >  3:682604/1‑149 (MQ=255)
                                                                                                                                  tatTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATATCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACa                 <  4:727948/149‑1 (MQ=255)
                                                                                                                                          cGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCa          <  2:628082/148‑1 (MQ=255)
                                                                                                                                           gggTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAat        <  3:171905/149‑1 (MQ=255)
                                                                                                                                                ttGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAatgcatg   >  4:678401/1‑149 (MQ=255)
                                                                                                                                                ttGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAatgcatg   >  4:762135/1‑149 (MQ=255)
                                                                                                                                                ttGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAatgcat    <  1:231071/148‑1 (MQ=255)
                                                                                                                                                  gAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAatgcatgc  <  1:182893/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAATTTAACAGCCCGCAATGCATGC  >  NC_000913/3844511‑3844804

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