Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,750,122 Δ1 bp 94.8% coding (365/1416 nt) tnaA ← tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,750,1220C→.94.8% 67.9 ‑1.3 ‑2.3290coding (365/1416 nt)tnaAtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  new base (116/159):  ref base (8/6):  total (125/165)
Fisher's exact test strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.28e-03

Only 100 of 298 total aligned reads displayed.

GATATAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGA  >  NC_007779/3750022‑3750218
                                                                                                    |                                                                                                
 atatAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGct                                                                                                 <  1:11256899/101‑2 (MQ=255)
   ataGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTtt                                                                                               >  1:4958748/1‑101 (MQ=255)
    taGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTttt                                                                                              <  1:8242624/101‑1 (MQ=255)
     aGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTtttt                                                                                             <  1:10229509/101‑1 (MQ=255)
     aGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTtttt                                                                                             <  1:2850018/101‑1 (MQ=255)
       aCGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTcc                                                                                          <  1:9452685/101‑1 (MQ=255)
        cGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCt                                                                                         >  1:3397613/1‑101 (MQ=255)
         gTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTg                                                                                        <  1:33023/101‑1 (MQ=255)
          ttACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGc                                                                                       <  1:4273118/101‑1 (MQ=255)
          ttACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGc                                                                                       <  1:2145156/101‑1 (MQ=255)
           tACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCt                                                                                      >  1:9014462/1‑101 (MQ=255)
           tACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCt                                                                                      >  1:8121837/1‑101 (MQ=255)
            aCGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTc                                                                                     >  1:7227482/1‑101 (MQ=255)
             cGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcg                                                                                    <  1:2235695/101‑1 (MQ=255)
             cGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcg                                                                                    <  1:75824/101‑1 (MQ=255)
              gCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcgc                                                                                   <  1:11301508/101‑1 (MQ=255)
               cACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcgcg                                                                                  >  1:11819059/1‑101 (MQ=255)
               cACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcgcg                                                                                  >  1:570776/1‑101 (MQ=255)
                aCGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGt                                                                                 >  1:1235681/1‑101 (MQ=255)
                 cGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtt                                                                                <  1:2637941/101‑1 (MQ=255)
                 cGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtt                                                                                >  1:4815763/1‑101 (MQ=255)
                     aCAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtttttt                                                                            <  1:6812834/101‑1 (MQ=255)
                     aCAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtttttt                                                                            >  1:706295/1‑101 (MQ=255)
                     aCAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtttttt                                                                            >  1:9804503/1‑101 (MQ=255)
                      cAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTa                                                                           <  1:6776030/101‑1 (MQ=255)
                       aGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTaa                                                                          <  1:5421996/101‑1 (MQ=255)
                       aGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTaa                                                                          >  1:11386892/1‑101 (MQ=255)
                        gCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCCGGCTTTTTCCTGCTCGCGTTTTTTAAt                                                                         >  1:2107328/1‑101 (MQ=255)
                        gCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAAt                                                                         >  1:9695906/1‑101 (MQ=255)
                          cGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCa                                                                       >  1:10110074/1‑101 (MQ=255)
                          cGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCa                                                                       >  1:2778939/1‑101 (MQ=255)
                           gTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAg                                                                      >  1:9645157/1‑101 (MQ=255)
                           gTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAg                                                                      >  1:637283/1‑101 (MQ=255)
                            ttGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGt                                                                     >  1:6255595/1‑101 (MQ=255)
                             tGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTa                                                                    <  1:11187800/101‑1 (MQ=255)
                             tGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTa                                                                    <  1:3244194/101‑1 (MQ=255)
                             tGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTa                                                                    <  1:1417206/101‑1 (MQ=255)
                               aTCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTAcc                                                                  >  1:2716555/1‑101 (MQ=255)
                               aTCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTAcc                                                                  <  1:4501222/101‑1 (MQ=255)
                                tCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCg                                                                 >  1:8814350/1‑101 (MQ=255)
                                tCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCg                                                                 <  1:10170330/101‑1 (MQ=255)
                                   ggCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGaa                                                              <  1:10468035/101‑1 (MQ=255)
                                   ggCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGaa                                                              <  1:2080613/101‑1 (MQ=255)
                                      tATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAatat                                                           <  1:11721307/101‑1 (MQ=255)
                                           cccTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGAtt                                                      <  1:8182256/101‑1 (MQ=255)
                                           cccTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGAtt                                                      <  1:6181061/101‑1 (MQ=255)
                                            ccTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGAttt                                                     <  1:3177060/101‑1 (MQ=255)
                                              tGCGTGGTATCAAAGAAATAGTTAGAGAACGCCTCCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGc                                                   <  1:11500928/101‑1 (MQ=255)
                                                 gTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctct                                                >  1:5325801/1‑101 (MQ=255)
                                                 gTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctct                                                <  1:8608931/101‑1 (MQ=255)
                                                 gTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctct                                                >  1:7839716/1‑101 (MQ=255)
                                                 gTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctct                                                <  1:11994945/101‑1 (MQ=255)
                                                 gTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctct                                                <  1:10941080/101‑1 (MQ=255)
                                                  tGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTg                                               <  1:9340806/101‑1 (MQ=255)
                                                   ggTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTgc                                              <  1:2526316/101‑1 (MQ=255)
                                                   ggTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTgc                                              >  1:5468374/1‑101 (MQ=255)
                                                    gTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTgcg                                             >  1:1898851/1‑101 (MQ=255)
                                                      aTCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGcc                                           <  1:10821307/101‑1 (MQ=255)
                                                      aTCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGcc                                           <  1:10338455/101‑1 (MQ=255)
                                                       tCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCa                                          >  1:4043269/1‑101 (MQ=255)
                                                        cAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCAc                                         >  1:5542078/1‑101 (MQ=255)
                                                        cAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCAc                                         >  1:9204221/1‑101 (MQ=255)
                                                         aaaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAGTATAGATTTGCTCTGCGCCACg                                        <  1:3040589/101‑1 (MQ=255)
                                                         aaaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACg                                        >  1:6965198/1‑101 (MQ=255)
                                                          aaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACgg                                       <  1:1644265/101‑1 (MQ=255)
                                                          aaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATCTGCTCTGCGCCACgg                                       >  1:6626024/1‑101 (MQ=255)
                                                           aGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGc                                      >  1:387262/1‑101 (MQ=255)
                                                              aaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCt                                   <  1:5623244/101‑1 (MQ=255)
                                                              aaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCt                                   <  1:5065654/101‑1 (MQ=255)
                                                                 aGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGt                                <  1:3893192/101‑1 (MQ=255)
                                                                 aGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGt                                >  1:1767758/1‑101 (MQ=255)
                                                                  gTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTg                               >  1:8633238/1‑101 (MQ=255)
                                                                  gTTAGAGAACGCCACCATTTTGCCGCGATCCAGGCTGTTACCTGCTCGCGTGTTTTAATCAGGACCGGACTATAGATTTGCTCTGCGCCACGGCCCtagcg                               >  1:12199823/1‑97 (MQ=255)
                                                                   ttAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGa                              >  1:9472527/1‑101 (MQ=255)
                                                                    tAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAg                             <  1:8912418/101‑1 (MQ=255)
                                                                    tAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAg                             <  1:5804158/101‑1 (MQ=255)
                                                                     agagAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGt                            <  1:1343667/101‑1 (MQ=255)
                                                                     agagAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTc                           >  1:3942804/1‑101 (MQ=255)
                                                                      gagaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTc                           <  1:668843/101‑1 (MQ=255)
                                                                      gagaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTc                           <  1:3606188/101‑1 (MQ=255)
                                                                       agaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCg                          >  1:2415434/1‑101 (MQ=255)
                                                                       agaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCg                          <  1:8457165/101‑1 (MQ=255)
                                                                        gaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCgg                         <  1:1989191/101‑1 (MQ=255)
                                                                         aaCGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGa                        >  1:11005506/1‑101 (MQ=255)
                                                                          aCGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGaa                       <  1:4419461/101‑1 (MQ=255)
                                                                           cGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAAt                      >  1:10050391/1‑101 (MQ=255)
                                                                            gCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATg                     <  1:4600743/101‑1 (MQ=255)
                                                                              caccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgt                   >  1:9757224/1‑101 (MQ=255)
                                                                               accaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgtg                  <  1:7802013/101‑1 (MQ=255)
                                                                                 caTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTa                >  1:7178564/1‑101 (MQ=255)
                                                                                  aTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAt               >  1:5966734/1‑101 (MQ=255)
                                                                                   ttttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAtt              <  1:7944352/101‑1 (MQ=255)
                                                                                    tttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTg             <  1:8364551/101‑1 (MQ=255)
                                                                                     ttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGa            <  1:3513267/101‑1 (MQ=255)
                                                                                        ctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATaa         <  1:8536482/101‑1 (MQ=255)
                                                                                           cGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCa      <  1:2472082/101‑1 (MQ=255)
                                                                                             aTCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaaa    <  1:9079371/101‑1 (MQ=255)
                                                                                              tCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAg   >  1:4667684/1‑101 (MQ=255)
                                                                                               ccAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGa  <  1:5579938/101‑1 (MQ=255)
                                                                                               ccAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGa  >  1:955026/1‑101 (MQ=255)
                                                                                                    |                                                                                                
GATATAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGA  >  NC_007779/3750022‑3750218

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