Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,326,067 +C 11.8% coding (724/1935 nt) ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,326,0671.→C11.8% 163.7 ‑1.0 ‑0.1280coding (724/1935 nt)ftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  new base (13/20):  ref base (119/128):  total (132/148)
Fisher's exact test strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.59e-01

Only 100 of 305 total aligned reads displayed.

AGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGC  >  NC_007779/3325967‑3326166
                                                                                                     |                                                                                                   
   gTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTttctt                                                                                                  <  1:10053652/101‑1 (MQ=255)
   gTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTttctt                                                                                                  <  1:1724609/101‑1 (MQ=255)
    tCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTttcttc                                                                                                 <  1:2323520/101‑1 (MQ=255)
     cTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCg                                                                                                <  1:8987810/101‑1 (MQ=255)
      tGTTGACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGc                                                                                               >  1:13519181/1‑101 (MQ=255)
      tGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGc                                                                                               <  1:9395591/101‑1 (MQ=255)
         tCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTg                                                                                            <  1:2651084/101‑1 (MQ=255)
         tCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTg                                                                                            <  1:9544954/101‑1 (MQ=255)
          cACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGt                                                                                           <  1:12897506/101‑1 (MQ=255)
            cGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTc                                                                                         >  1:12420038/1‑101 (MQ=255)
             gTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCg                                                                                        <  1:10651798/101‑1 (MQ=255)
              ttCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGa                                                                                       >  1:624894/1‑101 (MQ=255)
                cATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAAc                                                                                     <  1:4207353/101‑1 (MQ=255)
                cATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAAc                                                                                     <  1:12163622/101‑1 (MQ=255)
                  tCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACAt                                                                                   <  1:8089261/101‑1 (MQ=255)
                    gTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGt                                                                                 <  1:10393601/101‑1 (MQ=255)
                    gTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGt                                                                                 >  1:4046686/1‑101 (MQ=255)
                      gACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCa                                                                               <  1:14300102/101‑1 (MQ=255)
                        ccaccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACg                                                                             <  1:2898999/101‑1 (MQ=255)
                         caccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGa                                                                            >  1:4288458/1‑101 (MQ=255)
                           ccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAac                                                                          <  1:9739846/101‑1 (MQ=255)
                             gCCCAGACCAGCGCCACGCTGGCGTCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAacac                                                                        <  1:9226161/101‑1 (MQ=255)
                              cccAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACg                                                                       <  1:3982259/101‑1 (MQ=255)
                               ccAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACgg                                                                      >  1:3439922/1‑101 (MQ=255)
                                 aGACCAGCGCCACGCTGGCGGGCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGa                                                                    <  1:7596493/101‑1 (MQ=255)
                                 aGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGa                                                                    <  1:13812904/101‑1 (MQ=255)
                                 aGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGa                                                                    >  1:14021894/1‑101 (MQ=255)
                                  gACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGAt                                                                   <  1:11507879/101‑1 (MQ=255)
                                   aCCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATg                                                                  <  1:6968747/101‑1 (MQ=255)
                                    ccAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGc                                                                 >  1:6171211/1‑101 (MQ=255)
                                     cAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCa                                                                >  1:6819310/1‑101 (MQ=255)
                                      aGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCAc                                                               >  1:5958521/1‑101 (MQ=255)
                                        cgcCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCAccc                                                             <  1:516215/101‑1 (MQ=255)
                                         gcCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCca                                                            <  1:8536969/101‑1 (MQ=255)
                                           cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                                          <  1:3854978/101‑1 (MQ=255)
                                           cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                                          <  1:9901728/101‑1 (MQ=255)
                                           cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                                          <  1:11874903/101‑1 (MQ=255)
                                                 ggcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacg                                                    <  1:7259775/101‑1 (MQ=255)
                                                  gcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacga                                                   >  1:11653724/1‑101 (MQ=255)
                                                   cggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGaa                                                  >  1:10941099/1‑101 (MQ=255)
                                                    ggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAAc                                                 >  1:1884465/1‑101 (MQ=255)
                                                     gcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACa                                                <  1:14151841/101‑1 (MQ=255)
                                                      ccTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAt                                               >  1:524815/1‑101 (MQ=255)
                                                      ccTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAt                                               >  1:7792876/1‑101 (MQ=255)
                                                        tACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAttt                                             >  1:2814785/1‑101 (MQ=255)
                                                         aCGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTc                                            >  1:5697835/1‑101 (MQ=255)
                                                         aCGGCGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTc                                            <  1:6857395/101‑1 (MQ=255)
                                                          cGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCt                                           <  1:5889734/101‑1 (MQ=255)
                                                          cGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCt                                           <  1:4911894/101‑1 (MQ=255)
                                                             cGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACg                                        >  1:13156655/1‑101 (MQ=255)
                                                             cGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACg                                        >  1:9498850/1‑101 (MQ=255)
                                                              gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                       >  1:2756230/1‑101 (MQ=255)
                                                              gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                       <  1:4412515/101‑1 (MQ=255)
                                                              gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                                       >  1:8783185/1‑101 (MQ=255)
                                                               tCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGaa                                      >  1:3107486/1‑101 (MQ=255)
                                                                cGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTACTACGAAg                                     >  1:4678655/1‑101 (MQ=255)
                                                                 aatctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGt                                    >  1:7530728/8‑101 (MQ=255)
                                                                  atctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTc                                   <  1:8176941/95‑1 (MQ=255)
                                                                   tctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTcc                                  <  1:12367150/96‑2 (MQ=255)
                                                                    ttCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAg                                 >  1:13625093/1‑101 (MQ=255)
                                                                     tcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGa                                <  1:5024504/98‑1 (MQ=255)
                                                                     tcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGa                                <  1:10784059/98‑1 (MQ=255)
                                                                     tCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGa                                >  1:2009993/1‑101 (MQ=255)
                                                                     tCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGa                                >  1:6494560/1‑101 (MQ=255)
                                                                       aTCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCAACACCGACGAACATTTCTACGAAGTCAGAAc                              >  1:2507587/1‑101 (MQ=255)
                                                                        tCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAcc                             <  1:7374401/101‑1 (MQ=255)
                                                                         cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg                            >  1:3540600/1‑101 (MQ=255)
                                                                         cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg                            <  1:6599576/101‑1 (MQ=255)
                                                                          gATAAAGATGATGCACGGTCGCGCTTTCTTCGCCTGTTCGAACACCTCACCAACACGGGATGCACCCACACCGAAGAACATTTCTACGAAGTCAGAACCgg                           <  1:10518501/101‑1 (MQ=255)
                                                                           aTAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGga                          >  1:10222009/1‑101 (MQ=255)
                                                                           aTAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGga                          >  1:2984162/1‑101 (MQ=255)
                                                                            tgaaGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCgag                         <  1:6136118/99‑1 (MQ=255)
                                                                            tAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgag                         >  1:1074744/1‑101 (MQ=255)
                                                                            tAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgag                         >  1:4754170/1‑101 (MQ=255)
                                                                             aaaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgaga                        >  1:5311716/1‑101 (MQ=255)
                                                                               aGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATa                      >  1:11302741/1‑101 (MQ=255)
                                                                                gatgatGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAg                     >  1:9644037/1‑101 (MQ=255)
                                                                                gatgatGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAg                     <  1:3698094/101‑1 (MQ=255)
                                                                                  tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa                   >  1:3294840/1‑101 (MQ=255)
                                                                                  tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa                   <  1:7905063/101‑1 (MQ=255)
                                                                                  tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa                   >  1:12061652/1‑101 (MQ=255)
                                                                                  taatgCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAtggcg                   <  1:5518931/99‑5 (MQ=255)
                                                                                   gatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaa                  <  1:6912280/101‑1 (MQ=255)
                                                                                     tgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAg                <  1:2609605/101‑1 (MQ=255)
                                                                                      gCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCAACGAAGTCAGAACCGGAGATAGTAAAGa               <  1:13430555/101‑1 (MQ=255)
                                                                                         cGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACg            >  1:12683189/1‑101 (MQ=255)
                                                                                         cGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACg            >  1:6379414/1‑101 (MQ=255)
                                                                                         cGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACg            >  1:2576320/1‑101 (MQ=255)
                                                                                              ccGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACt       <  1:13276162/101‑1 (MQ=255)
                                                                                              ccGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACt       >  1:7649355/1‑101 (MQ=255)
                                                                                               cGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACtt      >  1:12247357/1‑101 (MQ=255)
                                                                                                gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt     >  1:11037458/1‑101 (MQ=255)
                                                                                                gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt     <  1:8433475/101‑1 (MQ=255)
                                                                                                gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt     >  1:11451821/1‑101 (MQ=255)
                                                                                                gCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACttt     >  1:6079173/1‑101 (MQ=255)
                                                                                                ccTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCGCCCACACCGACGAACATTTCTACGAAGGCAGAACCGGAGATAGTAAAGAACGGAACttt     <  1:1605874/100‑1 (MQ=255)
                                                                                                 cTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttc    <  1:2173597/101‑1 (MQ=255)
                                                                                                  tttCTTCGCCTGTTCGAACATGTCACGAACACGGGATGTACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcg   >  1:2107593/1‑101 (MQ=255)
                                                                                                  tttCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcg   <  1:7075957/101‑1 (MQ=255)
                                                                                                   ttcttcGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTttcgc  <  1:11121697/101‑1 (MQ=255)
                                                                                                     |                                                                                                   
AGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGC  >  NC_007779/3325967‑3326166

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