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.5% 659.4 / NA 204Q315R (CAG→CGG) acrDAcrD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (80/123);  minor base T (1/0);  total (81/123)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

CAGCCTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGG  >  NC_000913/2586468‑2586650
                                                                                          |                                                                                            
cagccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGt                                                                                            >  1:5464586/1‑93 (MQ=255)
 agccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTa                                                                                           >  1:4698724/1‑93 (MQ=255)
  gccTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAt                                                                                          <  2:407141/93‑1 (MQ=255)
    cTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAttt                                                                                        >  2:6241575/1‑93 (MQ=255)
    cTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTAttt                                                                                        <  2:2650358/93‑1 (MQ=255)
     tCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTc                                                                                       >  1:3890419/1‑93 (MQ=255)
     tCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTc                                                                                       >  1:3186284/1‑93 (MQ=255)
     tCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTc                                                                                       >  2:1905794/1‑93 (MQ=255)
     tCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTc                                                                                       <  1:2480834/93‑1 (MQ=255)
     tCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTc                                                                                       >  2:4031808/1‑93 (MQ=255)
     tCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGAGCGGTATTTc                                                                                       <  1:6279205/93‑1 (MQ=255)
      ccGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTcc                                                                                      <  1:2093093/93‑1 (MQ=255)
      ccGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTcc                                                                                      <  2:5211788/93‑1 (MQ=255)
            tGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATg                                                                                >  2:3004714/1‑93 (MQ=255)
             gggggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATgg                                                                               >  2:2132485/1‑93 (MQ=255)
              ggggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGa                                                                              <  2:2530216/93‑1 (MQ=255)
              ggggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGa                                                                              >  2:4720496/1‑93 (MQ=255)
               gggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGAc                                                                             >  1:2926810/1‑93 (MQ=255)
                ggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACt                                                                            <  2:2422358/93‑1 (MQ=255)
                ggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACt                                                                            <  1:3382011/93‑1 (MQ=255)
                ggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACt                                                                            <  2:741804/93‑1 (MQ=255)
                ggTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACt                                                                            <  2:4805332/93‑1 (MQ=255)
                 gTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTg                                                                           >  2:5399864/1‑93 (MQ=255)
                   aaaaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGa                                                                         <  2:3231891/93‑1 (MQ=255)
                    aaaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGaa                                                                        <  2:1238955/93‑1 (MQ=255)
                    aaaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGaa                                                                        >  1:429699/1‑93 (MQ=255)
                     aaCTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAAt                                                                       <  1:1115839/93‑1 (MQ=255)
                       cTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATAc                                                                     <  2:599475/93‑1 (MQ=255)
                        tGGACTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACa                                                                    >  2:4910951/1‑93 (MQ=255)
                            ccccGGCGCTAACGAAATGGCGACAGCTGAGCTGGTGCTAATTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGt                                                                <  2:437830/91‑1 (MQ=255)
                             tCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTg                                                               >  2:1301355/1‑93 (MQ=255)
                             tCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTg                                                               <  1:5669336/93‑1 (MQ=255)
                              ccGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTgg                                                              <  2:1067465/93‑1 (MQ=255)
                               cGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGc                                                             <  1:55108/93‑1 (MQ=255)
                                 gcgcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGt                                                           >  1:5724470/1‑93 (MQ=255)
                                  cgcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTa                                                          >  1:2618804/1‑93 (MQ=255)
                                   gcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTAt                                                         <  2:3904774/93‑1 (MQ=255)
                                   gcTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTAt                                                         >  1:1777658/1‑93 (MQ=255)
                                     tAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGa                                                       >  2:3932573/1‑93 (MQ=255)
                                     tAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGa                                                       <  2:4977759/93‑1 (MQ=255)
                                      aaCGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGaa                                                      >  2:6127668/1‑93 (MQ=255)
                                      aaCGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGaa                                                      <  2:1180618/93‑1 (MQ=255)
                                      aaCGAAATAGCGACAGCGGAGCTGGAGCTCAATCGTTTCGACCCGCCGGCGCGCGCATTCCCGCATGGACTCGGACACAAGGTGGCGTATGaa                                                      >  2:209481/1‑93 (MQ=255)
                                       aCGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGaaa                                                     >  2:2990440/1‑93 (MQ=255)
                                        cGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAac                                                    >  1:158846/1‑93 (MQ=255)
                                        cGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAac                                                    <  1:5881616/93‑1 (MQ=255)
                                         gAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAacc                                                   <  1:431161/93‑1 (MQ=255)
                                          aaaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAacca                                                  >  1:927647/1‑93 (MQ=255)
                                          aaaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAacca                                                  >  1:4526690/1‑93 (MQ=255)
                                          aaaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAacca                                                  <  1:363122/93‑1 (MQ=255)
                                          aaaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAacca                                                  >  1:2729381/1‑93 (MQ=255)
                                           aaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAaccac                                                 <  1:70652/93‑1 (MQ=255)
                                           aaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAaccac                                                 <  1:6136434/93‑1 (MQ=255)
                                           aaTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAaccac                                                 <  1:1000348/93‑1 (MQ=255)
                                            aTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAaccacc                                                <  2:2395503/93‑1 (MQ=255)
                                            aTGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAaccacc                                                <  1:5831817/93‑1 (MQ=255)
                                             tGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCt                                               >  2:3425245/1‑93 (MQ=255)
                                             tGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCt                                               <  2:4593608/93‑1 (MQ=255)
                                              ggCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTc                                              <  2:2850511/93‑1 (MQ=255)
                                                 gACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGtt                                           >  2:5558791/1‑93 (MQ=255)
                                                  aCAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGttt                                          <  2:3645585/93‑1 (MQ=255)
                                                      cGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTa                                      >  1:2022177/1‑93 (MQ=255)
                                                        gAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTaaa                                    >  1:2379364/1‑93 (MQ=255)
                                                          gCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGc                                  <  1:3529439/93‑1 (MQ=255)
                                                          gCTGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGc                                  <  2:2858782/93‑1 (MQ=255)
                                                            tGGTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCt                                >  2:1542038/1‑93 (MQ=255)
                                                              gTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTcc                              <  2:4353382/93‑1 (MQ=255)
                                                              gTGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTcc                              <  2:1360494/93‑1 (MQ=255)
                                                               tGCTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATCCAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCa                             <  2:3300430/93‑1 (MQ=255)
                                                                gCTCAATCGTCTCGTGGAGCTGGCGCGGTATTACCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCAt                            <  2:1995044/93‑1 (MQ=255)
                                                                 cTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCAtt                           <  1:2192192/93‑1 (MQ=255)
                                                                 cTCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCAtt                           <  2:5722375/93‑1 (MQ=255)
                                                                 cTCAATCGTCTCGACGAGCCAGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCAtt                           >  1:4374092/1‑93 (MQ=255)
                                                                  tCAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTg                          >  1:1319767/1‑93 (MQ=255)
                                                                   cAATCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGa                         <  2:2751490/93‑1 (MQ=255)
                                                                    aaTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGaa                        >  2:3760764/1‑93 (MQ=255)
                                                                    aaTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGaa                        <  1:821679/93‑1 (MQ=255)
                                                                    aaTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGaa                        >  1:4581570/1‑93 (MQ=255)
                                                                     aTCGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAg                       <  1:314713/93‑1 (MQ=255)
                                                                       cGTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGAc                     >  2:2286420/1‑93 (MQ=255)
                                                                        gTCTCGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACg                    <  1:1829901/93‑1 (MQ=255)
                                                                         tctcGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgt                   <  2:5154822/93‑1 (MQ=255)
                                                                         tctcGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgt                   <  2:2269101/93‑1 (MQ=255)
                                                                          ctcGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtg                  <  1:6236749/93‑1 (MQ=255)
                                                                           tcGACGAGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtgg                 <  1:283241/93‑1 (MQ=255)
                                                                             gacgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACgtggtg               <  1:4415671/93‑1 (MQ=255)
                                                                              acgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGa              >  2:5377544/1‑93 (MQ=255)
                                                                              acgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGa              >  1:141358/1‑93 (MQ=255)
                                                                               cgaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaa             <  2:6173626/93‑1 (MQ=255)
                                                                                gaGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaaa            >  1:374047/1‑93 (MQ=255)
                                                                                 aGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaaaa           <  1:5016069/93‑1 (MQ=255)
                                                                                 aGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaaaa           >  1:5140716/1‑93 (MQ=255)
                                                                                 aGCTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGaaaa           <  1:2560406/93‑1 (MQ=255)
                                                                                  gCTGGCGCTGTATTTCCCGCATGGACTGGTATACAAGGTGGCGTATGAAATCACCTTGTTTGTTAAAGCCTCCATTGTACACGTGGTGTAAAc          >  2:5471805/1‑93 (MQ=255)
                                                                                   cTGGCGCGGTATTTCCCGCATGGACTGGAATACAAGCTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACg         >  2:2172723/1‑93 (MQ=255)
                                                                                    tGGCGCGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgc        >  2:3021509/1‑93 (MQ=255)
                                                                                      gcgcGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAACCCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgctg      <  2:3127295/93‑1 (MQ=255)
                                                                                      gcgcGGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgctg      <  1:4989315/93‑1 (MQ=255)
                                                                                       cgcgGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACgctgc     <  2:767150/93‑1 (MQ=255)
                                                                                          aGTAATTCCCGCATGGGCTGGAACACAAGGTGGCGGATGAAACCCCCTCGTTTGTTAAAGACTCCATTGAAGACGTGGTGAACACGCTGCTgg  >  2:2481420/1‑93 (MQ=255)
                                                                                          |                                                                                            
CAGCCTCCGGGCTGGGGGTAAAACTGGCCTCCGGCGCTAACGAAATGGCGACAGCGGAGCTGGTGCTCAATCGTCTCGACGAGCTGGCGCAGTATTTCCCGCATGGACTGGAATACAAGGTGGCGTATGAAACCACCTCGTTTGTTAAAGCCTCCATTGAAGACGTGGTGAAAACGCTGCTGG  >  NC_000913/2586468‑2586650

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