Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 13,403 A→G 6.1% 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→G 6.1% 30.0 0.0 ‑0.1229K414R (AAA→AGA) dnaKchaperone Hsp70, co‑chaperone with DnaJ
Reads supporting (aligned to +/- strand):  new base (8/6):  ref base (115/100):  total (123/106)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.13e-01

Only 100 of 244 total aligned reads displayed.

GTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGT  >  NC_007779/13303‑13503
                                                                                                    |                                                                                                    
gtTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaa                                                                                                      <  1:3720262/101‑1 (MQ=255)
   cTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaaaa                                                                                                   >  1:14166630/1‑101 (MQ=255)
    tGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAac                                                                                                  >  1:2670663/1‑101 (MQ=255)
     gACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAaca                                                                                                 <  1:2583284/101‑1 (MQ=255)
       cTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAAcacc                                                                                               >  1:4369096/1‑101 (MQ=255)
       cTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAAcacc                                                                                               >  1:13721980/1‑101 (MQ=255)
        tGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAAcacca                                                                                              >  1:637037/1‑101 (MQ=255)
          gTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACt                                                                                            <  1:1319626/101‑1 (MQ=255)
          gTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACt                                                                                            <  1:10748144/101‑1 (MQ=255)
           tGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTa                                                                                           >  1:1103415/1‑101 (MQ=255)
            gACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAt                                                                                          >  1:11432596/1‑101 (MQ=255)
            gACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAt                                                                                          <  1:6574683/101‑1 (MQ=255)
              cGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATcc                                                                                        >  1:11337613/1‑101 (MQ=255)
              cGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATcc                                                                                        >  1:7415675/1‑101 (MQ=255)
               gTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATccc                                                                                       <  1:6507437/101‑1 (MQ=255)
                 caaggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGa                                                                                     >  1:1221415/5‑101 (MQ=255)
                 aaaaGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCTCTATCCCGa                                                                                     <  1:12644049/101‑1 (MQ=255)
                      aCGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAg                                                                                <  1:9406615/101‑1 (MQ=255)
                       cGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGc                                                                               <  1:9780702/101‑1 (MQ=255)
                         tgctgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgt                                                                             <  1:13595067/99‑3 (MQ=255)
                         tACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGcac                                                                             <  1:1692714/101‑1 (MQ=255)
                          gctgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgtt                                                                            <  1:11601029/100‑4 (MQ=255)
                          aCTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGcaca                                                                            <  1:5188014/101‑1 (MQ=255)
                           ctgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAg                                                                           >  1:7515959/1‑101 (MQ=255)
                           ctgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAg                                                                           <  1:9177022/101‑1 (MQ=255)
                           ctgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAg                                                                           >  1:7192022/1‑101 (MQ=255)
                           ctgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttc                                                                           <  1:1455531/101‑5 (MQ=255)
                             gctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGcc                                                                         >  1:13347899/1‑101 (MQ=255)
                             gctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttcgg                                                                         >  1:4513335/1‑95 (MQ=255)
                              ctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCa                                                                        <  1:9619543/101‑1 (MQ=255)
                               tgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAg                                                                       <  1:5451502/101‑1 (MQ=255)
                               tgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAg                                                                       >  1:2689383/1‑101 (MQ=255)
                               tgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttcggag                                                                       >  1:7818699/1‑93 (MQ=255)
                                tggttaCGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAgg                                                                      >  1:13546548/6‑101 (MQ=255)
                                gctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAgg                                                                      >  1:11709306/1‑101 (MQ=255)
                                    gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGtt                                                                  >  1:7980481/1‑101 (MQ=255)
                                     aCGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtc                                                                 >  1:9873677/1‑101 (MQ=255)
                                     aCGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtc                                                                 >  1:3369319/1‑101 (MQ=255)
                                      cGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtct                                                                <  1:7768360/101‑1 (MQ=255)
                                      cGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtct                                                                <  1:4250006/101‑1 (MQ=255)
                                       gTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtctc                                                               <  1:6853188/101‑1 (MQ=255)
                                        ttACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                              <  1:5668362/101‑1 (MQ=255)
                                        ttACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct                                                              >  1:6318838/1‑101 (MQ=255)
                                         tACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTa                                                             >  1:2131481/1‑101 (MQ=255)
                                         tACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTa                                                             >  1:13635822/1‑101 (MQ=255)
                                         tACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTa                                                             >  1:12837455/1‑101 (MQ=255)
                                          aCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTAc                                                            <  1:10037783/101‑1 (MQ=255)
                                            cccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCg                                                          >  1:10072331/1‑101 (MQ=255)
                                             ccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGc                                                         <  1:3005288/101‑1 (MQ=255)
                                              cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                        >  1:7087821/1‑101 (MQ=255)
                                              cGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                        <  1:5409387/101‑1 (MQ=255)
                                              cGCTGTCTCCGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCt                                                        >  1:12465630/1‑101 (MQ=255)
                                                cTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGa                                                      >  1:1913203/1‑101 (MQ=255)
                                                cTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGa                                                      >  1:5300143/1‑101 (MQ=255)
                                                 tGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGaa                                                     <  1:2229231/101‑1 (MQ=255)
                                                  gTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACGCCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAg                                                    <  1:1874535/101‑1 (MQ=255)
                                                  gTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAg                                                    <  1:12191455/101‑1 (MQ=255)
                                                   tctcTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGa                                                   >  1:10362033/1‑101 (MQ=255)
                                                   tatctGGGTATCGAAACCACGGGCGGTGTGATGAAGACGCTGATCCCGAAAAACGCCACTATCCCGACCAAGCTCAGCCAGGTGCTCTCTACCGCTGAAGa                                                   <  1:888346/99‑1 (MQ=255)
                                                       tGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAAc                                               >  1:9581522/1‑101 (MQ=255)
                                                         ggTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCa                                             <  1:12019327/101‑1 (MQ=255)
                                                         ggTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCa                                             <  1:11079124/101‑1 (MQ=255)
                                                           tATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGt                                           <  1:2436827/101‑1 (MQ=255)
                                                           tATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGt                                           >  1:3912804/1‑101 (MQ=255)
                                                            aTCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTc                                          >  1:5560987/1‑101 (MQ=255)
                                                             tCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCt                                         <  1:8236432/101‑1 (MQ=255)
                                                                aaaCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCg                                      >  1:11840085/1‑101 (MQ=255)
                                                                 aaCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCgg                                     >  1:5070243/1‑101 (MQ=255)
                                                                 aaCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCgg                                     >  1:3432602/1‑101 (MQ=255)
                                                                  aCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGt                                    >  1:9446914/1‑101 (MQ=255)
                                                                   ccATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTa                                   <  1:4776183/101‑1 (MQ=255)
                                                                   ccATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTa                                   >  1:9522195/1‑101 (MQ=255)
                                                                     aTGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAc                                 <  1:9651651/101‑1 (MQ=255)
                                                                         gCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatc                             <  1:4928535/101‑1 (MQ=255)
                                                                          cGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcc                            <  1:4659957/101‑1 (MQ=255)
                                                                          cGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcc                            >  1:9935490/1‑101 (MQ=255)
                                                                             tgtgATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCAtg                         >  1:8947951/1‑101 (MQ=255)
                                                                              gtgATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCAtgt                        <  1:8774226/101‑1 (MQ=255)
                                                                                 atgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAGCCATCCATGtgct                     <  1:1237144/101‑1 (MQ=255)
                                                                                 atgaCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgct                     <  1:14046370/101‑1 (MQ=255)
                                                                                   gacgacGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctgc                   >  1:7112276/1‑101 (MQ=255)
                                                                                   gacgacGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctgc                   >  1:5727809/1‑101 (MQ=255)
                                                                                     cgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAg                 >  1:2336781/1‑101 (MQ=255)
                                                                                     cgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAg                 <  1:2946873/101‑1 (MQ=255)
                                                                                     cgacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAg                 >  1:6443665/1‑101 (MQ=255)
                                                                                      gacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAgg                <  1:812653/101‑1 (MQ=255)
                                                                                      gacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAgg                <  1:10139934/101‑1 (MQ=255)
                                                                                       acgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAggg               >  1:13241373/1‑101 (MQ=255)
                                                                                       acgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAggg               >  1:12280971/1‑101 (MQ=255)
                                                                                        cgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGt              >  1:8361050/1‑101 (MQ=255)
                                                                                         gCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTg             >  1:5806693/1‑101 (MQ=255)
                                                                                           tGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGaa           >  1:3633576/1‑101 (MQ=255)
                                                                                           tGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGaa           >  1:5364070/1‑101 (MQ=255)
                                                                                             aTCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACg         >  1:11208468/1‑101 (MQ=255)
                                                                                              tCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGt        >  1:143577/1‑101 (MQ=255)
                                                                                               cgcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTa       >  1:10828732/1‑101 (MQ=255)
                                                                                               cgcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTa       >  1:10412811/1‑101 (MQ=255)
                                                                                                gcgAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTaa      <  1:5740867/101‑1 (MQ=255)
                                                                                                  gAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCGGGGGGACCGAaaag    >  1:765163/1‑100 (MQ=255)
                                                                                                    aaaaCACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGt  <  1:5047720/101‑1 (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 < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.