Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,449,851 G→A 13.6% A494V (GCG→GTG)  rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,449,8510G→A13.6% 94.9 ‑0.9 ‑0.0279A494V (GCG→GTG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  new base (16/22):  ref base (120/121):  total (136/143)
Fisher's exact test strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.90e-01

Only 100 of 341 total aligned reads displayed.

CAGCACCATGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGA  >  NC_007779/3449751‑3449950
                                                                                                    |                                                                                                   
cAGCACCATGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGgttcgccgttcg                                                                                                     >  1:1077881/1‑101 (MQ=255)
 aGCACCATGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGgttcgccgttcgc                                                                                                    >  1:13850501/1‑101 (MQ=255)
  gCACCATGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGgttcgccgttcgcc                                                                                                   <  1:2931822/101‑1 (MQ=255)
  gCACCATGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGgttcgccgttcgcc                                                                                                   <  1:12929598/101‑1 (MQ=255)
  gCACCATGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGgttcgccgttcgcc                                                                                                   <  1:11311169/101‑1 (MQ=255)
  gCACCATGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGgttcgccgttcacc                                                                                                   >  1:3854180/1‑101 (MQ=255)
   cACCATGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGttcgccgttcgccg                                                                                                  >  1:7525827/1‑101 (MQ=255)
    aCCATGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCgg                                                                                                 >  1:14389764/1‑101 (MQ=255)
     ccATGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCggg                                                                                                >  1:13751996/1‑101 (MQ=255)
      cATGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCgggg                                                                                               <  1:2717917/101‑1 (MQ=255)
        tGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCACCGGGGAc                                                                                             <  1:11431246/101‑1 (MQ=255)
         gCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTAAAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACa                                                                                            >  1:5954292/1‑101 (MQ=255)
           cTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAgg                                                                                          <  1:10522452/101‑1 (MQ=255)
             tCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTTAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGAt                                                                                        <  1:12056567/101‑1 (MQ=255)
              cGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGAtg                                                                                       >  1:5470218/1‑101 (MQ=255)
              cGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCACCGGGGACAGGAtg                                                                                       >  1:11899016/1‑101 (MQ=255)
              cGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCACCGGGGACAGGAtg                                                                                       <  1:736451/101‑1 (MQ=255)
               gCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGAtgt                                                                                      <  1:3276947/101‑1 (MQ=255)
                ccTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGAtgtt                                                                                     <  1:11163187/101‑1 (MQ=255)
                ccTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGAtgtt                                                                                     >  1:4967318/1‑101 (MQ=255)
                ccTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGACACGGAACGATGATTGGTTCGCCGTTCGCCGGGGCCAGGGtgtt                                                                                     >  1:4797477/1‑101 (MQ=255)
                 cTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGAtgttg                                                                                    <  1:10111049/101‑1 (MQ=255)
                 cTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGAtgttg                                                                                    >  1:4383004/1‑101 (MQ=255)
                 cTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCACCGGGGACAGGAtgttg                                                                                    >  1:4407632/1‑101 (MQ=255)
                    tGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGttgttg                                                                                 <  1:6436872/101‑1 (MQ=255)
                     ggCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTgg                                                                                >  1:13138542/1‑101 (MQ=255)
                      gCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGt                                                                               >  1:12540517/1‑101 (MQ=255)
                      gCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGt                                                                               >  1:957623/1‑101 (MQ=255)
                       cGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTa                                                                              >  1:6923246/1‑101 (MQ=255)
                       cGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTa                                                                              >  1:12142873/1‑101 (MQ=255)
                       cGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTa                                                                              <  1:10632900/101‑1 (MQ=255)
                       cGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTa                                                                              >  1:9316493/1‑101 (MQ=255)
                       cGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTa                                                                              <  1:11703425/101‑1 (MQ=255)
                        gTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAg                                                                             >  1:9676572/1‑101 (MQ=255)
                          tAACACAGTCACGGGTCATGTAGTACAGCCCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGAc                                                                           <  1:1642722/101‑1 (MQ=255)
                              acaGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGAcatca                                                                       <  1:4227315/101‑1 (MQ=255)
                                aGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACatcatc                                                                     >  1:4545417/1‑101 (MQ=255)
                                aGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACatcatc                                                                     >  1:3165621/1‑101 (MQ=255)
                                aGTCACGGGTCATGTAGTACAGACCCAGTACAACGACCTGAGACGGAACGATGATTGGTTCGCCGTTCACCGGGGACAGGATGTTGTTGGTAGACatcatc                                                                     <  1:11468066/101‑1 (MQ=255)
                                 gTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACAtcatca                                                                    <  1:9027972/101‑1 (MQ=255)
                                 gTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACAtcatca                                                                    <  1:9006426/101‑1 (MQ=255)
                                     cGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAgcgc                                                                <  1:13207726/101‑1 (MQ=255)
                                        gTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCAcg                                                             >  1:10744126/1‑101 (MQ=255)
                                         tCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCAcgc                                                            <  1:10946224/101‑1 (MQ=255)
                                         tCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCAcgc                                                            <  1:126905/101‑1 (MQ=255)
                                         tCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCACCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCAcgc                                                            <  1:12102525/101‑1 (MQ=255)
                                          cATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCAcgcg                                                           <  1:530853/101‑1 (MQ=255)
                                          cATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCAcgcg                                                           >  1:14008018/1‑101 (MQ=255)
                                          cATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCACCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCAcgcg                                                           >  1:2192804/1‑101 (MQ=255)
                                           aTGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACgcgc                                                          >  1:3191665/1‑101 (MQ=255)
                                            tGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCt                                                         <  1:896432/101‑1 (MQ=255)
                                             gtagtaCAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCtt                                                        >  1:10361293/1‑101 (MQ=255)
                                             gtagtaCAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCtt                                                        >  1:4106672/1‑101 (MQ=255)
                                                gtaCAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCa                                                     >  1:14496303/1‑101 (MQ=255)
                                                   cAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCt                                                  >  1:1263054/1‑101 (MQ=255)
                                                      aCCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTggg                                               >  1:8729115/1‑101 (MQ=255)
                                                       cccAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGc                                              <  1:6974286/101‑1 (MQ=255)
                                                       cccAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGc                                              >  1:8799800/1‑101 (MQ=255)
                                                        ccAGTACAACGTCCTTAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCt                                             <  1:10861853/101‑1 (MQ=255)
                                                        ccAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCt                                             <  1:8530858/101‑1 (MQ=255)
                                                         cAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCtt                                            >  1:8185980/1‑101 (MQ=255)
                                                          aGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTc                                           <  1:12827656/101‑1 (MQ=255)
                                                           gTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTcc                                          >  1:8251740/1‑101 (MQ=255)
                                                            tACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCa                                         <  1:3924848/101‑1 (MQ=255)
                                                            tACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCa                                         >  1:12240214/1‑101 (MQ=255)
                                                             aCAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAg                                        >  1:1358386/1‑101 (MQ=255)
                                                             aCAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCACCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAg                                        <  1:3014488/101‑1 (MQ=255)
                                                              cAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGc                                       >  1:6366453/1‑101 (MQ=255)
                                                                aCGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGt                                     >  1:8024664/1‑101 (MQ=255)
                                                                  gTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCa                                   >  1:6082566/1‑101 (MQ=255)
                                                                   tCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAg                                  <  1:6685633/101‑1 (MQ=255)
                                                                    ccTGAGACGGAACGATGATTGGTTCGCCGTTCACCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGc                                 <  1:7299304/101‑1 (MQ=255)
                                                                      tGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCgg                               >  1:8439883/1‑101 (MQ=255)
                                                                      tGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCgg                               <  1:8363653/101‑1 (MQ=255)
                                                                      tGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCgg                               >  1:7870206/1‑101 (MQ=255)
                                                                      tGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCgg                               >  1:11044113/1‑101 (MQ=255)
                                                                          aCGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACg                           <  1:9194733/101‑1 (MQ=255)
                                                                            ggAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTg                         <  1:8681030/101‑1 (MQ=255)
                                                                            ggAACGATGATTGGTTCGCCGTTCACCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTg                         >  1:3515535/1‑101 (MQ=255)
                                                                              aaCGATGATTGGTTCGCCGTTCACCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGaa                       <  1:7432990/101‑1 (MQ=255)
                                                                                  atgatTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGc                   >  1:4714477/1‑101 (MQ=255)
                                                                                   tgatTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGcc                  >  1:1047813/1‑101 (MQ=255)
                                                                                   tgatTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGcc                  >  1:12385357/1‑101 (MQ=255)
                                                                                      ttGGTTCGCCGTTCACCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATc               <  1:14285151/101‑1 (MQ=255)
                                                                                       tGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCt              <  1:1943368/101‑1 (MQ=255)
                                                                                        ggTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTg             <  1:7103860/101‑1 (MQ=255)
                                                                                           tcgccgttcgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTc          <  1:5502434/101‑1 (MQ=255)
                                                                                             gccgttcgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCAc        >  1:3718064/1‑101 (MQ=255)
                                                                                             gccgttcgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCAc        >  1:13081747/1‑101 (MQ=255)
                                                                                              ccgttcgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCAcc       <  1:561478/101‑1 (MQ=255)
                                                                                              ccgttcgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCAcc       >  1:10431697/1‑101 (MQ=255)
                                                                                               cgttcgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCa      <  1:2323277/101‑1 (MQ=255)
                                                                                               cgttcaccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCa      <  1:1826562/101‑1 (MQ=255)
                                                                                                gttcgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCAt     >  1:3948236/1‑101 (MQ=255)
                                                                                                gttcgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCAt     >  1:3653640/1‑101 (MQ=255)
                                                                                                gttcgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCAt     >  1:3315182/1‑101 (MQ=255)
                                                                                                gttcgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCAt     >  1:6809140/1‑101 (MQ=255)
                                                                                                  tcgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCg   <  1:13368721/101‑1 (MQ=255)
                                                                                                   cgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGa  <  1:1524871/101‑1 (MQ=255)
                                                                                                   cgccgGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGa  >  1:5194967/1‑101 (MQ=255)
                                                                                                    |                                                                                                   
CAGCACCATGCCTTCGCCTTTGGCGTTAACACAGTCACGGGTCATGTAGTACAGACCCAGTACAACGTCCTGAGACGGAACGATGATTGGTTCGCCGTTCGCCGGGGACAGGATGTTGTTGGTAGACATCATCAGCGCACGCGCTTCCAGCTGGGCTTCCAGCGTCAGCGGTACGTGAACAGCCATCTGGTCACCATCGA  >  NC_007779/3449751‑3449950

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.