Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,731,816 A→G 20.8% V338V (GTT→GTC clpB ← protein disaggregation chaperone

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077792,731,8160A→G20.8% 171.7 ‑2.6 ‑0.1275V338V (GTT→GTCclpBprotein disaggregation chaperone
Reads supporting (aligned to +/- strand):  new base (25/32):  ref base (125/92):  total (150/125)
Fisher's exact test strand distribution p-value = 7.33e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.09e-01

Only 100 of 301 total aligned reads displayed.

CAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCAC  >  NC_007779/2731716‑2731916
                                                                                                    |                                                                                                    
caacaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGca                                                                                                      >  1:6781352/1‑101 (MQ=255)
caacaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGca                                                                                                      <  1:12751241/101‑1 (MQ=255)
caacaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGca                                                                                                      >  1:13093935/1‑101 (MQ=255)
 aacaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaa                                                                                                     >  1:4152648/1‑101 (MQ=255)
  acaaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaac                                                                                                    <  1:11049248/101‑1 (MQ=255)
   caaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaaca                                                                                                   >  1:11126882/1‑101 (MQ=255)
    aaTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaacaa                                                                                                  <  1:10780410/101‑1 (MQ=255)
     aTTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCaacaaa                                                                                                 <  1:13582695/101‑1 (MQ=255)
      ttGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCaacaaac                                                                                                >  1:10376195/1‑101 (MQ=255)
        gCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAacac                                                                                              >  1:9461681/1‑101 (MQ=255)
        gCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAacac                                                                                              >  1:6168506/1‑101 (MQ=255)
        gCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAacac                                                                                              >  1:6309100/1‑101 (MQ=255)
          cGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACtt                                                                                            >  1:3243464/1‑101 (MQ=255)
          cGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCCGCAACAAACACtt                                                                                            <  1:12945544/101‑1 (MQ=255)
              tCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTg                                                                                        <  1:13795399/101‑1 (MQ=255)
              tCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTg                                                                                        >  1:4272446/1‑101 (MQ=255)
               cAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTgg                                                                                       >  1:3081063/1‑101 (MQ=255)
               cAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTgg                                                                                       <  1:244596/101‑1 (MQ=255)
                aGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGa                                                                                      >  1:1634462/1‑101 (MQ=255)
                 gTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaa                                                                                     <  1:4510965/101‑1 (MQ=255)
                  tAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaaa                                                                                    >  1:10233328/1‑101 (MQ=255)
                  tAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaaa                                                                                    <  1:14329200/101‑1 (MQ=255)
                  tAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaaa                                                                                    <  1:5075609/101‑1 (MQ=255)
                  tAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGaaa                                                                                    <  1:786517/101‑1 (MQ=255)
                   aaTTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAac                                                                                   >  1:1777397/1‑101 (MQ=255)
                    aTTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacg                                                                                  >  1:293250/1‑101 (MQ=255)
                     tttGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacga                                                                                 >  1:3350375/1‑101 (MQ=255)
                     tttGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacga                                                                                 <  1:10255658/101‑1 (MQ=255)
                      ttGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacgac                                                                                >  1:3929594/1‑101 (MQ=255)
                        gCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGt                                                                              >  1:14056125/1‑101 (MQ=255)
                          acaccgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTc                                                                            >  1:557566/6‑101 (MQ=255)
                            accgtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCa                                                                          <  1:7033976/98‑1 (MQ=255)
                              ggtggtGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAgc                                                                        <  1:2752183/101‑1 (MQ=255)
                               gtggtgGACCTCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAgcg                                                                       <  1:14788212/101‑1 (MQ=255)
                                tggtggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAgcgc                                                                      >  1:8557327/10‑101 (MQ=255)
                                tggtggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAgcgc                                                                      >  1:2130959/10‑101 (MQ=255)
                                tggtgCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAgcgc                                                                      <  1:15156776/101‑1 (MQ=255)
                                 ggtgCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGGgca                                                                     <  1:5898927/101‑1 (MQ=255)
                                  gtgCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCgcag                                                                    >  1:4947044/1‑101 (MQ=255)
                                   tggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCgcagc                                                                   <  1:6888834/95‑1 (MQ=255)
                                   tggacctCGTAGCATTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCgcagc                                                                   <  1:13206683/95‑1 (MQ=255)
                                    ggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCgcagca                                                                  >  1:12896221/6‑101 (MQ=255)
                                    ggacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCgcagca                                                                  >  1:10101709/6‑101 (MQ=255)
                                     gacctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCAt                                                                 >  1:3502562/5‑101 (MQ=255)
                                     cAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCAt                                                                 <  1:6678921/101‑1 (MQ=255)
                                      acctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATc                                                                <  1:4048185/98‑1 (MQ=255)
                                       cctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATcc                                                               >  1:14257313/3‑100 (MQ=255)
                                       aTTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCt                                                               >  1:3793773/1‑101 (MQ=255)
                                       aTTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCt                                                               >  1:4533461/1‑101 (MQ=255)
                                       aTTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCt                                                               <  1:10924969/101‑1 (MQ=255)
                                        ttCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtt                                                              <  1:1235999/101‑1 (MQ=255)
                                        ctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATcct                                                              <  1:6002077/100‑3 (MQ=255)
                                        ctCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATcct                                                              <  1:7351771/100‑3 (MQ=255)
                                         tCGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCctt                                                             <  1:8635194/101‑1 (MQ=255)
                                          cGTAGCGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCcttc                                                            >  1:4745324/1‑100 (MQ=255)
                                           gTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttt                                                           >  1:7593682/1‑101 (MQ=255)
                                            tagcGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTc                                                          >  1:9172028/4‑101 (MQ=255)
                                            tAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTc                                                          >  1:6967169/1‑101 (MQ=255)
                                             agcGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTcg                                                         >  1:9773228/3‑100 (MQ=255)
                                             aaCGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCa                                                         >  1:11214541/1‑101 (MQ=255)
                                              gcGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTcga                                                        >  1:776164/2‑99 (MQ=255)
                                               cGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAAt                                                       <  1:14644926/101‑1 (MQ=255)
                                               cGTTCCTTCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTCGAt                                                       <  1:7959267/101‑1 (MQ=255)
                                                gTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATg                                                      >  1:184981/1‑101 (MQ=255)
                                                 ttCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGt                                                     <  1:3387193/101‑1 (MQ=255)
                                                    tttCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACt                                                  <  1:8197039/101‑1 (MQ=255)
                                                    cttCAGGCCACGCAGAATGCCAATGGCGTCCTCAACGCTCGGCTCACCGACAAAGACCTGCTGGAAACGACGCTCCAGCGCAGCATCCTTCTCGATGTACt                                                  <  1:12272793/100‑1 (MQ=255)
                                                     ttCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTg                                                 <  1:12050749/101‑1 (MQ=255)
                                                           ccACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGata                                           <  1:8814168/101‑1 (MQ=255)
                                                            cACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGatat                                          <  1:5635458/101‑1 (MQ=255)
                                                            cACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGatat                                          >  1:1184558/1‑101 (MQ=255)
                                                             aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATAtt                                         >  1:2291641/1‑101 (MQ=255)
                                                              cGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTc                                        >  1:2190282/1‑101 (MQ=255)
                                                                cAGAATCGCAATGATATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGt                                      >  1:10334618/1‑101 (MQ=255)
                                                                 aGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                                     >  1:3866809/1‑101 (MQ=255)
                                                                  gAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCa                                    >  1:13754889/1‑101 (MQ=255)
                                                                    aTCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAg                                  <  1:10424333/101‑1 (MQ=255)
                                                                    aTCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAg                                  <  1:11616848/101‑1 (MQ=255)
                                                                    aTCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAg                                  >  1:11983070/1‑101 (MQ=255)
                                                                      cGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcg                                >  1:11388235/1‑101 (MQ=255)
                                                                        cAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtc                              <  1:13511617/101‑1 (MQ=255)
                                                                         aaTGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtcg                             >  1:3618518/1‑101 (MQ=255)
                                                                          aTGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtcgt                            <  1:43315/101‑1 (MQ=255)
                                                                            ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTgg                          <  1:5820641/101‑1 (MQ=255)
                                                                            ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTgg                          <  1:1039785/101‑1 (MQ=255)
                                                                              tATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCa                        >  1:9277159/1‑101 (MQ=255)
                                                                                tCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc                      >  1:13457927/1‑101 (MQ=255)
                                                                                 cTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                     <  1:12502357/101‑1 (MQ=255)
                                                                                 cTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                     <  1:13845726/101‑1 (MQ=255)
                                                                                  ttCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTa                    >  1:10973711/1‑101 (MQ=255)
                                                                                       cAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCATTGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAg               >  1:14169530/1‑101 (MQ=255)
                                                                                        aGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGt              >  1:8035121/1‑101 (MQ=255)
                                                                                        aGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGt              >  1:2601740/1‑101 (MQ=255)
                                                                                         gAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTg             <  1:9341405/101‑1 (MQ=255)
                                                                                          aaGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGc            >  1:10032978/1‑101 (MQ=255)
                                                                                               tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTc       >  1:4403406/1‑101 (MQ=255)
                                                                                               tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTc       >  1:9510746/1‑101 (MQ=255)
                                                                                                 ggCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcac     >  1:2899286/1‑101 (MQ=255)
                                                                                                 ggCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcac     <  1:13332497/101‑1 (MQ=255)
                                                                                                   caacaaACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGCCACCTACGCAGTGCAATTcacca   >  1:6394874/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
CAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCAC  >  NC_007779/2731716‑2731916

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.