Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 13,403 A→G 12.9% K414R (AAA→AGA)  dnaK → chaperone Hsp70, co‑chaperone with DnaJ

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_00777913,4030A→G12.9% 89.0 ‑2.3 ‑0.0262K414R (AAA→AGA) dnaKchaperone Hsp70, co‑chaperone with DnaJ
Reads supporting (aligned to +/- strand):  new base (21/13):  ref base (105/123):  total (126/136)
Fisher's exact test strand distribution p-value = 9.95e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.01e-01

Only 100 of 286 total aligned reads displayed.

GTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGT  >  NC_007779/13303‑13503
                                                                                                    |                                                                                                     
 ttCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaa                                                                                                      <  1:10055631/101‑1 (MQ=255)
 ttCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaa                                                                                                      <  1:6382793/101‑1 (MQ=255)
  tCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaaa                                                                                                     <  1:6179805/101‑1 (MQ=255)
    tGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAac                                                                                                   <  1:9297419/101‑1 (MQ=255)
    tGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAac                                                                                                   >  1:7631309/1‑101 (MQ=255)
     gACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAaca                                                                                                  <  1:11052334/101‑1 (MQ=255)
      gcTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTGCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAacac                                                                                                 <  1:3588026/100‑1 (MQ=255)
        cggTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAAcacca                                                                                               >  1:11913287/2‑101 (MQ=255)
          gTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGGGCAACACCACt                                                                                             >  1:4047689/1‑101 (MQ=255)
          gTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACt                                                                                             <  1:8011773/101‑1 (MQ=255)
           tGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTa                                                                                            <  1:10919644/101‑1 (MQ=255)
            gACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAt                                                                                           >  1:7306029/1‑101 (MQ=255)
             aCGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATc                                                                                          <  1:2327351/101‑1 (MQ=255)
             aCGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATc                                                                                          >  1:7795509/1‑101 (MQ=255)
               gTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATccc                                                                                        >  1:1978430/1‑101 (MQ=255)
                tAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCg                                                                                       <  1:4147142/101‑1 (MQ=255)
                tAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCg                                                                                       >  1:1049176/1‑101 (MQ=255)
                tAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCg                                                                                       <  1:1619694/101‑1 (MQ=255)
                 aaaaGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGa                                                                                      >  1:9031682/1‑101 (MQ=255)
                    aGACGTACTGCTGCTGTACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCa                                                                                   >  1:12469359/1‑101 (MQ=255)
                    aGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCa                                                                                   <  1:5047557/101‑1 (MQ=255)
                     gACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCaa                                                                                  <  1:12758494/101‑1 (MQ=255)
                     gACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCaa                                                                                  <  1:7566988/101‑1 (MQ=255)
                      aCGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAg                                                                                 >  1:13380766/1‑101 (MQ=255)
                       cGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGc                                                                               <  1:5901343/101‑1 (MQ=255)
                       cGTACTGCTGCTGGACGTCACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGc                                                                               <  1:11272800/101‑1 (MQ=255)
                        gtgctgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcg                                                                              >  1:12400602/4‑100 (MQ=255)
                        gtgctgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcg                                                                              <  1:5010280/98‑2 (MQ=255)
                        gTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGca                                                                              <  1:5334018/101‑1 (MQ=255)
                         tACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGcac                                                                             >  1:9243850/1‑101 (MQ=255)
                          gctgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgtt                                                                            >  1:2565541/2‑98 (MQ=255)
                            tgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGc                                                                          <  1:15105009/101‑1 (MQ=255)
                            tgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGc                                                                          <  1:3525415/101‑1 (MQ=255)
                              ctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCa                                                                        >  1:14437417/1‑101 (MQ=255)
                               tgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAg                                                                       >  1:14557492/1‑101 (MQ=255)
                               tgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAg                                                                       <  1:2401038/101‑1 (MQ=255)
                               tgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttctgag                                                                       >  1:8464595/1‑93 (MQ=255)
                               tgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttcggag                                                                       >  1:8035914/1‑93 (MQ=255)
                                gctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAgg                                                                      >  1:2783318/1‑101 (MQ=255)
                                gctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttcggaga                                                                      >  1:9172838/1‑92 (MQ=255)
                                 ctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgt                                                                     <  1:1330249/101‑1 (MQ=255)
                                 ctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgt                                                                     <  1:8896873/101‑1 (MQ=255)
                                   ggACGTTGCCCCGCTGTCTCTGGGTATCGAAACCAGGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgtgt                                                                   <  1:3945291/101‑1 (MQ=255)
                                   ggACGTTCCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgtgt                                                                   <  1:11480103/101‑1 (MQ=255)
                                    gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGtt                                                                  >  1:9535726/1‑101 (MQ=255)
                                    gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGtt                                                                  >  1:14762747/1‑101 (MQ=255)
                                     aCGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGGGGGTTGACGCCGCTGACGCGAAAAACCCCAATATCCCCACCCAGCCACACCCAGGTGTtc                                                                 >  1:5223523/1‑101 (MQ=255)
                                      cGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtct                                                                <  1:13941451/101‑1 (MQ=255)
                                        ttcccccGCTGTCTGGGGGTATAGAAACCATGGGGTGTGTGAGGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                              <  1:6824919/98‑1 (MQ=255)
                                        ttACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                              <  1:4113451/101‑1 (MQ=255)
                                         tACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTCTGTTGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTa                                                             <  1:10619326/101‑1 (MQ=255)
                                          aCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTAc                                                            <  1:9736830/101‑1 (MQ=255)
                                           ccccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTAcc                                                           >  1:13741995/1‑101 (MQ=255)
                                            cccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCg                                                          >  1:3556963/1‑101 (MQ=255)
                                              cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                        <  1:8973410/101‑1 (MQ=255)
                                              cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                        >  1:10225626/1‑101 (MQ=255)
                                              cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                        <  1:6614166/101‑1 (MQ=255)
                                                   tctcTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACGCCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGa                                                   <  1:10108929/101‑1 (MQ=255)
                                                    ctctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGAc                                                  <  1:7459250/101‑1 (MQ=255)
                                                    ctctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGAc                                                  >  1:6881220/1‑101 (MQ=255)
                                                    ctctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGAc                                                  >  1:4434496/1‑101 (MQ=255)
                                                    ctctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGAc                                                  <  1:5122354/101‑1 (MQ=255)
                                                     tctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACa                                                 >  1:10173480/1‑101 (MQ=255)
                                                           tATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGt                                           <  1:11207601/101‑1 (MQ=255)
                                                             tCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCt                                         >  1:6243661/1‑101 (MQ=255)
                                                                 aaCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCgg                                     >  1:7057484/1‑101 (MQ=255)
                                                                 aaCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCgg                                     <  1:1939455/101‑1 (MQ=255)
                                                                   ccATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTa                                   <  1:6100003/101‑1 (MQ=255)
                                                                   ccATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTa                                   <  1:14889407/101‑1 (MQ=255)
                                                                      tGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAcc                                <  1:4817368/101‑1 (MQ=255)
                                                                        ggCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccat                              <  1:12180131/101‑1 (MQ=255)
                                                                         gCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatc                             >  1:14149282/1‑101 (MQ=255)
                                                                           ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcca                           <  1:1549958/101‑1 (MQ=255)
                                                                           ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcca                           >  1:11109438/1‑101 (MQ=255)
                                                                           ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcca                           >  1:14452386/1‑101 (MQ=255)
                                                                              gtgATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCAtgt                        <  1:3435101/101‑1 (MQ=255)
                                                                               tgatgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCAtgtg                       >  1:12303179/1‑101 (MQ=255)
                                                                                gatgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgc                      <  1:8288275/101‑1 (MQ=255)
                                                                                 atgaCGCCGCTGATCGCGAAAACCACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTCCCGCTGAAGACAACCAGTCTGCGCTAACCATCCATGtgct                     <  1:1447114/101‑1 (MQ=255)
                                                                                   gacgacGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctgc                   <  1:399659/101‑1 (MQ=255)
                                                                                   gacgacGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctgc                   <  1:4336352/101‑1 (MQ=255)
                                                                                     cgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAg                 >  1:3745081/1‑101 (MQ=255)
                                                                                     cgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAg                 <  1:588362/101‑1 (MQ=255)
                                                                                         gCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTg             <  1:2688822/101‑1 (MQ=255)
                                                                                          cTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGa            >  1:3335030/1‑101 (MQ=255)
                                                                                           tGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGaa           >  1:9419460/1‑101 (MQ=255)
                                                                                             aTCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACg         <  1:8806758/101‑1 (MQ=255)
                                                                                               cgcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTa       >  1:8108645/1‑101 (MQ=255)
                                                                                                gcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaa      >  1:2235601/1‑101 (MQ=255)
                                                                                                gcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaa      <  1:3179325/101‑1 (MQ=255)
                                                                                                gcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaa      >  1:12849846/1‑101 (MQ=255)
                                                                                                gcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaa      <  1:13080671/101‑1 (MQ=255)
                                                                                                gcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaa      >  1:6260679/1‑101 (MQ=255)
                                                                                                 cgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaaa     <  1:7249245/101‑1 (MQ=255)
                                                                                                 cgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaaa     >  1:5503701/1‑101 (MQ=255)
                                                                                                  gAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAAc    >  1:796754/1‑101 (MQ=255)
                                                                                                  gAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAAc    >  1:4634959/1‑101 (MQ=255)
                                                                                                   aaaaaCACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACg   >  1:5623174/1‑101 (MQ=255)
                                                                                                    aaaaCACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGt  <  1:11626317/101‑1 (MQ=255)
                                                                                                    aaaaCACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGt  >  1:5603235/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
GTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGT  >  NC_007779/13303‑13503

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.