Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,586,558 A→G Q315R (CAG→CGG)  acrD → AcrD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,586,5580AG99.1% 683.2 / ‑3.3 213Q315R (CAG→CGG) acrDAcrD
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (99/112);  total (100/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 221 total aligned reads displayed.

CCAGCCTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGC  >  NC_000913/2586467‑2586647
                                                                                           |                                                                                         
  cgccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTa                                                                                        <  2:379577/92‑1 (MQ=255)
   gccTCCTGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAt                                                                                       <  1:4562558/93‑1 (MQ=255)
   gccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAt                                                                                       >  2:9884605/1‑93 (MQ=255)
   gccTCCGGGCTGGGGGTAAAACTGGCCGCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAt                                                                                       <  1:6037954/93‑1 (MQ=255)
    ccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAtt                                                                                      >  1:3098043/1‑93 (MQ=255)
     cTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTTAATCGTCTCGACGAGCTGGCGCGGTAttt                                                                                     >  1:5496501/1‑93 (MQ=255)
      tCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTc                                                                                    >  1:1063559/1‑93 (MQ=255)
      tCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTc                                                                                    >  1:6443541/1‑93 (MQ=255)
      tCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTc                                                                                    <  2:8382761/93‑1 (MQ=255)
         gggCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCg                                                                                 >  2:8135971/1‑93 (MQ=255)
          ggCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGc                                                                                >  1:9541598/1‑93 (MQ=255)
          ggCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGc                                                                                <  1:4358046/93‑1 (MQ=255)
           gCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCa                                                                               >  2:7973996/1‑93 (MQ=255)
           gCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCa                                                                               <  1:4505258/93‑1 (MQ=255)
           gCTGGGGGTAAAACTGGCCTCCGGCACTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCa                                                                               <  2:878726/93‑1 (MQ=255)
             tGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATg                                                                             >  2:6233968/1‑93 (MQ=255)
               gggaTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCAGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGa                                                                           <  1:8926022/93‑1 (MQ=255)
                gggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGAc                                                                          <  1:9389098/93‑1 (MQ=255)
                      aaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAAt                                                                    >  1:3617886/1‑93 (MQ=255)
                      aaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAAt                                                                    >  1:9557728/1‑93 (MQ=255)
                       aCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATa                                                                   <  1:8016012/93‑1 (MQ=255)
                        cTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATAc                                                                  <  2:779902/93‑1 (MQ=255)
                        cTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATAc                                                                  <  2:846015/93‑1 (MQ=255)
                         tGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACa                                                                 <  1:10021902/93‑1 (MQ=255)
                         tGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACa                                                                 >  2:887634/1‑93 (MQ=255)
                           gCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAg                                                               <  1:3410820/93‑1 (MQ=255)
                           gCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAg                                                               <  2:4625279/93‑1 (MQ=255)
                           gCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAg                                                               <  2:2884608/93‑1 (MQ=255)
                            ccTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAgg                                                              <  2:91046/93‑1 (MQ=255)
                            ccTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAgg                                                              <  2:4764655/93‑1 (MQ=255)
                            ccTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAgg                                                              <  1:9267278/93‑1 (MQ=255)
                             cTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCGCGCATGGACTGGAACACAAGGt                                                             <  2:7470327/93‑1 (MQ=255)
                             cTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGt                                                             >  2:6065185/1‑93 (MQ=255)
                             cTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGt                                                             <  1:2173847/93‑1 (MQ=255)
                              tCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTg                                                            >  2:9571359/1‑93 (MQ=255)
                              tCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTg                                                            <  2:1828465/93‑1 (MQ=255)
                              tCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTg                                                            <  1:7547164/93‑1 (MQ=255)
                              tCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTg                                                            >  1:2487992/1‑93 (MQ=255)
                               ccGGCGCTAACGAAATGGCGACCGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTgg                                                           >  1:9638585/1‑93 (MQ=255)
                               ccGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTgg                                                           <  2:3643290/93‑1 (MQ=255)
                                cGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGc                                                          <  2:2606007/93‑1 (MQ=255)
                                 ggCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCg                                                         <  1:7242673/93‑1 (MQ=255)
                                   cgcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTa                                                       <  2:6885631/93‑1 (MQ=255)
                                    gcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTAt                                                      >  2:4965602/1‑93 (MQ=255)
                                    gcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTAt                                                      >  2:2954390/1‑93 (MQ=255)
                                    gcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACCAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTAt                                                      >  2:7014561/1‑93 (MQ=255)
                                     cTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATg                                                     >  2:9507502/1‑93 (MQ=255)
                                      tAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGa                                                    <  1:8274953/93‑1 (MQ=255)
                                      tAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGa                                                    <  2:7533008/93‑1 (MQ=255)
                                      tAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGa                                                    >  1:2717838/1‑93 (MQ=255)
                                         cGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACGAGGTGGCGTATGAAac                                                 >  1:3674519/1‑93 (MQ=255)
                                         cGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAac                                                 <  2:2428895/93‑1 (MQ=255)
                                          gAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAacc                                                <  1:3537713/93‑1 (MQ=255)
                                           aaaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAacca                                               <  2:5466553/93‑1 (MQ=255)
                                           aaaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAacca                                               <  2:1748342/93‑1 (MQ=255)
                                           aaaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAacca                                               <  1:1324494/93‑1 (MQ=255)
                                            aaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAaccac                                              >  1:8359828/1‑93 (MQ=255)
                                            aaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAaccac                                              <  2:1409927/93‑1 (MQ=255)
                                               ggCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTc                                           >  2:9901688/1‑93 (MQ=255)
                                                  gACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATAGAAGGTGGCGTATGAAACCACCTCGtt                                        <  2:2860129/93‑1 (MQ=255)
                                                   aCAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGttt                                       >  2:7597466/1‑93 (MQ=255)
                                                    cAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTg                                      >  1:3714735/1‑93 (MQ=255)
                                                    cAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTg                                      <  2:6232772/93‑1 (MQ=255)
                                                     aGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGt                                     >  1:4602863/1‑93 (MQ=255)
                                                     aGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGt                                     >  1:8838681/1‑93 (MQ=255)
                                                           gCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGc                               >  2:3591052/1‑93 (MQ=255)
                                                           gCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGc                               <  1:7339887/93‑1 (MQ=255)
                                                           gCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGc                               >  1:984169/1‑93 (MQ=255)
                                                            cTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGcc                              >  1:2337559/1‑93 (MQ=255)
                                                            cTGGTGCTAAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTAGTTTGTTACAGcc                              >  1:4193361/1‑93 (MQ=255)
                                                             tGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCt                             >  1:6996849/1‑93 (MQ=255)
                                                              cgTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTc                            <  1:65819/92‑1 (MQ=255)
                                                                 gCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCAt                         >  1:3044393/1‑93 (MQ=255)
                                                                  cTCAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCAtt                        <  1:5472879/93‑1 (MQ=255)
                                                                   tCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTg                       >  1:8565876/1‑93 (MQ=255)
                                                                   tCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTg                       <  2:3266816/93‑1 (MQ=255)
                                                                    cAATCGTTTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCATCCATTGa                      >  2:2271561/1‑93 (MQ=255)
                                                                     aaTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGaa                     <  2:92585/93‑1 (MQ=255)
                                                                      aTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAg                    >  2:9437012/1‑93 (MQ=255)
                                                                      aTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAg                    <  2:6514272/93‑1 (MQ=255)
                                                                      aTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAg                    <  1:6916239/93‑1 (MQ=255)
                                                                      aTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAg                    >  2:3481968/1‑93 (MQ=255)
                                                                       tCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGa                   <  2:1654354/93‑1 (MQ=255)
                                                                       tCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGa                   >  1:6227773/1‑93 (MQ=255)
                                                                       tCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGa                   <  1:5303290/93‑1 (MQ=255)
                                                                       tCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGGTAAAGCCTCCATTGAAGa                   <  2:1183601/93‑1 (MQ=255)
                                                                          tctcGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgt                >  2:6747401/1‑93 (MQ=255)
                                                                            tcGAGGGGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGCCgtgg              <  2:1926576/93‑1 (MQ=255)
                                                                              gacgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtggtg            >  1:1569093/1‑93 (MQ=255)
                                                                               acgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGa           <  1:566142/93‑1 (MQ=255)
                                                                                cgaGCTGGCGCGGGATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtggtgca          >  2:1896119/1‑91 (MQ=255)
                                                                                 gaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAACCCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaaa         <  2:6676033/93‑1 (MQ=255)
                                                                                 gaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaaa         <  1:7791469/93‑1 (MQ=255)
                                                                                  aGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaaaa        <  2:8531732/93‑1 (MQ=255)
                                                                                    cTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACg      <  2:2245856/93‑1 (MQ=255)
                                                                                     tGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgc     <  1:6669150/93‑1 (MQ=255)
                                                                                     tGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgc     <  1:7852003/93‑1 (MQ=255)
                                                                                     tGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgc     >  2:5253837/1‑93 (MQ=255)
                                                                                     tGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgc     >  2:5110885/1‑93 (MQ=255)
                                                                                       gcgcGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgctg   <  2:9424690/93‑1 (MQ=255)
                                                                                           |                                                                                         
CCAGCCTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGC  >  NC_000913/2586467‑2586647

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: