Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 4,128,293 IS5 (+) +4 bp coding (100‑103/318 nt) metJ ← transcriptional repressor, S‑adenosylmethionine‑binding

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2290113NA (NA)51 (0.800) 49/298 0.5 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 = 4128296 0 (0.000)coding (100/318 nt) metJ transcriptional repressor, S‑adenosylmethionine‑binding
* ? NC_000913 3652036 =NA (NA)74 (1.160) 52/298 0.4 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 4128293 = 0 (0.000)coding (103/318 nt) metJ transcriptional repressor, S‑adenosylmethionine‑binding

GTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCC  >  NC_000913/3651886‑3652176
                                                                                                                                                      |                                                                                                                                            
gTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAgg                                                                                                                                              >  1:970268/1‑151 (MQ=255)
 tGAAAATAAAAACAGGATTAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACATGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAgg                                                                                                                                              <  2:568429/150‑1 (MQ=255)
 tGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCAAGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAgg                                                                                                                                              <  1:322221/150‑1 (MQ=255)
      atttaaaCAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAgg                                                                                                                                         >  2:60060/5‑150 (MQ=255)
        aaaaaCAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTg                                                                                                                                       <  2:108461/150‑1 (MQ=255)
          aaaCAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGa                                                                                                                                    <  1:1007897/151‑1 (MQ=255)
          aaaCAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGa                                                                                                                                    <  1:322311/151‑1 (MQ=255)
           aaCAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGaa                                                                                                                                   <  1:680693/151‑1 (MQ=255)
            acacgctGACAGCCTTTGACGGAACTTAATTTACTGATTTGAATCTGACAGGACAGTCCCTGTCCGGCTTCGTTCTAACCTTGTCAAGCTTTTATTGCGGAGATCTATTTAATATAAATTTCACTCCAGCCGATTTAGGGAAGGTGccagt                                                                                                                                  >  2:570296/7‑147 (MQ=255)
            aCAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAAt                                                                                                                                  <  1:134609/151‑1 (MQ=255)
            aCAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAAt                                                                                                                                  <  1:973624/151‑1 (MQ=255)
            aCAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAAt                                                                                                                                  <  2:680145/151‑1 (MQ=255)
                 aTGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCg                                                                                                                             >  2:470634/1‑151 (MQ=255)
                          tttGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAAttc                                                                                                                    >  1:343890/1‑151 (MQ=255)
                             gCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAAttcttc                                                                                                                 <  2:343840/151‑1 (MQ=255)
                              ccGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATtcttct                                                                                                                >  2:759418/1‑151 (MQ=255)
                                      aaCTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGAc                                                                                                        <  2:525030/151‑1 (MQ=255)
                                         tttATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTc                                                                                                      >  1:65368/1‑150 (MQ=255)
                                          ttATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTc                                                                                                      <  1:584020/149‑1 (MQ=255)
                                           tATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAg                                                                                                    <  2:167847/150‑1 (MQ=255)
                                           tATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAg                                                                                                    <  2:996929/150‑1 (MQ=255)
                                           tATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt                                                                                                   >  2:254722/1‑151 (MQ=255)
                                            aTTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt                                                                                                   <  2:246931/150‑1 (MQ=255)
                                             ttGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtc                                                                                                  >  1:692059/1‑150 (MQ=255)
                                               gATTTCAATATGACAGGACAGAGCCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcat                                                                                                <  1:810888/150‑1 (MQ=255)
                                                aTTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcattt                                                                                              >  1:777267/1‑151 (MQ=255)
                                                   tCAATATGACAGGACAGAACCTGTCCGGGTGCGTTCTAACCTTGTCAAGGTTTTATTTAGAAGATCTACTTAATATAAATTTCCCTGCCGCTAATTTAGGGAATGTGCGAATAAGCGGGGAAATTCTTCTCGACTGACTCATtcatttca                                                                                            <  2:733649/150‑1 (MQ=37)
                                                    cAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttcatt                                                                                          >  1:514476/1‑151 (MQ=255)
                                                          gacaggacagAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCa                                                                                    <  2:338933/151‑1 (MQ=255)
                                                           acaggacagAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCAt                                                                                   <  1:942427/151‑1 (MQ=255)
                                                           acaggacagAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCAt                                                                                   <  1:974386/151‑1 (MQ=255)
                                                            caggacagAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATg                                                                                  >  1:344425/1‑151 (MQ=255)
                                                             aggacagAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAAGTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGt                                                                                 <  1:557803/151‑1 (MQ=255)
                                                              ggacagAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGtt                                                                                <  1:307564/151‑1 (MQ=255)
                                                                 cagAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGa                                                                             <  1:665077/151‑1 (MQ=255)
                                                                  agAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAg                                                                            <  1:255541/151‑1 (MQ=255)
                                                                         gTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAttt                                                                     >  1:194831/1‑151 (MQ=255)
                                                                                  cgttcTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCg                                                            <  2:1001027/151‑1 (MQ=255)
                                                                                       tAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAt                                                       >  2:678453/1‑151 (MQ=255)
                                                                                        aaCCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAt                                                       >  2:295411/1‑150 (MQ=255)
                                                                                             tGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCtt                                                  <  2:123409/150‑1 (MQ=255)
                                                                                                cAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAAt                                              >  2:389989/1‑151 (MQ=255)
                                                                                                 aaGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATc                                             >  1:155790/1‑151 (MQ=255)
                                                                                                    gTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCa                                            >  2:419199/1‑149 (MQ=255)
                                                                                                    gTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                          <  2:143162/151‑1 (MQ=255)
                                                                                                         aTTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCGTTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAt                                      <  1:439373/150‑1 (MQ=255)
                                                                                                         aTTTCGAAGATATACTTAATATAAATTTCACTGACGCAAATTTAGGGTAGGTTCGAATAAGCGAAGAAATTCTTCTCGGCTGACTCAGTAATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTCAATGCCTTGAATCAGCCTAtt                                     >  2:422880/1‑151 (MQ=18)
                                                                                                              gAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAc                                >  1:371978/1‑151 (MQ=255)
                                                                                                                aGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCg                              >  1:711250/1‑151 (MQ=255)
                                                                                                                      aCTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttctt                        >  2:702381/1‑151 (MQ=255)
                                                                                                                        ttAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCg                      >  1:398075/1‑151 (MQ=255)
                                                                                                                          aaTATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGc                     <  1:208749/150‑1 (MQ=255)
                                                                                                                                aaTTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTaa              <  1:678539/151‑1 (MQ=40)
                                                                                                                                aaTTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTa               >  2:136641/1‑150 (MQ=255)
                                                                                                                                aaTTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTa               >  1:574112/1‑150 (MQ=255)
                                                                                                                                  tttCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAgg            >  2:410539/1‑151 (MQ=39)
                                                                                                                                   ttccatgaccgaaaTTTACGCATGGTACGAGTAAGCCAAGCAATTCTTGTCAGTTGACCCAGGCATTTAACTTCTTCATATTTAACCCGAGTATTTGTCGGGTTATTAACTTGTATCCGACTATTTAGACCGTTTCTTCGCCATTTAAGGc           >  2:816815/12‑151 (MQ=2)
                                                                                                                                     cACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg          <  2:609333/150‑1 (MQ=38)
                                                                                                                                     cACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg          <  2:691973/150‑1 (MQ=38)
                                                                                                                                     cACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGt         >  2:59268/1‑151 (MQ=37)
                                                                                                                                         gACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTa       <  2:228605/149‑1 (MQ=34)
                                                                                                                                         gACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATc     <  1:562610/151‑1 (MQ=34)
                                                                                                                                           cGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATccc   >  2:587/1‑151 (MQ=33)
                                                                                                                                              aaaTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc  >  1:631646/1‑149 (MQ=31)
                                                                                                                                                      |                                                                                                                                            
GTGAAAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCC  >  NC_000913/3651886‑3652176

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