Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 13,309 A→C 11.9% T383P (ACT→CCT)  dnaK → chaperone Hsp70, co‑chaperone with DnaJ

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_00777913,3090A→C11.9% 100.7 ‑0.3 ‑0.0277T383P (ACT→CCT) dnaKchaperone Hsp70, co‑chaperone with DnaJ
Reads supporting (aligned to +/- strand):  new base (14/19):  ref base (115/129):  total (129/148)
Fisher's exact test strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.96e-01

Only 100 of 281 total aligned reads displayed.

CAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACAC  >  NC_007779/13210‑13409
                                                                                                   |                                                                                                    
cAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGAc                                                                                                     >  1:4928667/1‑101 (MQ=255)
cAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGAc                                                                                                     >  1:264197/1‑101 (MQ=255)
 agaagaAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACt                                                                                                    >  1:2336342/1‑101 (MQ=255)
 agaagaAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACt                                                                                                    <  1:6193930/101‑1 (MQ=255)
 agaagaAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACt                                                                                                    >  1:13103877/1‑101 (MQ=255)
  gaagaaAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTg                                                                                                   <  1:7128452/101‑1 (MQ=255)
  gaagaaAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTg                                                                                                   >  1:3028338/1‑101 (MQ=255)
    agaaAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGt                                                                                                 >  1:1584306/1‑101 (MQ=255)
    agaaAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGt                                                                                                 <  1:3673577/101‑1 (MQ=255)
     gaaAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTg                                                                                                >  1:6456212/1‑101 (MQ=255)
     gaaAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTg                                                                                                >  1:6017484/1‑101 (MQ=255)
      aaaGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGa                                                                                               <  1:10451416/101‑1 (MQ=255)
          ttGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTa                                                                                           >  1:3189248/1‑101 (MQ=255)
          ttGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTa                                                                                           <  1:3380700/101‑1 (MQ=255)
            gCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTaaa                                                                                         <  1:2232319/101‑1 (MQ=255)
             cTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTaaaa                                                                                        <  1:7439825/101‑1 (MQ=255)
             cTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTaaaa                                                                                        <  1:6429731/101‑1 (MQ=255)
             cTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTaaaa                                                                                        >  1:216306/1‑101 (MQ=255)
             cTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTaaaa                                                                                        >  1:2961533/1‑101 (MQ=255)
             cTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTaaaa                                                                                        <  1:6065358/101‑1 (MQ=255)
             cTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTaaaa                                                                                        >  1:6595928/1‑101 (MQ=255)
               gAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGa                                                                                      >  1:9749668/1‑101 (MQ=255)
                  ttCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGt                                                                                   <  1:13993395/101‑1 (MQ=255)
                   tCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTa                                                                                  <  1:11702614/101‑1 (MQ=255)
                    cTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTAc                                                                                 <  1:8474643/101‑1 (MQ=255)
                       tGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTActgc                                                                              >  1:11413475/1‑101 (MQ=255)
                        ggTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTActgct                                                                             <  1:5114848/101‑1 (MQ=255)
                        ggTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTActgct                                                                             >  1:1907619/1‑101 (MQ=255)
                          tAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTActgctgc                                                                           >  1:12904751/1‑101 (MQ=255)
                          tAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTActgctgc                                                                           >  1:2426111/1‑101 (MQ=255)
                            aaGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTActgctgctg                                                                         <  1:4592699/101‑1 (MQ=255)
                             agagCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTgg                                                                        >  1:11348207/1‑101 (MQ=255)
                               agCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGAc                                                                      >  1:8408930/1‑101 (MQ=255)
                               agCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGAc                                                                      <  1:8159609/101‑1 (MQ=255)
                                 ccGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAGAAGACGTACTGCTGCTGGACGt                                                                    <  1:3772637/101‑1 (MQ=255)
                                   gcgTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTa                                                                  <  1:1350564/101‑1 (MQ=255)
                                   gcgTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTa                                                                  <  1:12994293/101‑1 (MQ=255)
                                    cgTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTAc                                                                 >  1:7554889/1‑101 (MQ=255)
                                     gTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTAcc                                                                <  1:5394894/101‑1 (MQ=255)
                                       aaaGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTAcccc                                                              <  1:4734563/101‑1 (MQ=255)
                                       aaaGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTAcccc                                                              >  1:5927955/1‑101 (MQ=255)
                                       aaaGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTAcccc                                                              <  1:7771222/101‑1 (MQ=255)
                                       aaaGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTAcccc                                                              <  1:3254698/101‑1 (MQ=255)
                                         aGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGc                                                            >  1:7940059/1‑101 (MQ=255)
                                         aGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGc                                                            >  1:8760063/1‑101 (MQ=255)
                                            ggtgaaCCCGGACGAAGCTGTAGCAATCGGTGCTGGTGGTCAGGGTGGTGTTCTGACTGGTGACATAAAAGACGTCCTGCTGCTGGACGTTACCCCGCTGt                                                         <  1:4851985/97‑1 (MQ=255)
                                            cGTTAACCCGGACGAGGTTGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGt                                                         <  1:9586602/101‑1 (MQ=255)
                                            cGTTAACCCGGACGAGGTTGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGt                                                         <  1:726139/101‑1 (MQ=255)
                                            cGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGt                                                         >  1:6883747/1‑101 (MQ=255)
                                             gTTAACCCGGACGAGGTTGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGtc                                                        <  1:10357763/101‑1 (MQ=255)
                                             gTTAACCCGGACGAGGTTGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGtc                                                        >  1:10669940/1‑101 (MQ=255)
                                                aaCCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTctct                                                     <  1:10751230/101‑1 (MQ=255)
                                                aaCCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTctct                                                     >  1:12470857/1‑101 (MQ=255)
                                                 aCCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTg                                                    <  1:1504022/101‑1 (MQ=255)
                                                  cccGGACGAGGTTGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTgg                                                   >  1:9259169/1‑101 (MQ=255)
                                                  cccGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTgg                                                   >  1:12739826/1‑101 (MQ=255)
                                                    cGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGt                                                 <  1:1677309/101‑1 (MQ=255)
                                                      gACGAGGTTGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTAt                                               >  1:6554387/1‑101 (MQ=255)
                                                      gACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTAt                                               <  1:421853/101‑1 (MQ=255)
                                                      gACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTAt                                               <  1:8236713/101‑1 (MQ=255)
                                                       aCGAGGTTGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTAtt                                              >  1:12219099/1‑100 (MQ=255)
                                                           aGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGaaa                                          <  1:3575592/101‑1 (MQ=255)
                                                            gttGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGGATTGAGAc                                         >  1:5559637/3‑101 (MQ=255)
                                                            gCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAAc                                         <  1:9869435/101‑1 (MQ=255)
                                                             cTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAAcc                                        >  1:14270624/1‑101 (MQ=255)
                                                              tGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCa                                       >  1:10034346/1‑101 (MQ=255)
                                                              tGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCa                                       >  1:8893782/1‑101 (MQ=255)
                                                                 aGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATgg                                    >  1:1866560/1‑101 (MQ=255)
                                                                   gagttggTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGc                                  <  1:11740870/96‑1 (MQ=255)
                                                                   cAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGc                                  <  1:3906441/101‑1 (MQ=255)
                                                                   cAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGc                                  <  1:952803/101‑1 (MQ=255)
                                                                    agttggTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCg                                 <  1:1013998/97‑1 (MQ=255)
                                                                    agttggTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCg                                 <  1:9282745/97‑1 (MQ=255)
                                                                    aaTCGTTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCg                                 <  1:10914067/101‑1 (MQ=255)
                                                                    aaTCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCg                                 <  1:11587731/101‑1 (MQ=255)
                                                                     aTCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTGCCCCGCTGTCTCTGGGTATCGAAACCATGGGCgg                                <  1:9420490/101‑1 (MQ=255)
                                                                     aTCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCgg                                >  1:13795525/1‑101 (MQ=255)
                                                                         gTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGtgtg                            >  1:13204128/1‑101 (MQ=255)
                                                                          tgctgcTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCggggggt                           >  1:5698791/1‑96 (MQ=255)
                                                                          tgctgcTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGtga                           >  1:5048874/1‑101 (MQ=255)
                                                                            ctgctgTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGtgatg                         >  1:5485314/1‑101 (MQ=255)
                                                                             tgctgTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGCGCGGTGtgatga                        >  1:1743376/1‑101 (MQ=255)
                                                                              gctgTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgac                       <  1:13598398/101‑1 (MQ=255)
                                                                               ctgTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgacg                      <  1:1482141/101‑1 (MQ=255)
                                                                               ctgTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgacg                      >  1:6855376/1‑101 (MQ=255)
                                                                               ctgTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgacg                      <  1:8326364/101‑1 (MQ=255)
                                                                               ctgTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgacg                      <  1:13296784/101‑1 (MQ=255)
                                                                                tgTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgacga                     >  1:11913458/1‑101 (MQ=255)
                                                                                 gTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACATACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgacgac                    <  1:10193052/101‑1 (MQ=255)
                                                                                  ttCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGacgacg                   >  1:1211717/1‑101 (MQ=255)
                                                                                  ttCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGacgacg                   <  1:1986991/101‑1 (MQ=255)
                                                                                    cAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCt                 <  1:14179262/101‑1 (MQ=255)
                                                                                     aGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCta                >  1:5273917/1‑100 (MQ=255)
                                                                                          ggtgTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATcgc           <  1:3122716/101‑1 (MQ=255)
                                                                                            tgtTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGa         <  1:12019380/101‑1 (MQ=255)
                                                                                              ttCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaa       >  1:9000334/1‑101 (MQ=255)
                                                                                               tCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaaa      <  1:1244024/101‑1 (MQ=255)
                                                                                                 tGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAac    >  1:9091326/1‑101 (MQ=255)
                                                                                                 tGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAac    >  1:9138817/1‑101 (MQ=255)
                                                                                                 tGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAac    >  1:12811591/1‑101 (MQ=255)
                                                                                                   |                                                                                                    
CAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACAC  >  NC_007779/13210‑13409

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