Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,750,123 C→T 4.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 4.4% 22.4 ‑0.2 ‑0.1297G122S (GGC→AGC) tnaAtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  new base (7/6):  ref base (137/145):  total (146/151)
Fisher's exact test strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.12e-01

Only 100 of 320 total aligned reads displayed.

ATAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGATATTT  >  NC_007779/3750025‑3750223
                                                                                                  |                                                                                                    
ataGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTtt                                                                                                    >  1:5794837/1‑101 (MQ=255)
ataGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTtt                                                                                                    >  1:1507493/1‑101 (MQ=255)
 taGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTttt                                                                                                   >  1:6483394/1‑101 (MQ=255)
   gACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTttttc                                                                                                 >  1:10864452/1‑100 (MQ=255)
     cGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCt                                                                                              <  1:6688230/101‑1 (MQ=255)
     cGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCt                                                                                              <  1:307564/101‑1 (MQ=255)
      gTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTg                                                                                             <  1:2481441/101‑1 (MQ=255)
        tACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCt                                                                                           <  1:13465411/101‑1 (MQ=255)
            cACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcgcg                                                                                       <  1:13206216/101‑1 (MQ=255)
             aCGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGt                                                                                      >  1:6168047/1‑101 (MQ=255)
             aCGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGt                                                                                      <  1:6026949/101‑1 (MQ=255)
              cGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtt                                                                                     >  1:12817572/1‑101 (MQ=255)
               ggTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGttt                                                                                    <  1:7078893/101‑1 (MQ=255)
               ggTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGttt                                                                                    >  1:10284564/1‑101 (MQ=255)
                gTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtttt                                                                                   >  1:9206946/1‑101 (MQ=255)
                  aCAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtttttt                                                                                 <  1:5234576/101‑1 (MQ=255)
                  aCAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtttttt                                                                                 <  1:672497/101‑1 (MQ=255)
                  aCAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtttttt                                                                                 >  1:3447835/1‑101 (MQ=255)
                   cAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTa                                                                                >  1:13350110/1‑101 (MQ=255)
                       cGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCa                                                                            >  1:5011020/1‑101 (MQ=255)
                       cGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCa                                                                            <  1:7535600/101‑1 (MQ=255)
                       cGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCa                                                                            >  1:9474172/1‑101 (MQ=255)
                       cGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCa                                                                            <  1:12004189/101‑1 (MQ=255)
                       cGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCa                                                                            >  1:11647318/1‑101 (MQ=255)
                        gTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAg                                                                           <  1:6252381/101‑1 (MQ=255)
                          tGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTa                                                                         <  1:2352233/101‑1 (MQ=255)
                            aTCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTAcc                                                                       >  1:1268081/1‑101 (MQ=255)
                            aTCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTAcc                                                                       >  1:12420879/1‑101 (MQ=255)
                              cTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCgg                                                                     <  1:1724587/101‑1 (MQ=255)
                               tGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGa                                                                    >  1:790540/1‑101 (MQ=255)
                                 gCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAat                                                                  >  1:12198053/1‑101 (MQ=255)
                                 gCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAat                                                                  <  1:4454615/101‑1 (MQ=255)
                                 gCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAat                                                                  <  1:2956517/101‑1 (MQ=255)
                                 gCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAat                                                                  >  1:10405041/1‑101 (MQ=255)
                                  cTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAata                                                                 <  1:7203583/101‑1 (MQ=255)
                                  cTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAata                                                                 <  1:10801245/101‑1 (MQ=255)
                                   tATGGCCCTGCGTGGGATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAatat                                                                <  1:9740704/101‑1 (MQ=255)
                                     tgtcccTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAg                                                              <  1:6395266/98‑1 (MQ=255)
                                     tGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAg                                                              <  1:4601256/101‑1 (MQ=255)
                                      ggCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGa                                                             >  1:1836681/1‑101 (MQ=255)
                                      ggCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGa                                                             <  1:5270800/101‑1 (MQ=255)
                                      ggCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGa                                                             <  1:10183401/101‑1 (MQ=255)
                                           tGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGct                                                       >  1:409532/1‑101 (MQ=255)
                                            gCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGct                                                       <  1:8157078/101‑1 (MQ=255)
                                             cGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctc                                                      >  1:858220/1‑101 (MQ=255)
                                             cGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTCTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctc                                                      >  1:11263391/1‑101 (MQ=255)
                                              gTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctct                                                     <  1:11788565/101‑1 (MQ=255)
                                              gTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctct                                                     <  1:822689/101‑1 (MQ=255)
                                               tGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTg                                                    >  1:11032825/1‑101 (MQ=255)
                                                gggaTCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGGATATAGATTTGCTCTgc                                                   >  1:11170105/4‑101 (MQ=255)
                                                 gTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTgcg                                                  >  1:8023781/1‑101 (MQ=255)
                                                  tATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTgcgc                                                 >  1:13794219/1‑101 (MQ=255)
                                                     cAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCAc                                              <  1:7678934/101‑1 (MQ=255)
                                                     cAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCAc                                              >  1:10043721/1‑101 (MQ=255)
                                                     cAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCAc                                              >  1:8804901/1‑101 (MQ=255)
                                                       aaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACgg                                            >  1:5396805/1‑101 (MQ=255)
                                                       aaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACgg                                            >  1:12755464/1‑101 (MQ=255)
                                                        aGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGc                                           <  1:139630/101‑1 (MQ=255)
                                                        aGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGc                                           <  1:14334955/101‑1 (MQ=255)
                                                         gAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGcc                                          >  1:1930129/1‑101 (MQ=255)
                                                           aaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCt                                        <  1:6881588/101‑1 (MQ=255)
                                                             tAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTgg                                      <  1:4293313/101‑1 (MQ=255)
                                                              aGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGt                                     <  1:2783213/101‑1 (MQ=255)
                                                               gTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTg                                    <  1:3037524/101‑1 (MQ=255)
                                                               gTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTg                                    <  1:4845777/101‑1 (MQ=255)
                                                                ttAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGa                                   >  1:13011470/1‑101 (MQ=255)
                                                                 tAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAg                                  >  1:4796788/1‑101 (MQ=255)
                                                                  agagAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGt                                 <  1:13587628/101‑1 (MQ=255)
                                                                  agagAACGCCACCAGTTTGCTGCGATCCAGGCTTTTTCGTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGt                                 <  1:10622733/101‑1 (MQ=255)
                                                                   gagaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTc                                <  1:262704/101‑1 (MQ=255)
                                                                     gaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCgg                              >  1:11947812/1‑101 (MQ=255)
                                                                     gaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCgg                              <  1:4915152/101‑1 (MQ=255)
                                                                     gaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCgg                              <  1:3643250/101‑1 (MQ=255)
                                                                       acaccaccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGaa                            >  1:11199591/4‑101 (MQ=255)
                                                                        cGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAAt                           <  1:14439824/101‑1 (MQ=255)
                                                                         gCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATg                          <  1:5904737/101‑1 (MQ=255)
                                                                         gCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATg                          <  1:8960978/101‑1 (MQ=255)
                                                                               aTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAt                    >  1:798646/1‑101 (MQ=255)
                                                                               aTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAt                    <  1:7628615/101‑1 (MQ=255)
                                                                                ttttGCTTGGATTCAGGCTTTTTTCTTGCTGCGGTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAtt                   <  1:9141648/101‑1 (MQ=255)
                                                                                 tttGGTGCGATTCAGGGTTTTTCCTGCTCGCGTTTTTTAATCAGTTACGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTg                  <  1:5058098/101‑1 (MQ=255)
                                                                                 tttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTg                  >  1:7159108/1‑101 (MQ=255)
                                                                                   tgctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGAt                <  1:3195253/101‑1 (MQ=255)
                                                                                   tgctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCAGAATGGTGTATTGAt                >  1:9872898/1‑101 (MQ=255)
                                                                                    gctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATa               <  1:2238002/101‑1 (MQ=255)
                                                                                     ctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATaa              <  1:3735719/101‑1 (MQ=255)
                                                                                     ctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATaa              <  1:2873485/101‑1 (MQ=255)
                                                                                      tgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAAc             >  1:3370818/1‑101 (MQ=255)
                                                                                         gATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaa          >  1:3478656/1‑101 (MQ=255)
                                                                                         gATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaa          >  1:9530563/1‑101 (MQ=255)
                                                                                          aTCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaaa         >  1:5854107/1‑101 (MQ=255)
                                                                                           tCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAg        <  1:12172103/101‑1 (MQ=255)
                                                                                            ccAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGa       <  1:1656786/101‑1 (MQ=255)
                                                                                            ccAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGa       >  1:10604445/1‑101 (MQ=255)
                                                                                            ccAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGa       <  1:10085681/101‑1 (MQ=255)
                                                                                              caggcTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGat      <  1:2107042/97‑1 (MQ=255)
                                                                                                ggcTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGatat    <  1:6812090/99‑1 (MQ=255)
                                                                                                ggcTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGatat    <  1:14106725/99‑1 (MQ=255)
                                                                                                ggcTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGatat    <  1:9354502/99‑1 (MQ=255)
                                                                                                  cTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGATAttt  <  1:14281345/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.