Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 13,402 A→C 13.0% K414Q (AAA→CAA)  dnaK → chaperone Hsp70, co‑chaperone with DnaJ

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_00777913,4020A→C13.0% 113.0 ‑2.3 ‑0.1262K414Q (AAA→CAA) 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 = 7.20e-01

Only 100 of 285 total aligned reads displayed.

TGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACG  >  NC_007779/13302‑13502
                                                                                                    |                                                                                                    
 gtTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaa                                                                                                     >  1:4410381/1‑101 (MQ=255)
 gtTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaa                                                                                                     <  1:390984/101‑1 (MQ=255)
  ttCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaa                                                                                                    <  1:10055631/101‑1 (MQ=255)
   tCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaaa                                                                                                   <  1:6179805/101‑1 (MQ=255)
     tGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAac                                                                                                 >  1:7631309/1‑101 (MQ=255)
     tGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAac                                                                                                 <  1:9297419/101‑1 (MQ=255)
      gACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAaca                                                                                                <  1:5479941/101‑1 (MQ=255)
       gcTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTGCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAacac                                                                                               <  1:3588026/100‑1 (MQ=255)
           gTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGGGCAACACCACt                                                                                           >  1:4047689/1‑101 (MQ=255)
           gTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACt                                                                                           <  1:8011773/101‑1 (MQ=255)
            tGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTa                                                                                          >  1:13562223/1‑101 (MQ=255)
            tGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTa                                                                                          >  1:3531127/1‑101 (MQ=255)
              aCGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATc                                                                                        >  1:7795509/1‑101 (MQ=255)
               ggtcaaggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATcc                                                                                       <  1:8090745/94‑1 (MQ=255)
               cGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATcc                                                                                       <  1:15161181/101‑1 (MQ=255)
                gtcaaggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATccc                                                                                      <  1:12624311/95‑1 (MQ=255)
                 tAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCg                                                                                     <  1:1619694/101‑1 (MQ=255)
                 tAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCg                                                                                     >  1:1049176/1‑101 (MQ=255)
                   aaaGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGAc                                                                                   >  1:5606378/1‑101 (MQ=255)
                    aggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGAcc                                                                                  >  1:4584375/3‑101 (MQ=255)
                     aGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCa                                                                                 <  1:5047557/101‑1 (MQ=255)
                      gACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCaa                                                                                <  1:7566988/101‑1 (MQ=255)
                      gACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCaa                                                                                <  1:3387621/101‑1 (MQ=255)
                        cGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGc                                                                              <  1:5638532/101‑1 (MQ=255)
                        cGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGc                                                                              <  1:5901343/101‑1 (MQ=255)
                        cGTACTGCTGCTGGACGTCACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGc                                                                              <  1:11272800/101‑1 (MQ=255)
                         gtgctgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcg                                                                             >  1:2346420/4‑100 (MQ=255)
                         gtgctgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcg                                                                             <  1:5010280/98‑2 (MQ=255)
                         gTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGca                                                                             >  1:10916866/1‑101 (MQ=255)
                          tACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGcac                                                                            >  1:9243850/1‑101 (MQ=255)
                          tACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGcac                                                                            >  1:5208303/1‑101 (MQ=255)
                           aCTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGcaca                                                                           >  1:2238369/1‑101 (MQ=255)
                             tgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGc                                                                         >  1:13820549/1‑101 (MQ=255)
                             tgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttcg                                                                         <  1:14905944/101‑6 (MQ=255)
                               ctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttcgga                                                                       >  1:14343553/1‑94 (MQ=255)
                                tgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAg                                                                      >  1:414462/1‑101 (MQ=255)
                                tgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttctgag                                                                      >  1:8464595/1‑93 (MQ=255)
                                 gctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAgg                                                                     >  1:2783318/1‑101 (MQ=255)
                                 gctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttcggaga                                                                     >  1:9172838/1‑92 (MQ=255)
                                  ctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgt                                                                    <  1:8896873/101‑1 (MQ=255)
                                  ctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgt                                                                    <  1:1330249/101‑1 (MQ=255)
                                    ggACGTTCCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgtgt                                                                  <  1:11480103/101‑1 (MQ=255)
                                     gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGtt                                                                 >  1:5114144/1‑101 (MQ=255)
                                     gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGtt                                                                 >  1:9535726/1‑101 (MQ=255)
                                       cGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtct                                                               <  1:9106974/101‑1 (MQ=255)
                                        gTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtctc                                                              >  1:6275763/1‑101 (MQ=255)
                                         ttcccccGCTGTCTGGGGGTATAGAAACCATGGGGTGTGTGAGGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                             <  1:6824919/98‑1 (MQ=255)
                                         ttACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                             >  1:14439635/1‑101 (MQ=255)
                                         ttACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                             <  1:4113451/101‑1 (MQ=255)
                                          tACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTa                                                            >  1:4241979/1‑101 (MQ=255)
                                          tACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTCTGTTGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTa                                                            <  1:10619326/101‑1 (MQ=255)
                                           aCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTAc                                                           <  1:9736830/101‑1 (MQ=255)
                                             cccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCg                                                         >  1:3556963/1‑101 (MQ=255)
                                               cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                       <  1:8973410/101‑1 (MQ=255)
                                               cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                       <  1:3205873/101‑1 (MQ=255)
                                               cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                       >  1:10225626/1‑101 (MQ=255)
                                               cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                       <  1:6614166/101‑1 (MQ=255)
                                                   gTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAg                                                   <  1:5257088/101‑1 (MQ=255)
                                                    tctcTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACGCCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGa                                                  <  1:10108929/101‑1 (MQ=255)
                                                     ctctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGAc                                                 >  1:6881220/1‑101 (MQ=255)
                                                      tctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACa                                                >  1:10173480/1‑101 (MQ=255)
                                                        tGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAAc                                              <  1:1469809/101‑1 (MQ=255)
                                                         gggTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAAcc                                             >  1:2125207/1‑101 (MQ=255)
                                                            tATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGt                                          <  1:11207601/101‑1 (MQ=255)
                                                            tATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGt                                          <  1:3437946/101‑1 (MQ=255)
                                                              tCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCt                                        >  1:6243661/1‑101 (MQ=255)
                                                                 aaaCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCg                                     >  1:10985865/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)
                                                                       tGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAcc                               <  1:4817368/101‑1 (MQ=255)
                                                                         ggCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccat                             <  1:12180131/101‑1 (MQ=255)
                                                                         ggCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccat                             <  1:1282302/101‑1 (MQ=255)
                                                                            ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcca                          >  1:11109438/1‑101 (MQ=255)
                                                                            ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcca                          >  1:11762618/1‑101 (MQ=255)
                                                                            ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcca                          >  1:14452386/1‑101 (MQ=255)
                                                                              tgtgATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCAtg                        >  1:2431582/1‑101 (MQ=255)
                                                                                 gatgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgc                     <  1:8288275/101‑1 (MQ=255)
                                                                                  atgaCGCCGCTGATCGCGAAAACCACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTCCCGCTGAAGACAACCAGTCTGCGCTAACCATCCATGtgct                    <  1:1447114/101‑1 (MQ=255)
                                                                                   tgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctg                   >  1:2579362/1‑101 (MQ=255)
                                                                                    gacgacGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctgc                  <  1:399659/101‑1 (MQ=255)
                                                                                    gacgacGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctgc                  <  1:12409759/101‑1 (MQ=255)
                                                                                     acgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCa                 >  1:12307171/1‑101 (MQ=255)
                                                                                     acgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCa                 <  1:7500614/101‑1 (MQ=255)
                                                                                      cgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAg                >  1:3745081/1‑101 (MQ=255)
                                                                                      cgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAg                >  1:1550098/1‑101 (MQ=255)
                                                                                      cgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAg                <  1:588362/101‑1 (MQ=255)
                                                                                        acgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAggg              <  1:7338457/101‑1 (MQ=255)
                                                                                          gCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTg            <  1:2688822/101‑1 (MQ=255)
                                                                                            tGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGaa          <  1:6478791/101‑1 (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:13976369/101‑1 (MQ=255)
                                                                                                  cgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaaa    <  1:7249245/101‑1 (MQ=255)
                                                                                                  cgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaaa    >  1:1484562/1‑101 (MQ=255)
                                                                                                   gAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAAc   >  1:796754/1‑101 (MQ=255)
                                                                                                   gAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAAc   >  1:12068758/1‑101 (MQ=255)
                                                                                                    aaaaaCACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACg  <  1:12873351/101‑1 (MQ=255)
                                                                                                    aaaaaCACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACg  <  1:13193421/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
TGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACG  >  NC_007779/13302‑13502

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