Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,750,122 Δ1 bp 95.4% 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→.95.4% 60.2 ‑2.4 ‑1.9285coding (365/1416 nt)tnaAtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  new base (134/138):  ref base (3/10):  total (137/148)
Fisher's exact test strand distribution p-value = 8.80e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.18e-02

Only 100 of 294 total aligned reads displayed.

GATATAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGA  >  NC_007779/3750022‑3750218
                                                                                                    |                                                                                                
gATATAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGc                                                                                                  >  1:2898066/1‑101 (MQ=255)
gATATAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGc                                                                                                  <  1:2402725/101‑1 (MQ=255)
   ataGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTtt                                                                                               <  1:8023633/101‑1 (MQ=255)
     aGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTtttt                                                                                             >  1:10039392/1‑101 (MQ=255)
     aGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTtttt                                                                                             <  1:9311384/101‑1 (MQ=255)
      gACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTttttc                                                                                            <  1:2539453/101‑2 (MQ=255)
       aCGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTcc                                                                                          <  1:10576727/101‑1 (MQ=255)
         gTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTg                                                                                        <  1:1153045/101‑1 (MQ=255)
         gTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTg                                                                                        <  1:532513/101‑1 (MQ=255)
               cACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcgcg                                                                                  <  1:2207726/101‑1 (MQ=255)
               cACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTcgcg                                                                                  >  1:4950231/1‑101 (MQ=255)
                aCGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGt                                                                                 >  1:7251704/1‑101 (MQ=255)
                aCGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGt                                                                                 <  1:4429184/101‑1 (MQ=255)
                  ggTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGttt                                                                               <  1:3001242/101‑1 (MQ=255)
                  ggTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGttt                                                                               >  1:11965893/1‑101 (MQ=255)
                  ggTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTGTTTCCTGCGCGCGttt                                                                               >  1:7843894/1‑101 (MQ=255)
                   gTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtttt                                                                              <  1:1411502/101‑1 (MQ=255)
                    tACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGttttt                                                                             <  1:9915200/101‑1 (MQ=255)
                     aCAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtttttt                                                                            >  1:1489791/1‑101 (MQ=255)
                     aCAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGtttttt                                                                            <  1:11591298/101‑1 (MQ=255)
                      cAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTa                                                                           >  1:3118388/1‑101 (MQ=255)
                      cAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTa                                                                           <  1:8065191/101‑1 (MQ=255)
                      cAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTa                                                                           >  1:1929342/1‑101 (MQ=255)
                           gTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAg                                                                      >  1:5123236/1‑101 (MQ=255)
                           gTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAg                                                                      >  1:6533452/1‑101 (MQ=255)
                             tGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTa                                                                    >  1:3700566/1‑101 (MQ=255)
                             tGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTa                                                                    <  1:6129897/101‑1 (MQ=255)
                                tCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCg                                                                 <  1:4403555/101‑1 (MQ=255)
                                 cTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCgg                                                                <  1:7276052/101‑1 (MQ=255)
                                 cTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCgg                                                                <  1:7223372/101‑1 (MQ=255)
                                 cTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCgg                                                                <  1:2701193/101‑1 (MQ=255)
                                  tGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGa                                                               >  1:896565/1‑101 (MQ=255)
                                    gCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAat                                                             <  1:7507546/101‑1 (MQ=255)
                                     cTATGGCCCTGCGTTGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAata                                                            <  1:1284593/101‑1 (MQ=255)
                                      tATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAatat                                                           >  1:10175346/1‑101 (MQ=255)
                                      tATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAatat                                                           >  1:4572088/1‑101 (MQ=255)
                                      tATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAatat                                                           >  1:10349072/1‑101 (MQ=255)
                                       aTGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAAtata                                                          <  1:11206659/101‑1 (MQ=255)
                                         ggCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGa                                                        <  1:4251344/101‑1 (MQ=255)
                                          gCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGAt                                                       >  1:10679986/1‑101 (MQ=255)
                                          gCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGAt                                                       >  1:4835840/1‑101 (MQ=255)
                                          gCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGAt                                                       <  1:582188/101‑1 (MQ=255)
                                           cccTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGAtt                                                      <  1:11449953/101‑1 (MQ=255)
                                            ccTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGAttt                                                     <  1:3817426/101‑1 (MQ=255)
                                             cTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTg                                                    >  1:8324741/1‑101 (MQ=255)
                                             cTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTg                                                    <  1:2443277/101‑1 (MQ=255)
                                               gCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGcn                                                  <  1:4370059/101‑2 (MQ=255)
                                                cGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctc                                                 <  1:785690/101‑1 (MQ=255)
                                                cGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGctc                                                 >  1:4847770/1‑101 (MQ=255)
                                                    gTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTgcg                                             <  1:4663862/101‑1 (MQ=255)
                                                      aTCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGcc                                           <  1:1781209/101‑1 (MQ=255)
                                                      aTCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGcc                                           >  1:2780025/1‑101 (MQ=255)
                                                         aaaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACg                                        >  1:11253687/1‑101 (MQ=255)
                                                          aaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACgg                                       <  1:3368708/101‑1 (MQ=255)
                                                          aaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACgg                                       <  1:5139638/101‑1 (MQ=255)
                                                          aaGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACgg                                       >  1:10926903/1‑101 (MQ=255)
                                                           aGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGc                                      >  1:10087832/1‑101 (MQ=255)
                                                            gAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGcc                                     >  1:6285429/1‑101 (MQ=255)
                                                             aaaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGccc                                    >  1:12126575/1‑101 (MQ=255)
                                                              aaTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCt                                   >  1:9071042/1‑101 (MQ=255)
                                                               aTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTg                                  >  1:2669272/1‑101 (MQ=255)
                                                               aTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTg                                  >  1:6481678/1‑101 (MQ=255)
                                                               aTAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTg                                  <  1:5916220/101‑1 (MQ=255)
                                                                tAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTgg                                 <  1:9693952/101‑1 (MQ=255)
                                                                 aGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGt                                >  1:11652298/1‑101 (MQ=255)
                                                                 aGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGt                                >  1:2066209/1‑101 (MQ=255)
                                                                 aGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGt                                <  1:8291857/101‑1 (MQ=255)
                                                                   ttAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGa                              >  1:3236400/1‑101 (MQ=255)
                                                                    tAGAGAACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAg                             >  1:329887/1‑101 (MQ=255)
                                                                      gagaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTc                           >  1:2182493/1‑101 (MQ=255)
                                                                      gagaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTc                           <  1:6031997/101‑1 (MQ=255)
                                                                      gagaACGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTc                           <  1:10773738/101‑1 (MQ=255)
                                                                         aaCGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGa                        <  1:281327/101‑1 (MQ=255)
                                                                          aCGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGaa                       <  1:7061810/101‑1 (MQ=255)
                                                                           cGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAAt                      <  1:5535628/101‑1 (MQ=255)
                                                                           cGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAAt                      <  1:7621932/101‑1 (MQ=255)
                                                                           cGCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAAt                      <  1:9507562/101‑1 (MQ=255)
                                                                            gCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATg                     >  1:123980/1‑101 (MQ=255)
                                                                            gCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATg                     <  1:39482/101‑1 (MQ=255)
                                                                            gCCACCATTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATg                     <  1:350059/101‑1 (MQ=255)
                                                                             ccaccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATgg                    >  1:6720700/1‑101 (MQ=255)
                                                                              caccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgt                   >  1:8581246/1‑101 (MQ=255)
                                                                              caccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgt                   >  1:4019994/1‑101 (MQ=255)
                                                                              caccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGCCGGAATGgt                   >  1:6889203/1‑101 (MQ=255)
                                                                               accaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgtg                  <  1:586376/101‑1 (MQ=255)
                                                                                ccaTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGgtgt                 <  1:1542128/101‑1 (MQ=255)
                                                                                  aTTTTGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAt               <  1:8168705/101‑1 (MQ=255)
                                                                                   ttttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTAtt              <  1:5679936/101‑1 (MQ=255)
                                                                                    tttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTg             <  1:2288090/101‑1 (MQ=255)
                                                                                     ttGCTGCGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGa            <  1:9469124/101‑1 (MQ=255)
                                                                                      tgctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGAt           <  1:4068519/101‑1 (MQ=255)
                                                                                        ctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGGCCTGGTGAGTCGGAATGGTGTATTGATaa         >  1:8473644/1‑101 (MQ=255)
                                                                                        ctgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATaa         <  1:4712757/101‑1 (MQ=255)
                                                                                         tgcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATTGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAAc        >  1:10005890/1‑101 (MQ=255)
                                                                                          gcGATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAAcc       >  1:8787510/1‑101 (MQ=255)
                                                                                            gATCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaa     <  1:11162206/101‑1 (MQ=255)
                                                                                             aTCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaaa    <  1:7431266/101‑1 (MQ=255)
                                                                                             aTCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaaa    <  1:1521086/101‑1 (MQ=255)
                                                                                             aTCCAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCaaa    <  1:11889914/101‑1 (MQ=255)
                                                                                               ccAGGCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGa  <  1:1591714/101‑1 (MQ=255)
                                                                                                    |                                                                                                
GATATAGACGTTACGCACGGTACAGCCGTTGATCTGGCTATGGCCCTGCGTGGTATCAAAGAAATAGTTAGAGAACGCCACCATTTTGCTGCGATCCAGGCCTTTTTCCTGCTCGCGTTTTTTAATCAGTACCGGAATATAGATTTGCTCTGCGCCACGGCCCTGGTGAGTCGGAATGGTGTATTGATAACCAAAGA  >  NC_007779/3750022‑3750218

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.