Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,326,050 G→A 12.5% I247I (ATC→ATT ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,326,0500G→A12.5% 78.9 ‑0.3 ‑2.4256I247I (ATC→ATTftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  new base (13/19):  ref base (102/122):  total (115/141)
Fisher's exact test strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.22e-03

Only 100 of 283 total aligned reads displayed.

CAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAA  >  NC_007779/3325950‑3326150
                                                                                                    |                                                                                                     
 aaCCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAgatga                                                                                                      >  1:13305906/1‑101 (MQ=255)
  aCCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAgatgat                                                                                                     <  1:1478906/101‑1 (MQ=255)
    cAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGc                                                                                                   >  1:8907600/1‑101 (MQ=255)
         tCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGt                                                                                              <  1:2873704/101‑1 (MQ=255)
         tCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGt                                                                                              <  1:13370537/101‑1 (MQ=255)
         tCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGt                                                                                              <  1:7838951/101‑1 (MQ=255)
          cTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTg                                                                                             <  1:3280396/101‑1 (MQ=255)
           tGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGc                                                                                            >  1:10023522/1‑101 (MQ=255)
              ttCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGc                                                                                         <  1:4260379/101‑1 (MQ=255)
              ttCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGc                                                                                         >  1:6758434/1‑101 (MQ=255)
              ttCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGc                                                                                         >  1:13190748/1‑101 (MQ=255)
               tCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCt                                                                                        <  1:9858502/101‑1 (MQ=255)
                cAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCtt                                                                                       <  1:7427594/101‑1 (MQ=255)
                    gTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTttctt                                                                                  <  1:5341405/101‑1 (MQ=255)
                    gTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTttctt                                                                                  <  1:1724609/101‑1 (MQ=255)
                    gTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTttctt                                                                                  <  1:10053652/101‑1 (MQ=255)
                    gTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTttctt                                                                                  >  1:7535708/1‑101 (MQ=255)
                      cTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCg                                                                                <  1:8987810/101‑1 (MQ=255)
                       tGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGc                                                                               <  1:9395591/101‑1 (MQ=255)
                          tCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTg                                                                            <  1:2651084/101‑1 (MQ=255)
                          tCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTg                                                                            >  1:1609514/1‑101 (MQ=255)
                          tCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTg                                                                            >  1:4888607/1‑101 (MQ=255)
                          tCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTg                                                                            <  1:9544954/101‑1 (MQ=255)
                            aCGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGtt                                                                          >  1:5921767/1‑101 (MQ=255)
                            aCGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGtt                                                                          <  1:2078596/101‑1 (MQ=255)
                             cGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTc                                                                         >  1:12420038/1‑101 (MQ=255)
                               ttCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGa                                                                       >  1:7788490/1‑101 (MQ=255)
                                tCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGaa                                                                      <  1:5796223/101‑1 (MQ=255)
                                 cATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAAc                                                                     >  1:10966708/1‑101 (MQ=255)
                                 cATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAAc                                                                     <  1:12163622/101‑1 (MQ=255)
                                  aTCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACa                                                                    <  1:4934907/101‑1 (MQ=255)
                                   tCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACAt                                                                   <  1:8089261/101‑1 (MQ=255)
                                     gTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGt                                                                 >  1:4046686/1‑101 (MQ=255)
                                          caccGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGa                                                            >  1:4288458/1‑101 (MQ=255)
                                           accGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGaa                                                           >  1:2805352/1‑101 (MQ=255)
                                             cGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAaca                                                         <  1:6950007/101‑1 (MQ=255)
                                              gCCCAGACCAGCGCCACGCTGGCGTCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAacac                                                        <  1:9226161/101‑1 (MQ=255)
                                               cccAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACg                                                       >  1:14199231/1‑101 (MQ=255)
                                               cccAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACg                                                       <  1:3982259/101‑1 (MQ=255)
                                                ccAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACgg                                                      >  1:10530136/1‑101 (MQ=255)
                                                  aGACCAGCGCCACGCTGGCGGGCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGa                                                    <  1:7596493/101‑1 (MQ=255)
                                                  aGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGa                                                    <  1:13812904/101‑1 (MQ=255)
                                                   gACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGAt                                                   <  1:11507879/101‑1 (MQ=255)
                                                     ccAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGc                                                 >  1:6171211/1‑101 (MQ=255)
                                                      cAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCAACGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCa                                                >  1:13899819/1‑101 (MQ=255)
                                                         cgcCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCAccc                                             <  1:516215/101‑1 (MQ=255)
                                                          gcCACGCTGGCGGCCTCCGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCCCCca                                            <  1:4467038/101‑1 (MQ=255)
                                                          gcCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCca                                            <  1:8536969/101‑1 (MQ=255)
                                                           ccACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcac                                           <  1:5112210/101‑1 (MQ=255)
                                                           ccACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcac                                           <  1:2384702/101‑1 (MQ=255)
                                                            cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                          <  1:9901728/101‑1 (MQ=255)
                                                            cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                          <  1:11874903/101‑1 (MQ=255)
                                                            cACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCcaca                                          <  1:3854978/101‑1 (MQ=255)
                                                             aCGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCacac                                         <  1:1161842/101‑1 (MQ=255)
                                                               gCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcg                                       <  1:3437955/101‑1 (MQ=255)
                                                               gCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcg                                       <  1:11086068/101‑1 (MQ=255)
                                                                cTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcga                                      >  1:4807811/1‑101 (MQ=255)
                                                                cTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcga                                      <  1:10667842/101‑1 (MQ=255)
                                                                  ggcggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACcgacg                                    <  1:11358308/101‑1 (MQ=255)
                                                                    cggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGaa                                  >  1:1945928/1‑101 (MQ=255)
                                                                     ggcCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAAc                                 <  1:6015184/101‑1 (MQ=255)
                                                                       ccTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAt                               <  1:2567303/101‑1 (MQ=255)
                                                                       ccTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAt                               >  1:524815/1‑101 (MQ=255)
                                                                        cTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAtt                              >  1:7351312/1‑101 (MQ=255)
                                                                         tACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACAttt                             <  1:3570722/101‑1 (MQ=255)
                                                                              cGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACg                        >  1:9498850/1‑101 (MQ=255)
                                                                              cGTCAATCTCGTCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACg                        >  1:4107251/1‑101 (MQ=255)
                                                                               gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                       >  1:2207098/1‑101 (MQ=255)
                                                                               gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                       >  1:2756230/1‑101 (MQ=255)
                                                                               gTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGa                       <  1:12660797/101‑1 (MQ=255)
                                                                                tcaatctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGaa                      <  1:8453802/92‑1 (MQ=255)
                                                                                tCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGaa                      >  1:3107486/1‑101 (MQ=255)
                                                                                tCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGaa                      >  1:2605634/1‑101 (MQ=255)
                                                                                  aatctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGt                    >  1:10414826/8‑101 (MQ=255)
                                                                                   atctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTc                   <  1:8176941/95‑1 (MQ=255)
                                                                                    tttCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCa                  <  1:14150035/101‑1 (MQ=255)
                                                                                     ttCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAg                 >  1:13625093/1‑101 (MQ=255)
                                                                                     ctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTccg                 >  1:10809357/5‑99 (MQ=255)
                                                                                     ctcgtCGATGAAGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTccg                 <  1:11677713/97‑3 (MQ=255)
                                                                                      tCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGa                >  1:6494560/1‑101 (MQ=255)
                                                                                         tCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAAcc             <  1:7066761/101‑1 (MQ=255)
                                                                                          cGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCg            <  1:13494102/101‑1 (MQ=255)
                                                                                           gATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCgg           <  1:12915467/101‑1 (MQ=255)
                                                                                            aTAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGga          >  1:10222009/1‑101 (MQ=255)
                                                                                             tgaaGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCgag         <  1:6136118/99‑1 (MQ=255)
                                                                                             tAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgag         >  1:1074744/1‑101 (MQ=255)
                                                                                             tAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgag         >  1:2949972/1‑101 (MQ=255)
                                                                                              aaaGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGgaga        <  1:1218520/101‑1 (MQ=255)
                                                                                               aaGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAt       <  1:6340959/101‑1 (MQ=255)
                                                                                               aaGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAt       <  1:5608202/101‑1 (MQ=255)
                                                                                                aGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATa      >  1:11302741/1‑101 (MQ=255)
                                                                                                aGATAATGCACGGGCTGTTTTCCTTGGCCTGATTAAACAAGTCACGGACACGAGAGGCACCGACACCGACGAACATTTCGACGAAGTCCGAACCCGAGAtg      <  1:6566409/101‑2 (MQ=255)
                                                                                                 gatgatGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAg     >  1:9644037/1‑101 (MQ=255)
                                                                                                 gatgatGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAg     <  1:972200/101‑1 (MQ=255)
                                                                                                 gatgatGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAg     >  1:3474352/1‑101 (MQ=255)
                                                                                                  atgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGt    <  1:4737395/101‑1 (MQ=255)
                                                                                                   tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa   >  1:12061652/1‑101 (MQ=255)
                                                                                                   tgatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTa   >  1:6847116/1‑101 (MQ=255)
                                                                                                    gatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaa  >  1:723459/1‑101 (MQ=255)
                                                                                                    gatgCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTaa  <  1:6912280/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
CAACCAGCATCTGGTTCAGAGTCTGTTCACGTTCATCGTGACCACCGCCCAGACCAGCGCCACGCTGGCGGCCTACGGCGTCGATTTCATCGATAAAGATGATGCACGGTGCCGCTTTCTTCGCCTGTTCGAACATGTCACGAACACGGGATGCACCCACACCGACGAACATTTCTACGAAGTCAGAACCGGAGATAGTAA  >  NC_007779/3325950‑3326150

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