Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 3,890,075 IS5 (–) +4 bp 100% coding (1346‑1349/1416 nt) tnaA → tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2290113NA (NA)97 (0.950) 60/278 NT 100% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 3890075 = 0 (0.000)coding (1346/1416 nt) tnaA tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
* ? NC_000913 3652036 =NA (NA)98 (0.960) 65/278 NT 100% noncoding (1195/1195 nt) IS5 repeat region
?NC_000913 = 3890078 0 (0.000)coding (1349/1416 nt) tnaA tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

AAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCA  >  NC_000913/3651890‑3652177
                                                                                                                                                  |                                                                                                                                             
aaaTAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGaa                                                                                                                                            <  1:1178591/150‑1 (MQ=255)
aaaTAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGaa                                                                                                                                            >  2:694082/1‑150 (MQ=255)
        aCAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGaa                                                                                                                                    >  1:236405/1‑150 (MQ=255)
          aGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATa                                                                                                                                  >  2:1177988/1‑150 (MQ=255)
            gATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAg                                                                                                                                >  2:1512064/1‑150 (MQ=255)
              tGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCg                                                                                                                              <  1:1465201/150‑1 (MQ=255)
               gAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCgg                                                                                                                             >  2:598606/1‑150 (MQ=255)
                    tCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGaaa                                                                                                                        >  1:251904/1‑150 (MQ=255)
                      tttGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAAtt                                                                                                                      >  1:398486/1‑150 (MQ=255)
                          ccGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAAttcttc                                                                                                                  >  1:1726533/1‑150 (MQ=255)
                          ccGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAAttcttc                                                                                                                  >  1:184071/1‑150 (MQ=255)
                          ccGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAAttcttc                                                                                                                  <  2:124477/150‑1 (MQ=255)
                          ccGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGGCAAGGTTGTATTTCGAAGATAGACTGAATATAAATTTCACTGACGCAAATTTAGGGAAGTTGCGAATCAGCGGGGAAAttctta                                                                                                                  >  1:184048/1‑149 (MQ=255)
                           cGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATtcttct                                                                                                                 <  1:168197/150‑1 (MQ=255)
                             gAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCg                                                                                                               <  2:969157/150‑1 (MQ=255)
                              aaCCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCgg                                                                                                              <  1:679047/150‑1 (MQ=255)
                                ccAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCt                                                                                                            <  2:225128/150‑1 (MQ=255)
                                ccAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCt                                                                                                            <  2:225242/150‑1 (MQ=255)
                                  aaCTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGa                                                                                                          >  2:194958/1‑150 (MQ=255)
                                   aCTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGAc                                                                                                         <  2:1393592/150‑1 (MQ=255)
                                    cTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACt                                                                                                        >  2:1593570/1‑150 (MQ=255)
                                    cTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACt                                                                                                        >  2:325664/1‑150 (MQ=255)
                                        aTTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAATGTGAGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt                                                                                                    <  1:990074/150‑1 (MQ=255)
                                        aTTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt                                                                                                    >  1:1036379/1‑150 (MQ=255)
                                                    atGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttct                                                                                        <  1:739363/150‑1 (MQ=255)
                                                    atGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttct                                                                                        <  1:653682/150‑1 (MQ=255)
                                                    atGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttct                                                                                        <  1:739362/150‑1 (MQ=255)
                                                     tGACAGGACAGAACCTGTCCGGGTTCTTTCTAACCTTGTCTATGTTTTATTTCGAAGATCTAATTAATATAAATTTCTCTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttctt                                                                                       <  1:521373/150‑1 (MQ=37)
                                                     tGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGTAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttctt                                                                                       <  1:59392/150‑1 (MQ=255)
                                                     tGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttctt                                                                                       <  2:169647/150‑1 (MQ=255)
                                                     tGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttctt                                                                                       <  2:1217981/150‑1 (MQ=255)
                                                      gacaggacagAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttcttc                                                                                      >  1:29533/1‑150 (MQ=255)
                                                      gacaggacagAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATttcttc                                                                                      >  2:475239/1‑150 (MQ=255)
                                                       acaggacatACCTTGTCCGGGTTGTTTCTACCCTTGTCATGGTTTTATTTCGATGATCTACTTATTATAAATTTCACTTACGAAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCa                                                                                     <  2:864433/150‑1 (MQ=21)
                                                       acaggacagAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCa                                                                                     >  2:286022/1‑150 (MQ=255)
                                                               gAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAg                                                                             >  1:1291151/1‑150 (MQ=255)
                                                               gAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAg                                                                             >  1:846426/1‑150 (MQ=255)
                                                               gAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAg                                                                             >  1:1291358/1‑150 (MQ=255)
                                                               gAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAg                                                                             <  2:398486/150‑1 (MQ=255)
                                                                aaCCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGc                                                                            >  2:1395994/1‑150 (MQ=255)
                                                                aaCCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGc                                                                            >  1:365659/1‑150 (MQ=255)
                                                                   cTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGa                                                                         >  2:1527110/1‑150 (MQ=255)
                                                                         gggTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                                                                   <  2:56269/150‑1 (MQ=255)
                                                                         gggTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttttt                                                                   >  2:960310/1‑150 (MQ=255)
                                                                           gttcgttcTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtct                                                                 >  2:502938/1‑150 (MQ=255)
                                                                           gttcgttcTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtct                                                                 <  2:373207/150‑1 (MQ=255)
                                                                            ttcgttcTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTtctc                                                                <  1:1288889/150‑1 (MQ=255)
                                                                             tcgttcTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTcc                                                               >  1:618893/1‑150 (MQ=255)
                                                                             tcgttcTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTcc                                                               <  1:326869/150‑1 (MQ=255)
                                                                              cgttcTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTccc                                                              <  2:236405/150‑1 (MQ=255)
                                                                                ttcTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGt                                                            >  2:863348/1‑150 (MQ=255)
                                                                                 tcTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTa                                                           >  1:1273823/1‑150 (MQ=255)
                                                                                  cTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaa                                                          >  2:912814/1‑150 (MQ=255)
                                                                                  cTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaa                                                          >  1:528985/1‑150 (MQ=255)
                                                                                  cTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaa                                                          >  1:531138/1‑150 (MQ=255)
                                                                                   tAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGTGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaaa                                                         >  1:1157544/1‑150 (MQ=255)
                                                                                    aaCCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAAt                                                        <  2:579180/150‑1 (MQ=255)
                                                                                      ccTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAATTGc                                                      >  2:79339/1‑150 (MQ=255)
                                                                                      ccTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc                                                      >  2:1423892/1‑150 (MQ=255)
                                                                                      ccTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc                                                      <  2:400158/150‑1 (MQ=255)
                                                                                        ttGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCt                                                    >  2:578140/1‑150 (MQ=255)
                                                                                               ggTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCa                                             <  1:694082/150‑1 (MQ=255)
                                                                                                 ttttATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGc                                           <  2:251904/150‑1 (MQ=255)
                                                                                                     aTTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAt                                       >  2:914180/1‑150 (MQ=255)
                                                                                                         cGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAg                                   <  2:1270881/150‑1 (MQ=255)
                                                                                                         cGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAg                                   <  1:627417/150‑1 (MQ=255)
                                                                                                          gAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGa                                  <  2:942310/150‑1 (MQ=255)
                                                                                                                   cTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttctt                         >  1:312748/1‑150 (MQ=255)
                                                                                                                   cTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttctt                         >  1:312756/1‑150 (MQ=255)
                                                                                                                     tAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCg                       <  2:1291358/150‑1 (MQ=255)
                                                                                                                     tAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCg                       <  1:607855/150‑1 (MQ=255)
                                                                                                                     tAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCg                       <  2:1291151/150‑1 (MQ=255)
                                                                                                                       atatAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGcc                     <  1:598606/150‑1 (MQ=255)
                                                                                                                         ataAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAt                   >  2:721749/1‑150 (MQ=255)
                                                                                                                          taAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAtt                  >  2:1126536/1‑150 (MQ=255)
                                                                                                                              tttCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAg              <  1:902490/150‑1 (MQ=40)
                                                                                                                              tttCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAg              <  1:502938/150‑1 (MQ=40)
                                                                                                                                 cACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg           <  1:897464/150‑1 (MQ=38)
                                                                                                                                 cACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg           <  1:706834/150‑1 (MQ=38)
                                                                                                                                 cACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg           <  1:806008/150‑1 (MQ=38)
                                                                                                                                 cACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg           <  2:1139085/150‑1 (MQ=38)
                                                                                                                                 cACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg           <  1:1279442/150‑1 (MQ=38)
                                                                                                                                 cACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg           <  1:1214135/150‑1 (MQ=38)
                                                                                                                                 cACTCACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg           <  2:710142/150‑1 (MQ=25)
                                                                                                                                   cTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGtt         <  1:194958/150‑1 (MQ=37)
                                                                                                                                    tGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTa        >  2:934083/1‑150 (MQ=35)
                                                                                                                                     gACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTAt       >  2:968378/1‑150 (MQ=35)
                                                                                                                                          aaaTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCa  <  2:645951/150‑1 (MQ=31)
                                                                                                                                          aaaTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCa  <  1:279550/150‑1 (MQ=31)
                                                                                                                                          aaaTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCa  <  1:759014/150‑1 (MQ=31)
                                                                                                                                                  |                                                                                                                                             
AAATAAAAACAGGATGAAAGTCTTTGCCGGAACCAACTTTATTGATTTCAATATGACAGGACAGAACCTGTCCGGGTTCGTTCTAACCTTGTCAAGGTTTTATTTCGAAGATCTACTTAATATAAATTTCACTGACGCAAATTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCA  >  NC_000913/3651890‑3652177

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