Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,213,026 C→T 2.1% intergenic (+270/+231) 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,0260C→T 2.1% 5.4 ‑0.4 ‑0.1289intergenic (+270/+231)mcrA/elbA5‑methylcytosine‑specific restriction endonuclease B/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/2):  ref base (152/130):  total (157/132)
Fisher's exact test strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.68e-01

Only 100 of 295 total aligned reads displayed.

GTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTT  >  NC_007779/1212926‑1213126
                                                                                                    |                                                                                                    
gTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcc                                                                                                      >  1:15363073/1‑101 (MQ=255)
  tataTGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTg                                                                                                    >  1:12683429/1‑101 (MQ=255)
  tataTGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTg                                                                                                    <  1:14560439/101‑1 (MQ=255)
   atatGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGa                                                                                                   >  1:6338842/1‑101 (MQ=255)
    tatGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGAt                                                                                                  >  1:4191075/1‑101 (MQ=255)
    tatGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGAt                                                                                                  <  1:9188630/101‑1 (MQ=255)
     atGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGCCCTGAtt                                                                                                 <  1:1547489/101‑1 (MQ=255)
     atGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGAtt                                                                                                 >  1:7608657/1‑101 (MQ=255)
        cAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGtt                                                                                              >  1:4561943/1‑101 (MQ=255)
        cAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGtt                                                                                              >  1:11232476/1‑101 (MQ=255)
        cAAATTAGTAAAATTATCTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGtt                                                                                              <  1:10634202/101‑1 (MQ=255)
         aaaTTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTa                                                                                             >  1:12071308/1‑101 (MQ=255)
         aaaTTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTa                                                                                             >  1:3196074/1‑101 (MQ=255)
           aTTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaaa                                                                                           <  1:3998430/101‑1 (MQ=255)
           aTTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaaa                                                                                           >  1:12314150/1‑101 (MQ=255)
             tAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTATTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGGTTGTTAAAgg                                                                                         >  1:9253131/1‑101 (MQ=255)
             tAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                                                         >  1:2131055/1‑101 (MQ=255)
              aGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGt                                                                                        <  1:9122415/101‑1 (MQ=255)
              aGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGt                                                                                        <  1:7173150/101‑1 (MQ=255)
               gTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                                                       >  1:5659151/1‑101 (MQ=255)
                tAAAATTATGTTGTTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGNGCCATATGGGTTGGACTGAAGCGNNTGACCTGATTGTNAAAGGTAt                                                                                      >  1:12956090/1‑101 (MQ=255)
                tAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTAt                                                                                      <  1:4459018/101‑1 (MQ=255)
                tAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTAt                                                                                      <  1:11464105/101‑1 (MQ=255)
                   aaTTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGa                                                                                   <  1:4897600/101‑1 (MQ=255)
                    aTTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGaa                                                                                  <  1:13107825/101‑1 (MQ=255)
                     ttATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGAAAg                                                                                 >  1:6574501/1‑101 (MQ=255)
                     ttATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGCTAAAGGTATGGAAg                                                                                 >  1:7351969/1‑101 (MQ=255)
                       aTGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgc                                                                               >  1:3330601/1‑101 (MQ=255)
                       aTGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgc                                                                               <  1:14478311/101‑1 (MQ=255)
                         gTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcgc                                                                             >  1:2536779/1‑101 (MQ=255)
                          ttGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGca                                                                            <  1:749771/101‑1 (MQ=255)
                           tGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaa                                                                           >  1:9509481/1‑101 (MQ=255)
                            gCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaat                                                                          <  1:10892864/101‑1 (MQ=255)
                             cTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatc                                                                         >  1:3065688/1‑101 (MQ=255)
                             cTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatc                                                                         >  1:4054660/1‑101 (MQ=255)
                             cTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCACCATATGGGTTGGACTGAAGCGGCTGACCTGCTTGTTAAAGGTATGGAAGGCGcaatc                                                                         >  1:4638581/1‑101 (MQ=255)
                              tATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatca                                                                        >  1:13987218/1‑101 (MQ=255)
                              tATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatca                                                                        >  1:8233651/1‑101 (MQ=255)
                                tGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatcaat                                                                      <  1:2277678/101‑1 (MQ=255)
                                tGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatcaat                                                                      >  1:8834772/1‑101 (MQ=255)
                                tGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatcaat                                                                      <  1:12780684/101‑1 (MQ=255)
                                tGTCAGATAGTTACGACTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatcaat                                                                      >  1:10274144/1‑101 (MQ=255)
                                   cAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGcc                                                                   >  1:13744247/1‑101 (MQ=255)
                                    aGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGCATGGAAGGAGCAATCAATGCCa                                                                  >  1:8871931/1‑101 (MQ=255)
                                     gATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCaa                                                                 <  1:7730792/101‑1 (MQ=255)
                                     gATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCaa                                                                 <  1:13882730/101‑1 (MQ=255)
                                      aTAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAg                                                                >  1:3871899/1‑101 (MQ=255)
                                      aTAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAg                                                                >  1:14121192/1‑101 (MQ=255)
                                       tAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGa                                                               <  1:14758272/101‑1 (MQ=255)
                                       tAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGa                                                               <  1:5786850/101‑1 (MQ=255)
                                        aGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGAc                                                              <  1:9626123/101‑1 (MQ=255)
                                         gTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGAcc                                                             >  1:4945224/1‑101 (MQ=255)
                                          ttACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCg                                                            >  1:1468626/1‑101 (MQ=255)
                                           tACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGt                                                           >  1:7435937/1‑101 (MQ=255)
                                            aCGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTa                                                          >  1:5865735/1‑101 (MQ=255)
                                             cGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTaa                                                         >  1:1889619/1‑101 (MQ=255)
                                                tttAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACtt                                                      >  1:5164600/1‑101 (MQ=255)
                                                tttAGTCATCTAACTAATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAcct                                                      >  1:9485899/1‑99 (MQ=14)
                                                 ttAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTa                                                     <  1:2000489/101‑1 (MQ=255)
                                                   aGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATg                                                   >  1:3716264/1‑101 (MQ=255)
                                                   aGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATg                                                   >  1:6937669/1‑101 (MQ=255)
                                                     tCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGAc                                                 <  1:6638163/101‑1 (MQ=255)
                                                      cATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACt                                                <  1:13284756/101‑1 (MQ=255)
                                                      cATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACt                                                <  1:206900/101‑1 (MQ=255)
                                                      cATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCACAATCAATGCCAAGACCGTAACTTATGACt                                                >  1:7896891/1‑101 (MQ=255)
                                                       gctgagatgaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACtt                                               >  1:14038431/10‑101 (MQ=255)
                                                       aTCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACtt                                               >  1:8645712/1‑101 (MQ=255)
                                                       aTCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACtt                                               >  1:11859234/1‑101 (MQ=255)
                                                         tgagatgaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCg                                             >  1:12012879/8‑101 (MQ=255)
                                                         ctaactaaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCg                                             >  1:4398283/1‑101 (MQ=255)
                                                           agatgaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaa                                           >  1:8081053/6‑101 (MQ=255)
                                                           aactaaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaa                                           >  1:9749070/1‑101 (MQ=255)
                                                            actaaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGTCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAAc                                          >  1:10697523/1‑101 (MQ=255)
                                                            actaaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAAc                                          <  1:501539/101‑1 (MQ=255)
                                                                     gcgcCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATgg                                 <  1:13624338/101‑1 (MQ=255)
                                                                      cgcCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGa                                >  1:8323326/1‑101 (MQ=255)
                                                                      cgcCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATAGa                                <  1:14376653/101‑1 (MQ=255)
                                                                       gcCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGaa                               <  1:15095529/101‑1 (MQ=255)
                                                                        ccATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAg                              <  1:14997728/101‑1 (MQ=255)
                                                                          atatGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGgc                            >  1:10044063/1‑101 (MQ=255)
                                                                             tGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCt                         >  1:11666929/1‑101 (MQ=255)
                                                                             tGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCt                         <  1:13443471/101‑1 (MQ=255)
                                                                             tGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCt                         <  1:6083945/101‑1 (MQ=255)
                                                                             tGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCt                         >  1:3628111/1‑101 (MQ=255)
                                                                              gggTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTTGAACGTCTGATGGAAGGCGCTa                        <  1:10461241/101‑1 (MQ=255)
                                                                              gggTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTa                        >  1:14188971/1‑101 (MQ=255)
                                                                               ggTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTaa                       >  1:10399095/1‑101 (MQ=255)
                                                                                 ttGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAgc                     <  1:10051615/101‑1 (MQ=255)
                                                                                 ttGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAgc                     <  1:15170834/101‑1 (MQ=255)
                                                                                  tGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAgct                    <  1:704448/101‑1 (MQ=255)
                                                                                   ggACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAgctg                   <  1:13699234/101‑1 (MQ=255)
                                                                                     aCTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAgctgct                 <  1:471403/101‑1 (MQ=255)
                                                                                     aCTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAgctgct                 <  1:1543098/101‑1 (MQ=255)
                                                                                        gAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGGCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGaa              <  1:13187232/101‑1 (MQ=255)
                                                                                        gAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGaa              <  1:12188547/101‑1 (MQ=255)
                                                                                         aaGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAAAGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGaaa             >  1:9009860/1‑101 (MQ=255)
                                                                                           gCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATg           >  1:1717085/1‑101 (MQ=255)
                                                                                            cGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGt          <  1:290493/101‑1 (MQ=255)
                                                                                              gCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTc        >  1:8495028/1‑101 (MQ=255)
                                                                                               cTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCa       >  1:5439696/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
GTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTT  >  NC_007779/1212926‑1213126

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.