Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,750,123 C→T 2.4% G122S (GGC→AGC)  tnaA ← tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,750,1230C→T 2.4% 12.7 ‑1.7 ‑0.2368G122S (GGC→AGC) tnaAtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  new base (6/3):  ref base (155/204):  total (161/207)
Fisher's exact test strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.81e-01

Only 100 of 396 total aligned reads displayed.

ATAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGATATTT  >  NC_007779/3750025‑3750223
                                                                                                  |                                                                                                    
ataGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTtt                                                                                                    >  1:3235618/1‑101 (MQ=255)
ataGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTtt                                                                                                    <  1:6778056/101‑1 (MQ=255)
 taGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTttt                                                                                                   >  1:16761355/1‑101 (MQ=255)
  aGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTtttt                                                                                                  >  1:15717564/1‑101 (MQ=255)
     cGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCt                                                                                              >  1:1651486/1‑101 (MQ=255)
     cGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCt                                                                                              <  1:11019746/101‑1 (MQ=255)
     cGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCt                                                                                              <  1:1965855/101‑1 (MQ=255)
      gTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTg                                                                                             >  1:6064970/1‑101 (MQ=255)
        tACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCt                                                                                           <  1:15968915/101‑1 (MQ=255)
        tACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCt                                                                                           >  1:2739765/1‑101 (MQ=255)
        tACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCt                                                                                           <  1:10182806/101‑1 (MQ=255)
          cGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcg                                                                                         <  1:17575220/101‑1 (MQ=255)
           gCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcgc                                                                                        >  1:5896639/1‑101 (MQ=255)
           gCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcgc                                                                                        >  1:15336082/1‑101 (MQ=255)
               ggTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGttt                                                                                    >  1:13323232/1‑101 (MQ=255)
                   cAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTa                                                                                <  1:14383195/101‑1 (MQ=255)
                      ccGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATc                                                                             >  1:17178076/1‑101 (MQ=255)
                        gTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAg                                                                           <  1:6267728/101‑1 (MQ=255)
                             tCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCg                                                                      >  1:7988135/1‑101 (MQ=255)
                              cTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCgg                                                                     <  1:2938950/101‑1 (MQ=255)
                               tGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGa                                                                    <  1:6874955/101‑1 (MQ=255)
                                    aTGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAAtata                                                               <  1:14004769/101‑1 (MQ=255)
                                     tGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAg                                                              >  1:16152336/1‑101 (MQ=255)
                                      ggCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGa                                                             >  1:12600274/1‑101 (MQ=255)
                                       gCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGAt                                                            >  1:2474937/1‑101 (MQ=255)
                                       gCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCATGCTCGCGTTTTTTAATCAGTACCGGAATATAGAt                                                            >  1:4516444/1‑101 (MQ=255)
                                        cccTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGAtt                                                           <  1:256789/101‑1 (MQ=255)
                                         ccTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGAttt                                                          <  1:8096778/101‑1 (MQ=255)
                                         ccTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGAttt                                                          <  1:1198994/101‑1 (MQ=255)
                                          cTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTg                                                         >  1:5465497/1‑101 (MQ=255)
                                           tGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGc                                                        <  1:17629863/101‑1 (MQ=255)
                                            gCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGct                                                       <  1:12871493/101‑1 (MQ=255)
                                            gCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGct                                                       <  1:9117892/101‑1 (MQ=255)
                                            gCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGct                                                       <  1:11873521/101‑1 (MQ=255)
                                              gTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctct                                                     <  1:7466967/101‑1 (MQ=255)
                                              gTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctct                                                     >  1:12677398/1‑101 (MQ=255)
                                               tGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTg                                                    >  1:153650/1‑101 (MQ=255)
                                               tGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTg                                                    >  1:3349117/1‑101 (MQ=255)
                                               tGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTg                                                    >  1:4848521/1‑101 (MQ=255)
                                                ggTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTgc                                                   <  1:12161412/101‑1 (MQ=255)
                                                 gTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTgcg                                                  <  1:13062175/101‑1 (MQ=255)
                                                 gTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTgcg                                                  <  1:9665860/101‑1 (MQ=255)
                                                  tATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTgcgc                                                 <  1:10691505/101‑1 (MQ=255)
                                                    tCAAAGAAATAGTTAGAGAACTCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCa                                               <  1:1809571/101‑1 (MQ=255)
                                                    tCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCa                                               >  1:17447525/1‑101 (MQ=255)
                                                    tCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCa                                               <  1:17796777/101‑1 (MQ=255)
                                                     cAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCAc                                              <  1:2210672/101‑1 (MQ=255)
                                                         gAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGcc                                          <  1:4620123/101‑1 (MQ=255)
                                                          aaaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGccc                                         >  1:5643001/1‑101 (MQ=255)
                                                          aaaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGccc                                         >  1:18236021/1‑101 (MQ=255)
                                                           aaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCt                                        <  1:5013652/101‑1 (MQ=255)
                                                           aaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCt                                        <  1:4027782/101‑1 (MQ=255)
                                                           aaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCt                                        <  1:15058314/101‑1 (MQ=255)
                                                           aaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCt                                        <  1:755057/101‑1 (MQ=255)
                                                           aaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCt                                        <  1:11309926/101‑1 (MQ=255)
                                                           aaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCGTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCt                                        <  1:5234499/101‑1 (MQ=255)
                                                            aTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTg                                       >  1:4196997/1‑101 (MQ=255)
                                                              aGTTAGAGAACGCCACCATTTTTCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTggg                                     >  1:15499105/1‑100 (MQ=255)
                                                              aGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGTTCGCGTTTTTTAATCAGTACCGGAATATAGATTTCCTCTGCGCCACGGCCCTGGt                                     <  1:11193889/101‑1 (MQ=255)
                                                              aGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGt                                     <  1:17691075/101‑1 (MQ=255)
                                                               gTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTg                                    <  1:6019801/101‑1 (MQ=255)
                                                               gTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTg                                    <  1:9735178/101‑1 (MQ=255)
                                                                ttAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGa                                   <  1:14684295/101‑1 (MQ=255)
                                                                ttAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGa                                   >  1:13717217/1‑101 (MQ=255)
                                                                 tAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAg                                  <  1:17906139/101‑1 (MQ=255)
                                                                   gagaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTc                                >  1:1042478/1‑101 (MQ=255)
                                                                      aaCGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGa                             >  1:17227307/1‑101 (MQ=255)
                                                                       aCGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGaa                            >  1:12389909/1‑101 (MQ=255)
                                                                       aCGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGaa                            <  1:8762668/101‑1 (MQ=255)
                                                                        cGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAAt                           >  1:8575319/1‑101 (MQ=255)
                                                                         gCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATg                          >  1:12089168/1‑101 (MQ=255)
                                                                          ccaccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATgg                         >  1:9271059/1‑101 (MQ=255)
                                                                           cacctTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgt                        >  1:8440751/1‑101 (MQ=255)
                                                                           caccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgt                        <  1:3502389/101‑1 (MQ=255)
                                                                            accaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgtg                       <  1:5336347/101‑1 (MQ=255)
                                                                             ccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgtgt                      <  1:7134014/101‑1 (MQ=255)
                                                                               aTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAt                    >  1:4407749/1‑101 (MQ=255)
                                                                                ttttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAtt                   >  1:3022019/1‑101 (MQ=255)
                                                                                atttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGTCCCTGGTGATTCGGAATGGTGTAtt                   <  1:14631568/100‑1 (MQ=255)
                                                                                  ttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGa                 <  1:9019727/101‑1 (MQ=255)
                                                                                   tgctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGAt                <  1:11559952/101‑1 (MQ=255)
                                                                                   tgctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGAt                <  1:1064785/101‑1 (MQ=255)
                                                                                    gctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATa               <  1:8243629/101‑1 (MQ=255)
                                                                                     ctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATaa              <  1:14239731/101‑1 (MQ=255)
                                                                                      tgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAAc             <  1:9399802/101‑1 (MQ=255)
                                                                                       gcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAAcc            >  1:7833526/1‑101 (MQ=255)
                                                                                        cGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCa           >  1:13204909/1‑101 (MQ=255)
                                                                                         gATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaa          >  1:13852224/1‑101 (MQ=255)
                                                                                         gATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaa          >  1:4773408/1‑101 (MQ=255)
                                                                                          aTCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaaa         <  1:3964038/101‑1 (MQ=255)
                                                                                           tCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAg        <  1:6583629/101‑1 (MQ=255)
                                                                                           tCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAg        >  1:5536624/1‑101 (MQ=255)
                                                                                              caggcTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCCCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGat      <  1:4232387/97‑1 (MQ=255)
                                                                                               aggcTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGata     >  1:490912/4‑101 (MQ=255)
                                                                                                 gcTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGATAtt   <  1:3700876/100‑1 (MQ=255)
                                                                                                 gcTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGATAtt   <  1:10046090/100‑1 (MQ=255)
                                                                                                  cTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGATAttt  <  1:6347242/101‑1 (MQ=255)
                                                                                                  cTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGATAttt  >  1:17541846/1‑101 (MQ=255)
                                                                                                  cTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGATAttt  >  1:2860527/1‑101 (MQ=255)
                                                                                                  cTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGATAttt  <  1:3629786/101‑1 (MQ=255)
                                                                                                  |                                                                                                    
ATAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGATATTT  >  NC_007779/3750025‑3750223

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