Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,213,001 T→C 1.6% intergenic (+245/+256) mcrA → / ← elbA 5‑methylcytosine‑specific restriction endonuclease B/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,213,0010T→C 1.6% 7.1 0.0 ‑0.4365intergenic (+245/+256)mcrA/elbA5‑methylcytosine‑specific restriction endonuclease B/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/2):  ref base (217/142):  total (221/144)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.70e-01

Only 100 of 383 total aligned reads displayed.

TTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGC  >  NC_007779/1212901‑1213100
                                                                                                    |                                                                                                   
  ataGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCatat                                                                                                   >  1:18511538/1‑101 (MQ=255)
    aGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGTGCCATATgg                                                                                                 >  1:18096697/1‑101 (MQ=255)
    aGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATgg                                                                                                 <  1:15877483/101‑1 (MQ=255)
      aaaTGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGt                                                                                               <  1:2655850/101‑1 (MQ=255)
        aTGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTg                                                                                             >  1:9363621/1‑101 (MQ=255)
        aTGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTg                                                                                             >  1:3392914/1‑101 (MQ=255)
           cTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGAc                                                                                          >  1:2964602/1‑101 (MQ=255)
           cTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGAc                                                                                          >  1:17069775/1‑101 (MQ=255)
             aaGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTg                                                                                        >  1:4915897/1‑101 (MQ=255)
             aaGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTg                                                                                        >  1:12795628/1‑101 (MQ=255)
              aGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGa                                                                                       <  1:14326299/101‑1 (MQ=255)
               gggCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGaa                                                                                      >  1:8415119/1‑101 (MQ=255)
               gggCTAATGAGTTAGATGCAAATTAGCAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTATTGCTGCGCCATATGGGTTGGACTGaa                                                                                      <  1:13433921/101‑1 (MQ=255)
                 gCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGc                                                                                    <  1:14430174/101‑1 (MQ=255)
                 gCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGc                                                                                    >  1:7796204/1‑101 (MQ=255)
                 gCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGc                                                                                    <  1:463109/101‑1 (MQ=255)
                 gCTAATGAGTTATATGCAAACTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGc                                                                                    <  1:16632945/101‑1 (MQ=255)
                  cTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCg                                                                                   >  1:13171792/1‑101 (MQ=255)
                   tAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCgg                                                                                  >  1:16926876/1‑101 (MQ=255)
                     aTGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCt                                                                                >  1:7014948/1‑101 (MQ=255)
                      tGAGTTATATGCAAAGTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTg                                                                               <  1:12583967/101‑1 (MQ=255)
                       gAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGa                                                                              >  1:8508324/1‑101 (MQ=255)
                           tataTGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTg                                                                          >  1:5506115/1‑101 (MQ=255)
                           tataTGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTg                                                                          >  1:12291355/1‑101 (MQ=255)
                            atatGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGa                                                                         >  1:6834677/1‑101 (MQ=255)
                              atGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGAtt                                                                       >  1:668901/1‑101 (MQ=255)
                               tGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTg                                                                      >  1:15707816/1‑101 (MQ=255)
                                gCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGTGGCTGACCTGATTGt                                                                     >  1:4495279/1‑101 (MQ=255)
                                gCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGt                                                                     <  1:5794195/101‑1 (MQ=255)
                                 cAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGtt                                                                    <  1:15495748/101‑1 (MQ=255)
                                     ttAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAg                                                                >  1:15082841/1‑101 (MQ=255)
                                     ttAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAg                                                                <  1:15342409/101‑1 (MQ=255)
                                      tAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                               >  1:124846/1‑101 (MQ=255)
                                       aGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGt                                                              >  1:11521587/1‑101 (MQ=255)
                                        gTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                             <  1:224623/101‑1 (MQ=255)
                                        gTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                             <  1:11750531/101‑1 (MQ=255)
                                         tAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTAt                                                            >  1:9717579/1‑101 (MQ=255)
                                         tAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTAt                                                            >  1:10375243/1‑101 (MQ=255)
                                          aaaaTTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATg                                                           >  1:8824508/1‑101 (MQ=255)
                                          aaaaTTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATg                                                           >  1:3919351/1‑101 (MQ=255)
                                           aaaTTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATgg                                                          <  1:14795885/101‑1 (MQ=255)
                                            aaTTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGa                                                         >  1:14195178/1‑101 (MQ=255)
                                            aaTTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGCATGGa                                                         >  1:5136983/1‑101 (MQ=255)
                                            aaTTATGTTGCTATGTCAGATAGTTACGATTTAGTAATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGAATGGa                                                         >  1:4338485/1‑101 (MQ=255)
                                             aTTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGaa                                                        >  1:9224283/1‑101 (MQ=255)
                                               tATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAgg                                                      >  1:3800987/1‑101 (MQ=255)
                                                 tGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcg                                                    >  1:7674461/1‑101 (MQ=255)
                                                 tGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcg                                                    >  1:13471757/1‑101 (MQ=255)
                                                 tGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcg                                                    <  1:14670932/101‑1 (MQ=255)
                                                 tGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcg                                                    <  1:13291307/101‑1 (MQ=255)
                                                     gCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaat                                                <  1:16318045/101‑1 (MQ=255)
                                                      cTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatc                                               >  1:16419185/1‑101 (MQ=255)
                                                      cTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatc                                               <  1:13736908/101‑1 (MQ=255)
                                                      cTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatc                                               >  1:10268568/1‑101 (MQ=255)
                                                       taagTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatca                                              >  1:17529767/4‑101 (MQ=255)
                                                       tATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatca                                              >  1:5011232/1‑101 (MQ=255)
                                                       tATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatca                                              <  1:18355571/101‑1 (MQ=255)
                                                        aTGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGATTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatcaa                                             >  1:11336227/1‑101 (MQ=255)
                                                        aTGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatcaa                                             >  1:6299506/1‑101 (MQ=255)
                                                        aTGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatcaa                                             >  1:10769935/1‑101 (MQ=255)
                                                          gTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATg                                           <  1:10672902/101‑1 (MQ=255)
                                                           tCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGc                                          <  1:11966486/101‑1 (MQ=255)
                                                            cAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATGTGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGcc                                         >  1:13982082/1‑101 (MQ=255)
                                                             aGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCa                                        >  1:2786132/1‑101 (MQ=255)
                                                             aGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCa                                        >  1:17182306/1‑101 (MQ=255)
                                                             aGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCa                                        >  1:12107696/1‑101 (MQ=255)
                                                              gATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCaa                                       >  1:2730806/1‑101 (MQ=255)
                                                              gATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCaa                                       <  1:3298390/101‑1 (MQ=255)
                                                              gATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCaa                                       >  1:16108635/1‑101 (MQ=255)
                                                              gATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCaa                                       <  1:10155793/101‑1 (MQ=255)
                                                               aTAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAg                                      >  1:8210743/1‑101 (MQ=255)
                                                               aTAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAg                                      >  1:5974238/1‑101 (MQ=255)
                                                                tAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGa                                     >  1:5703063/1‑101 (MQ=255)
                                                                tAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATAAATGCCaaaa                                     >  1:18255158/1‑99 (MQ=255)
                                                                  gTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGAcc                                   <  1:7535727/101‑1 (MQ=255)
                                                                  gTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGAcc                                   <  1:11819906/101‑1 (MQ=255)
                                                                    tACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGt                                 >  1:6562923/1‑101 (MQ=255)
                                                                     aCGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTa                                >  1:4050010/1‑101 (MQ=255)
                                                                      cGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTaa                               >  1:11177141/1‑101 (MQ=255)
                                                                      cGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTaa                               >  1:1302004/1‑101 (MQ=255)
                                                                       gATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAc                              >  1:5262786/1‑101 (MQ=255)
                                                                       gATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAc                              >  1:14938055/1‑101 (MQ=255)
                                                                       gATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAc                              >  1:7212004/1‑101 (MQ=255)
                                                                        aTTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTATCt                             <  1:10967846/101‑1 (MQ=255)
                                                                         tttAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACtt                            <  1:12374780/101‑1 (MQ=255)
                                                                             gTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAAGCGCAATCAATGCCAAGACCGTAACTTATGa                        >  1:2864143/1‑101 (MQ=255)
                                                                              tCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGAc                       >  1:7386088/1‑101 (MQ=255)
                                                                              tCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGAc                       >  1:17681116/1‑101 (MQ=255)
                                                                              tCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGAc                       <  1:17351436/101‑1 (MQ=255)
                                                                               cATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACt                      >  1:8724856/1‑101 (MQ=255)
                                                                                aTCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACtt                     >  1:2300071/1‑101 (MQ=255)
                                                                                aTCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACtt                     <  1:10018198/101‑1 (MQ=255)
                                                                                  ctaactaaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCg                   >  1:17817692/1‑101 (MQ=255)
                                                                                       taaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGt              >  1:2523521/1‑101 (MQ=255)
                                                                                       taaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGt              <  1:17931920/101‑1 (MQ=255)
                                                                                           gctgcGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGa          >  1:13847149/1‑101 (MQ=255)
                                                                                              gcgcCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATgg       <  1:2600620/101‑1 (MQ=255)
                                                                                               cgcCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGa      >  1:9040074/1‑101 (MQ=255)
                                                                                               cgcCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGa      >  1:9825055/1‑101 (MQ=255)
                                                                                                 ccATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAg    <  1:10525946/101‑1 (MQ=255)
                                                                                                    |                                                                                                   
TTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGC  >  NC_007779/1212901‑1213100

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