Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,802,070 A→G 7.2% K689K (AAA→AAG)  nrdE → ribonucleoside‑diphosphate reductase 2 subunit alpha

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077792,802,0700A→G 7.2% 33.5 ‑0.2 ‑0.6208K689K (AAA→AAG) nrdEribonucleoside‑diphosphate reductase 2 subunit alpha
Reads supporting (aligned to +/- strand):  new base (8/7):  ref base (110/83):  total (118/90)
Fisher's exact test strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.64e-01

Only 100 of 219 total aligned reads displayed.

ATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTATGAAACTCTCAC  >  NC_007779/2801970‑2802170
                                                                                                    |                                                                                                    
aTGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCaaa                                                                                                      <  1:7402589/101‑1 (MQ=255)
aTGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGCATCaaa                                                                                                      >  1:5888523/1‑101 (MQ=255)
   tCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACg                                                                                                   >  1:10031013/1‑101 (MQ=255)
   tCGATCAGGGGCTGTCGCTGACGATTTTTTTCCCCGATACCGCCACCACTCGCGATATGAACAAAGCGCAAATTTAAGCCTGGCGCAAGGGCCTCaaaaac                                                                                                   >  1:216482/1‑99 (MQ=255)
    cGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCTCCACCACTAGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGc                                                                                                  >  1:13898123/1‑101 (MQ=255)
    cGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGc                                                                                                  <  1:4986391/101‑1 (MQ=255)
     gATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGct                                                                                                 >  1:11527465/1‑101 (MQ=255)
     gATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGct                                                                                                 >  1:10288265/1‑101 (MQ=255)
      aTCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGctc                                                                                                <  1:3569811/101‑1 (MQ=255)
      aTCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGctc                                                                                                >  1:8546569/1‑101 (MQ=255)
      aTCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGctc                                                                                                <  1:841048/101‑1 (MQ=255)
       tCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGctct                                                                                               <  1:3747381/101‑1 (MQ=255)
       tCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGctct                                                                                               <  1:5432824/101‑1 (MQ=255)
       tCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGctct                                                                                               >  1:7583158/1‑101 (MQ=255)
          ggggCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt                                                                                            <  1:4583075/101‑1 (MQ=255)
          ggggCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt                                                                                            <  1:976041/101‑1 (MQ=255)
          ggggCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt                                                                                            >  1:13349486/1‑101 (MQ=255)
          ggggCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt                                                                                            >  1:6167782/1‑101 (MQ=255)
               tGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATc                                                                                       >  1:2991253/1‑101 (MQ=255)
               tGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATc                                                                                       >  1:11729868/1‑101 (MQ=255)
                  cGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGc                                                                                    <  1:9246740/101‑1 (MQ=255)
                  cGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGc                                                                                    >  1:10932241/1‑101 (MQ=255)
                  cGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGc                                                                                    >  1:9002158/1‑101 (MQ=255)
                   gCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGcc                                                                                   >  1:12655387/1‑101 (MQ=255)
                     tGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTg                                                                                 <  1:5744840/101‑1 (MQ=255)
                        cGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGt                                                                              <  1:6591667/101‑1 (MQ=255)
                        cGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGt                                                                              <  1:6505256/101‑1 (MQ=255)
                        cGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGt                                                                              <  1:6044021/101‑1 (MQ=255)
                         gCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTc                                                                             >  1:10662165/1‑101 (MQ=255)
                           tttttttCCCCGATACCGCCACCTCGCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAg                                                                           <  1:11611284/101‑1 (MQ=255)
                           tttttttCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAg                                                                           <  1:469099/101‑1 (MQ=255)
                            ttttttCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGa                                                                          <  1:1149223/101‑1 (MQ=255)
                                ttCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGgc                                                                      <  1:4546196/101‑1 (MQ=255)
                                ttCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGgc                                                                      >  1:11130702/1‑101 (MQ=255)
                                 tCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGgcg                                                                     >  1:7826344/1‑101 (MQ=255)
                                     cGATACCGCCACCGCTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTgg                                                                 <  1:10424499/101‑1 (MQ=255)
                                     cGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTgg                                                                 <  1:6378625/101‑1 (MQ=255)
                                      gATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGa                                                                >  1:4769476/1‑101 (MQ=255)
                                      gATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGa                                                                >  1:4783723/1‑101 (MQ=255)
                                       aTACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGaa                                                               >  1:8309592/1‑101 (MQ=255)
                                       aTACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGaa                                                               >  1:1252192/1‑101 (MQ=255)
                                        tACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAg                                                              >  1:9959054/1‑101 (MQ=255)
                                        tACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAg                                                              <  1:1816983/101‑1 (MQ=255)
                                         aCCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAgg                                                             >  1:13655901/1‑101 (MQ=255)
                                           cGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCa                                                           >  1:12046518/1‑101 (MQ=255)
                                            gCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCAc                                                          <  1:2579799/101‑1 (MQ=255)
                                               accacTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGa                                                       <  1:4301736/101‑1 (MQ=255)
                                                gcacTCGTGACCTCAACCGTGCGCAGATTTACGCATGGCGCAAGGGTATTAAGACTCTGTACTACATCCGCCTGCGTCAGGTTGCGCTGGAGGGTACTGag                                                      <  1:2064650/100‑2 (MQ=255)
                                                gcacTCGTGACCTCAACCGTGCGCAGATTTACGCATGGCGCAAGGGTATTAAGACTCTGTACTACATCCGCCTGCGTCAGGTTGCGCTGGAGGGTACTGag                                                      <  1:6498344/100‑2 (MQ=255)
                                                      gcgATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGaa                                                >  1:7128318/1‑101 (MQ=255)
                                                           atCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTg                                           >  1:12500858/1‑101 (MQ=255)
                                                           atCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTg                                           <  1:5677897/101‑1 (MQ=255)
                                                           atCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTg                                           <  1:9453299/101‑1 (MQ=255)
                                                            tCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGc                                          <  1:8916588/101‑1 (MQ=255)
                                                              aaccgtgcgcAGATTTACGCATGGCGCAAGGGTATTAAGACTCTGTACTACATCCGCCTGCGTCAGGTTGCGCTGGAGGGTACTGAGGTCGAGGGTTGCGt                                        >  1:3439310/7‑101 (MQ=255)
                                                              aacaaAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGt                                        <  1:2275302/101‑1 (MQ=255)
                                                              aacaaAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGt                                        <  1:14323834/101‑1 (MQ=255)
                                                               acaaAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGtc                                       <  1:9724733/101‑1 (MQ=255)
                                                               acaaAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGtc                                       >  1:9643782/1‑101 (MQ=255)
                                                                 aaaGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGtctc                                     >  1:2807237/1‑101 (MQ=255)
                                                                  aaGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTcc                                    >  1:11449413/1‑101 (MQ=255)
                                                                  aaGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTcc                                    >  1:4059518/1‑101 (MQ=255)
                                                                   aGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCt                                   <  1:1919889/101‑1 (MQ=255)
                                                                   aGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCt                                   >  1:445231/1‑101 (MQ=255)
                                                                   aGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCt                                   >  1:4411270/1‑101 (MQ=255)
                                                                    gcgcAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCtg                                  <  1:409720/101‑1 (MQ=255)
                                                                      gcAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCtgtg                                >  1:458464/1‑101 (MQ=255)
                                                                      gcAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCtgtg                                >  1:106205/1‑101 (MQ=255)
                                                                        aGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCa                              <  1:13284992/101‑1 (MQ=255)
                                                                        aGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCa                              <  1:13089279/101‑1 (MQ=255)
                                                                          aTTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCCCTGAAATTGAAGGCTGCGTCTACTGTGCACt                            <  1:13527260/101‑1 (MQ=255)
                                                                            ttACGCATGGCGCAAGGGTATTAAGACTCTGTACTACATCCGCCTGCGTCAGGTTGCGCTGGAGGGTACTGAGGTCGAGGGTTGCGTATCCtgcatgcttt                          >  1:11916386/1‑93 (MQ=255)
                                                                             tACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTa                         >  1:11189369/1‑101 (MQ=255)
                                                                             tACGCATGGCGCAAGGGTATTAAGACTCTGTACTACATCCGCCTGCGTCAGGTTGCGCTGGAGGGTACTGAGGTCGAGGGTTGCGTATCCtgcatgcttta                         >  1:7277460/1‑92 (MQ=255)
                                                                               cGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAg                       <  1:13983007/101‑1 (MQ=255)
                                                                               cGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAg                       >  1:8201023/1‑101 (MQ=255)
                                                                               cGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAg                       >  1:12029269/1‑101 (MQ=255)
                                                                                gCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAgg                      >  1:13366014/1‑101 (MQ=255)
                                                                                cccTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAgg                      >  1:13240742/2‑101 (MQ=255)
                                                                                  cTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGaa                    <  1:5861340/101‑1 (MQ=255)
                                                                                   tGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAat                   >  1:12964509/1‑101 (MQ=255)
                                                                                   tGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAat                   >  1:7477891/1‑101 (MQ=255)
                                                                                   tGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAat                   >  1:10852244/1‑101 (MQ=255)
                                                                                   tGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAat                   >  1:5133093/1‑101 (MQ=255)
                                                                                    ggCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAata                  <  1:388628/101‑1 (MQ=255)
                                                                                       gcAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCt               >  1:9132676/1‑101 (MQ=255)
                                                                                         aaGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTAt             >  1:10180392/1‑101 (MQ=255)
                                                                                         aaGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTAt             >  1:3682272/1‑101 (MQ=255)
                                                                                         aaGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTAt             <  1:9062119/101‑1 (MQ=255)
                                                                                         aaGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTAt             >  1:5237555/1‑101 (MQ=255)
                                                                                          aGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTATg            >  1:5216234/1‑101 (MQ=255)
                                                                                          aGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTATg            >  1:3615615/1‑101 (MQ=255)
                                                                                           gggTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTATGa           >  1:4441924/1‑101 (MQ=255)
                                                                                              tATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTATGAAAc        >  1:11851162/1‑101 (MQ=255)
                                                                                              tATCAAAACGCTATATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTATGAAAc        <  1:12660127/101‑1 (MQ=255)
                                                                                               aTCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTATGAAAct       <  1:11994090/101‑1 (MQ=255)
                                                                                                 cAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTATGAAActct     <  1:11635643/101‑1 (MQ=255)
                                                                                                 cAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTATGAAActct     >  1:4750767/1‑101 (MQ=255)
                                                                                                 cAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTATGAAActct     <  1:7320649/101‑1 (MQ=255)
                                                                                                  aaaaCGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTATGAAACtctc    >  1:1687932/1‑101 (MQ=255)
                                                                                                    |                                                                                                    
ATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGCGTCTCCTGTGCACTTTAAGGAATATCTATGAAACTCTCAC  >  NC_007779/2801970‑2802170

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.