Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,750,122 Δ1 bp 96.9% 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→.96.9% 33.3 ‑0.7 ‑1.6254coding (365/1416 nt)tnaAtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  new base (109/137):  ref base (5/3):  total (114/140)
Fisher's exact test strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.35e-02

Only 100 of 261 total aligned reads displayed.

GATATAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAG  >  NC_007779/3750022‑3750217
                                                                                                    |                                                                                               
 atatAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGct                                                                                                >  1:948153/1‑100 (MQ=255)
  tataGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGctt                                                                                               <  1:6714597/101‑3 (MQ=255)
    taGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTttt                                                                                             <  1:8964112/101‑1 (MQ=255)
      gACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTttttc                                                                                           >  1:9768866/1‑100 (MQ=255)
       aCGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTcc                                                                                         >  1:8626925/1‑101 (MQ=255)
        cGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCt                                                                                        >  1:4237882/1‑101 (MQ=255)
        cGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCt                                                                                        <  1:12561369/101‑1 (MQ=255)
         gTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTg                                                                                       >  1:5641914/1‑101 (MQ=255)
         gTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTg                                                                                       >  1:5866116/1‑101 (MQ=255)
          ttACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGc                                                                                      <  1:6163136/101‑1 (MQ=255)
          ttACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGc                                                                                      <  1:10448820/101‑1 (MQ=255)
            aCGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTc                                                                                    <  1:11184745/101‑1 (MQ=255)
            aCGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTc                                                                                    <  1:12603376/101‑1 (MQ=255)
               cACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcgcg                                                                                 <  1:8922034/101‑1 (MQ=255)
               cACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcgcg                                                                                 <  1:7078591/101‑1 (MQ=255)
                aCGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGt                                                                                >  1:9345874/1‑101 (MQ=255)
                aCGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGt                                                                                >  1:4493970/1‑101 (MQ=255)
                aCGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGt                                                                                <  1:8815460/101‑1 (MQ=255)
                 cGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtt                                                                               <  1:325870/101‑1 (MQ=255)
                    tACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGttttt                                                                            >  1:6654840/1‑101 (MQ=255)
                      cggCCGTTGATCTGGTTATGGCCTTGGGTGGTATAAAAAAAATAGTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTaa                                                                         <  1:9667110/99‑1 (MQ=255)
                      cAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTa                                                                          >  1:1651758/1‑101 (MQ=255)
                       aGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTaa                                                                         >  1:6864577/1‑101 (MQ=255)
                       aGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTaa                                                                         >  1:7727682/1‑101 (MQ=255)
                       aGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTaa                                                                         <  1:453098/101‑1 (MQ=255)
                         ccGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATc                                                                       <  1:11994373/101‑1 (MQ=255)
                           gTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAg                                                                     <  1:8464681/101‑1 (MQ=255)
                           gTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAg                                                                     <  1:790105/101‑1 (MQ=255)
                               aTCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTAcc                                                                 <  1:5418250/101‑1 (MQ=255)
                                  tGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGa                                                              <  1:9108686/101‑1 (MQ=255)
                                    gCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAat                                                            <  1:875837/101‑1 (MQ=255)
                                     cTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAata                                                           <  1:6266625/101‑1 (MQ=255)
                                       aTGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAAtata                                                         <  1:8531186/101‑1 (MQ=255)
                                        tGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAg                                                        <  1:1295928/101‑1 (MQ=255)
                                        tGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAg                                                        <  1:6122104/101‑1 (MQ=255)
                                             cTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTg                                                   <  1:6949214/101‑1 (MQ=255)
                                             cTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTg                                                   <  1:12841490/101‑1 (MQ=255)
                                             cTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGAGTTTTTTAATCAGTACCGGAATATAGATTTg                                                   <  1:2461959/101‑1 (MQ=255)
                                               gCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGct                                                 >  1:11598542/1‑101 (MQ=255)
                                                 gTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctct                                               >  1:12025506/1‑101 (MQ=255)
                                                 gTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctct                                               <  1:4221575/101‑1 (MQ=255)
                                                  tGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTg                                              >  1:3391805/1‑101 (MQ=255)
                                                     tATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTgcgc                                           <  1:5503262/101‑1 (MQ=255)
                                                     tATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTgcgc                                           <  1:2160677/101‑1 (MQ=255)
                                                      aTCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGcc                                          >  1:7523482/1‑101 (MQ=255)
                                                      aTCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGcc                                          <  1:9417533/101‑1 (MQ=255)
                                                       tCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCa                                         >  1:2023013/1‑101 (MQ=255)
                                                        cAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCAc                                        >  1:2600078/1‑101 (MQ=255)
                                                         aaaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACg                                       >  1:7353199/1‑101 (MQ=255)
                                                          aaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACgg                                      >  1:4084979/1‑101 (MQ=255)
                                                          aaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACgg                                      <  1:3814028/101‑1 (MQ=255)
                                                           aGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGc                                     >  1:12425944/1‑101 (MQ=255)
                                                           aGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGc                                     >  1:3177152/1‑101 (MQ=255)
                                                            gAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGcc                                    >  1:12797991/1‑101 (MQ=255)
                                                              aaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCt                                  <  1:12145787/101‑1 (MQ=255)
                                                                  gTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTg                              >  1:3892017/1‑101 (MQ=255)
                                                                  gTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTg                              >  1:8194070/1‑101 (MQ=255)
                                                                   ttAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGa                             <  1:6605776/101‑1 (MQ=255)
                                                                   ttAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGa                             <  1:1735480/101‑1 (MQ=255)
                                                                   ttAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGa                             >  1:5286290/1‑101 (MQ=255)
                                                                     agagAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGt                           <  1:5076628/101‑1 (MQ=255)
                                                                     agagAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTACTCTGCGCCACGGCCCTGGTGAGt                           <  1:6298580/101‑1 (MQ=255)
                                                                      gagaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTc                          >  1:11306915/1‑101 (MQ=255)
                                                                       agaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCg                         >  1:9384279/1‑101 (MQ=255)
                                                                       agaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCg                         >  1:9518672/1‑101 (MQ=255)
                                                                       agaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCg                         <  1:589256/101‑1 (MQ=255)
                                                                        gaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCgg                        >  1:10161584/1‑101 (MQ=255)
                                                                        gaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCgg                        >  1:11728675/1‑101 (MQ=255)
                                                                        gaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCgg                        <  1:681857/101‑1 (MQ=255)
                                                                         aaCGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGa                       <  1:10126973/101‑1 (MQ=255)
                                                                         aaCGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGa                       >  1:5541935/1‑101 (MQ=255)
                                                                          aCGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGaa                      <  1:898705/101‑1 (MQ=255)
                                                                          aCGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGaa                      <  1:2352780/101‑1 (MQ=255)
                                                                             ccaccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATgg                   <  1:10185881/101‑1 (MQ=255)
                                                                                ccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgtgt                <  1:9714969/101‑1 (MQ=255)
                                                                                ccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgtgt                <  1:8384875/101‑1 (MQ=255)
                                                                                ccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgtgt                <  1:5186326/101‑1 (MQ=255)
                                                                                 caTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTa               >  1:7975705/1‑101 (MQ=255)
                                                                                 caTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTa               <  1:8242107/101‑1 (MQ=255)
                                                                                 caTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTa               >  1:3600452/1‑101 (MQ=255)
                                                                                  aTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAt              >  1:9148433/1‑101 (MQ=255)
                                                                                  aTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAt              <  1:6361413/101‑1 (MQ=255)
                                                                                  aTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAt              <  1:5830864/101‑1 (MQ=255)
                                                                                  aTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGGGTAt              >  1:9845056/1‑101 (MQ=255)
                                                                                   ttttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAtt             >  1:11420189/1‑101 (MQ=255)
                                                                                   ttttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAtt             >  1:6516603/1‑101 (MQ=255)
                                                                                   ttttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAtt             >  1:6011381/1‑101 (MQ=255)
                                                                                   ttttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAtt             <  1:1811250/101‑1 (MQ=255)
                                                                                    tttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTg            >  1:552274/1‑101 (MQ=255)
                                                                                    tttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTg            >  1:2515786/1‑101 (MQ=255)
                                                                                    tttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTg            <  1:10023427/101‑1 (MQ=255)
                                                                                     ttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGa           <  1:10684640/101‑1 (MQ=255)
                                                                                       gctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATa         <  1:1424043/101‑1 (MQ=255)
                                                                                         tgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAAc       >  1:3020357/1‑101 (MQ=255)
                                                                                           cGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCa     <  1:7075195/101‑1 (MQ=255)
                                                                                           cGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCa     >  1:10866986/1‑101 (MQ=255)
                                                                                            gATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaa    >  1:8038581/1‑101 (MQ=255)
                                                                                            gATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaa    >  1:7147665/1‑101 (MQ=255)
                                                                                            gATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaa    <  1:12240188/101‑1 (MQ=255)
                                                                                              tCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAg  >  1:2271809/1‑101 (MQ=255)
                                                                                                    |                                                                                               
GATATAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAG  >  NC_007779/3750022‑3750217

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