Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,213,028 G→A 2.5% intergenic (+272/+229) 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,0280G→A 2.5% 9.0 ‑0.8 ‑0.1274intergenic (+272/+229)mcrA/elbA5‑methylcytosine‑specific restriction endonuclease B/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (3/4):  ref base (153/114):  total (156/118)
Fisher's exact test strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.55e-01

Only 100 of 281 total aligned reads displayed.

TATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTG  >  NC_007779/1212928‑1213128
                                                                                                    |                                                                                                    
tataTGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTg                                                                                                      >  1:13510029/1‑101 (MQ=255)
 atatGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGa                                                                                                     >  1:972886/1‑101 (MQ=255)
 atatGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGa                                                                                                     >  1:8030111/1‑101 (MQ=255)
  tatGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGAt                                                                                                    <  1:7264673/101‑1 (MQ=255)
   atGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGAtt                                                                                                   <  1:14441130/101‑1 (MQ=255)
   atGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGAtt                                                                                                   <  1:14170476/101‑1 (MQ=255)
    tGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTg                                                                                                  <  1:1501414/101‑1 (MQ=255)
       aaaTTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTa                                                                                               >  1:2995311/1‑101 (MQ=255)
        aaTTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaa                                                                                              >  1:8453825/1‑101 (MQ=255)
          ttAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAg                                                                                            >  1:11657193/1‑101 (MQ=255)
             gTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                                                         <  1:14138670/101‑1 (MQ=255)
             gTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTa                                                                                         <  1:11688674/101‑1 (MQ=255)
              tAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTAt                                                                                        <  1:9209893/101‑1 (MQ=255)
               aaaaTTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATg                                                                                       <  1:4171658/101‑1 (MQ=255)
                 aaTTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGCTTGTTAAAGGTATGGa                                                                                     >  1:6160884/1‑101 (MQ=255)
                    tATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAgg                                                                                  <  1:7404486/101‑1 (MQ=255)
                     aTGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgc                                                                                 >  1:12735517/1‑101 (MQ=255)
                     aTGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgc                                                                                 <  1:11865288/101‑1 (MQ=18)
                      tGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcg                                                                                <  1:8412142/101‑1 (MQ=255)
                       gTTGCTATGTCAGATAGTTACGATTTAGTCGTCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcgc                                                                               >  1:13905142/1‑101 (MQ=255)
                       gTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcgc                                                                               >  1:13653961/1‑101 (MQ=255)
                        ttGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGca                                                                              >  1:7778375/1‑101 (MQ=255)
                         tGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaa                                                                             >  1:14462950/1‑101 (MQ=255)
                         tGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaa                                                                             >  1:12080110/1‑101 (MQ=255)
                         tGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaa                                                                             >  1:6073461/1‑101 (MQ=255)
                         tGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGcaa                                                                             <  1:3485376/101‑1 (MQ=17)
                           cTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatc                                                                           >  1:633124/1‑101 (MQ=255)
                            tATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatca                                                                          >  1:620385/1‑101 (MQ=255)
                            tATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatca                                                                          <  1:2151107/101‑1 (MQ=255)
                            tATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatca                                                                          >  1:3732614/1‑101 (MQ=255)
                            tATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatca                                                                          >  1:10913915/1‑101 (MQ=255)
                              tGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGcaatcaat                                                                        >  1:7201848/1‑101 (MQ=255)
                               gtaaGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATg                                                                       >  1:11182659/4‑101 (MQ=255)
                               gTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATg                                                                       >  1:2303570/1‑101 (MQ=255)
                               gTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATg                                                                       >  1:10481972/1‑101 (MQ=255)
                                tCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGc                                                                      <  1:167017/101‑1 (MQ=255)
                                tCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGc                                                                      >  1:13219317/1‑101 (MQ=255)
                                   gATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCaa                                                                   >  1:2225233/1‑101 (MQ=255)
                                     tAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGa                                                                 >  1:10226266/1‑101 (MQ=255)
                                     tAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGa                                                                 >  1:3977198/1‑101 (MQ=255)
                                       gTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGAcc                                                               <  1:11793379/101‑1 (MQ=255)
                                       gTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTCACCTGATTGTTAAAGGTATGGAAGGCGCAATCAAGGCCAAGAcc                                                               >  1:10596483/1‑101 (MQ=255)
                                        ttACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCg                                                              >  1:4082795/1‑101 (MQ=255)
                                        ttACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCg                                                              <  1:8884391/101‑1 (MQ=255)
                                        ttACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGAAATAAATGCCAAGAACg                                                              >  1:5788914/1‑101 (MQ=255)
                                         tACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTCAAGGTATGGAAGGCGCAATCAATGCCAAGACCGt                                                             >  1:9130260/1‑101 (MQ=255)
                                           cGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTaa                                                           >  1:10154153/1‑101 (MQ=255)
                                           cGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTaa                                                           <  1:11408380/101‑1 (MQ=255)
                                             aTTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACt                                                         >  1:9857588/1‑101 (MQ=255)
                                             aTTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACt                                                         <  1:12195233/101‑1 (MQ=255)
                                              tttAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACtt                                                        <  1:14041581/101‑1 (MQ=255)
                                              tttAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACtt                                                        >  1:10772016/1‑101 (MQ=255)
                                               ttAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTa                                                       <  1:7823674/101‑1 (MQ=255)
                                                tAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTAt                                                      >  1:2450174/1‑101 (MQ=255)
                                                tAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTAt                                                      <  1:5467258/101‑1 (MQ=255)
                                                 aGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATg                                                     >  1:6516021/1‑101 (MQ=255)
                                                   tCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGAc                                                   <  1:5003447/101‑1 (MQ=255)
                                                      tCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTc                                                <  1:9509940/101‑1 (MQ=255)
                                                      tCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTc                                                >  1:10866631/1‑101 (MQ=255)
                                                      tCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTc                                                <  1:12548575/101‑1 (MQ=255)
                                                      tCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTc                                                <  1:4492099/101‑1 (MQ=255)
                                                         aactaaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaa                                             >  1:6771425/1‑101 (MQ=255)
                                                         aactaaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaa                                             >  1:2889561/1‑101 (MQ=255)
                                                          actaaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAAc                                            >  1:8225731/1‑101 (MQ=255)
                                                          actaaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAAc                                            >  1:2782480/1‑101 (MQ=255)
                                                            taaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGt                                          >  1:1731906/1‑101 (MQ=255)
                                                             aaTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTc                                         <  1:1111409/101‑1 (MQ=255)
                                                              aTGCTGCGCCATATGAGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCt                                        >  1:10994327/1‑101 (MQ=255)
                                                               tgctgcGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTTGAACGTCTg                                       >  1:3234336/1‑101 (MQ=255)
                                                               tgctgcGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTg                                       >  1:9053096/1‑101 (MQ=255)
                                                               tgctgcGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTg                                       >  1:521836/1‑101 (MQ=255)
                                                                gctgcGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGa                                      >  1:2044443/1‑101 (MQ=255)
                                                                 ctgcGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGAt                                     <  1:12463136/101‑1 (MQ=255)
                                                                  tgcGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATg                                    >  1:11429516/1‑101 (MQ=255)
                                                                     gcCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGTaa                                 <  1:13770868/101‑1 (MQ=255)
                                                                         tatGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGgcg                             <  1:5887555/101‑1 (MQ=255)
                                                                         tatGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGgcg                             <  1:7847704/101‑1 (MQ=255)
                                                                          atGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGgcgc                            >  1:14319738/1‑101 (MQ=255)
                                                                          atGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGgcgc                            >  1:6937227/1‑101 (MQ=255)
                                                                            gggTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTa                          >  1:1301012/1‑101 (MQ=255)
                                                                             ggTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTaa                         >  1:11205914/1‑101 (MQ=255)
                                                                              gTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAg                        >  1:433522/1‑101 (MQ=255)
                                                                                 ggACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAgctg                     >  1:10033269/1‑101 (MQ=255)
                                                                                  gACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCCATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTCAgctgc                    >  1:3128082/1‑101 (MQ=255)
                                                                                   aCTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAgctgct                   <  1:860828/101‑1 (MQ=255)
                                                                                     tGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGa                 >  1:4609762/1‑101 (MQ=255)
                                                                                        aGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAAt              <  1:4453367/101‑1 (MQ=255)
                                                                                        aGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAAt              <  1:10734355/101‑1 (MQ=255)
                                                                                        aGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAAt              <  1:3352581/101‑1 (MQ=255)
                                                                                           ggCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGtt           >  1:8500252/1‑101 (MQ=255)
                                                                                           ggCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGtt           >  1:3584651/1‑101 (MQ=255)
                                                                                           ggCTGACCTGATTGGTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGtt           <  1:10536820/101‑1 (MQ=255)
                                                                                               gACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCaga       >  1:4124038/1‑101 (MQ=255)
                                                                                                 ccTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGt     <  1:12863857/101‑1 (MQ=255)
                                                                                                  cTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGtt    >  1:5376361/1‑101 (MQ=255)
                                                                                                  cTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGtt    >  1:13340/1‑101 (MQ=255)
                                                                                                   tGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGttt   <  1:12311773/101‑1 (MQ=255)
                                                                                                   tGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGttt   >  1:6985618/1‑101 (MQ=255)
                                                                                                   tGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGttt   <  1:865591/101‑1 (MQ=255)
                                                                                                   tGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGttt   >  1:13176390/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
TATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTG  >  NC_007779/1212928‑1213128

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