Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,731,822 C→G 21.0% V336V (GTG→GTC clpB ← protein disaggregation chaperone

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077792,731,8220C→G21.0% 188.1 ‑2.6 ‑0.1272V336V (GTG→GTCclpBprotein disaggregation chaperone
Reads supporting (aligned to +/- strand):  new base (25/32):  ref base (123/92):  total (148/124)
Fisher's exact test strand distribution p-value = 7.51e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.96e-01

Only 100 of 295 total aligned reads displayed.

TTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCA  >  NC_007779/2731722‑2731922
                                                                                                    |                                                                                                    
  gCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAacac                                                                                                    >  1:6309100/1‑101 (MQ=255)
  gCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAacac                                                                                                    >  1:14716088/1‑101 (MQ=255)
  gCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAacac                                                                                                    >  1:6168506/1‑101 (MQ=255)
  gCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAacac                                                                                                    <  1:6616171/101‑1 (MQ=255)
      ggTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTc                                                                                                >  1:9418963/1‑101 (MQ=255)
         cAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTgg                                                                                             <  1:244596/101‑1 (MQ=255)
          aGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGa                                                                                            >  1:1634462/1‑101 (MQ=255)
           gTAATTTGCCCATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaa                                                                                           <  1:3782802/101‑1 (MQ=255)
            tAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaaa                                                                                          >  1:6971277/1‑101 (MQ=255)
            tAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaaa                                                                                          <  1:786517/101‑1 (MQ=255)
              aTTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacg                                                                                        <  1:1423698/101‑1 (MQ=255)
              aTTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacg                                                                                        >  1:6721607/1‑101 (MQ=255)
                   aacaccgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAacgacgct                                                                                   <  1:12898653/95‑3 (MQ=255)
                    acaccgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTc                                                                                  >  1:557566/6‑101 (MQ=255)
                    acaTGGTGGTTCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTc                                                                                  <  1:338479/101‑1 (MQ=255)
                      aTGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCa                                                                                <  1:13713379/101‑1 (MQ=255)
                       tggtggTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAg                                                                               <  1:10254947/101‑1 (MQ=255)
                         gttgtgCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAgcg                                                                             >  1:13985040/4‑101 (MQ=255)
                         gtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAgcg                                                                             >  1:11855173/1‑101 (MQ=255)
                         gtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAgcg                                                                             <  1:331213/101‑1 (MQ=255)
                          tggtggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAgcgc                                                                            >  1:2130959/10‑101 (MQ=255)
                          tggtggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAgcgc                                                                            <  1:13832998/92‑1 (MQ=255)
                           ggtggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCgca                                                                           <  1:1041801/93‑1 (MQ=255)
                           ggtgCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGGgca                                                                           <  1:5898927/101‑1 (MQ=255)
                           ggtgCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCgca                                                                           >  1:4531180/1‑101 (MQ=255)
                            gtgCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCgcag                                                                          >  1:4947044/1‑101 (MQ=255)
                             tggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCgcagc                                                                         <  1:1347839/95‑1 (MQ=255)
                             tgCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCgcagc                                                                         >  1:3974981/1‑101 (MQ=255)
                                aaTTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATc                                                                      >  1:13312001/1‑101 (MQ=255)
                                 cctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATcc                                                                     >  1:8970068/3‑100 (MQ=255)
                                 cctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATcc                                                                     >  1:2854144/3‑100 (MQ=255)
                                 aTTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCt                                                                     <  1:10924969/101‑1 (MQ=255)
                                 aTTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCt                                                                     <  1:13144125/101‑1 (MQ=255)
                                  ttCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtt                                                                    >  1:9830855/1‑101 (MQ=255)
                                  ctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATcct                                                                    <  1:6002077/100‑3 (MQ=255)
                                  ctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATcct                                                                    <  1:7351771/100‑3 (MQ=255)
                                   tCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCctt                                                                   <  1:6848798/101‑1 (MQ=255)
                                   tCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttt                                                                   <  1:9649888/101‑1 (MQ=255)
                                    cGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCcttc                                                                  >  1:4745324/1‑100 (MQ=255)
                                    cGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCcttc                                                                  <  1:6914556/101‑2 (MQ=255)
                                    cGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtttt                                                                  >  1:5212049/1‑101 (MQ=255)
                                     gTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCcttct                                                                 <  1:638473/101‑3 (MQ=255)
                                     gTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttt                                                                 >  1:12373207/1‑101 (MQ=255)
                                     gTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttt                                                                 >  1:7593682/1‑101 (MQ=255)
                                     gTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttt                                                                 <  1:2703656/101‑1 (MQ=255)
                                      tagcGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTc                                                                <  1:11187773/98‑1 (MQ=255)
                                       aaCGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCa                                                               <  1:10396981/101‑1 (MQ=255)
                                        gcGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTcga                                                              >  1:776164/2‑99 (MQ=255)
                                        gcGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTcga                                                              >  1:11782132/2‑99 (MQ=255)
                                         cGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAAt                                                             <  1:14644926/101‑1 (MQ=255)
                                         cGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTCGAt                                                             <  1:7959267/101‑1 (MQ=255)
                                         cGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTCGAt                                                             >  1:4316948/1‑101 (MQ=255)
                                          gTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACCGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATg                                                            <  1:10267922/101‑1 (MQ=255)
                                          gTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCACTg                                                            <  1:10353314/101‑1 (MQ=255)
                                          gTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATg                                                            >  1:10208592/1‑101 (MQ=255)
                                          gTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATg                                                            >  1:3243268/1‑101 (MQ=255)
                                          gTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATg                                                            <  1:10757473/101‑1 (MQ=255)
                                          gTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATg                                                            >  1:184981/1‑101 (MQ=255)
                                            tCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTa                                                          <  1:438918/101‑1 (MQ=255)
                                              tttCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACt                                                        <  1:11577086/101‑1 (MQ=255)
                                              tttCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACt                                                        <  1:8197039/101‑1 (MQ=255)
                                              cttCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTCGATGTACt                                                        <  1:1296704/100‑1 (MQ=255)
                                              cttCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTCGATGTACt                                                        <  1:3924809/100‑1 (MQ=255)
                                               ttCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTg                                                       >  1:13445061/1‑101 (MQ=255)
                                                tCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTgg                                                      >  1:3461143/1‑101 (MQ=255)
                                                    gCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGat                                                  >  1:10092543/1‑101 (MQ=255)
                                                     gcACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGata                                                 <  1:12092862/100‑1 (MQ=255)
                                                      cACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGatat                                                >  1:12519251/1‑101 (MQ=255)
                                                      cACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGatat                                                <  1:5635458/101‑1 (MQ=255)
                                                       aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATAtt                                               >  1:3843859/1‑101 (MQ=255)
                                                       aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATAtt                                               >  1:2291641/1‑101 (MQ=255)
                                                        cGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTc                                              >  1:1429720/1‑101 (MQ=255)
                                                        cGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTc                                              >  1:2190282/1‑101 (MQ=255)
                                                          cAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGt                                            >  1:13538667/1‑101 (MQ=255)
                                                          cAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGt                                            >  1:4497180/1‑101 (MQ=255)
                                                           aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                                           >  1:1781903/1‑101 (MQ=255)
                                                           aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                                           <  1:11249491/101‑1 (MQ=255)
                                                                  cAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtc                                    >  1:8678350/1‑101 (MQ=255)
                                                                     tGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                 >  1:291049/1‑101 (MQ=255)
                                                                      ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTgg                                <  1:5820641/101‑1 (MQ=255)
                                                                       gTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGc                               <  1:3562160/101‑1 (MQ=255)
                                                                         aTCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAc                             >  1:7042504/1‑101 (MQ=255)
                                                                          tCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc                            >  1:9508413/1‑101 (MQ=255)
                                                                            ttCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTa                          >  1:4119137/1‑101 (MQ=255)
                                                                            ttCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTa                          >  1:10973711/1‑101 (MQ=255)
                                                                                  aGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGt                    >  1:8035121/1‑101 (MQ=255)
                                                                                  aGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGt                    >  1:2601740/1‑101 (MQ=255)
                                                                                   gAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTg                   <  1:3054910/101‑1 (MQ=255)
                                                                                    aaGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGc                  >  1:10032978/1‑101 (MQ=255)
                                                                                     aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCa                 <  1:1275782/101‑1 (MQ=255)
                                                                                       gCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATGTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAAt               >  1:15068505/1‑101 (MQ=255)
                                                                                       gCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAAt               >  1:923373/1‑101 (MQ=255)
                                                                                       gCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAAt               >  1:12005794/1‑101 (MQ=255)
                                                                                          cGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTca            >  1:14130426/1‑101 (MQ=255)
                                                                                          cGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTca            <  1:9314927/101‑1 (MQ=255)
                                                                                           ggCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcac           <  1:13794634/101‑1 (MQ=255)
                                                                                            gCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcacc          <  1:8583731/101‑1 (MQ=255)
                                                                                            gCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcacc          <  1:12333402/101‑1 (MQ=255)
                                                                                                caaacaCTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCAcgc      >  1:4272066/1‑101 (MQ=255)
                                                                                                    cacTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCa  <  1:4070945/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
TTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCA  >  NC_007779/2731722‑2731922

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