Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,377,792 A→T 16.4% intergenic (‑70/+19) blc ← / ← ampC outer membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,377,7920AT16.4% 508.8 / 88.5 207intergenic (‑70/+19)blc/ampCouter membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (93/80);  new base T (10/24);  total (103/104)
Fisher's exact test for biased strand distribution p-value = 1.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-02

Only 100 of 213 total aligned reads displayed.

CGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGAT  >  NC_000913/4377693‑4377891
                                                                                                    |                                                                                                   
cGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGaaaa                                                                                                     <  1:699746/100‑1 (MQ=255)
 gCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAAt                                                                                                    >  1:598558/1‑100 (MQ=255)
 gCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAAt                                                                                                    >  2:280727/1‑100 (MQ=255)
  cTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAAtt                                                                                                   <  2:1510817/100‑1 (MQ=255)
   tGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTc                                                                                                  <  1:1487296/100‑1 (MQ=255)
     cGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCgg                                                                                                >  1:4366256/1‑100 (MQ=255)
      ggCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGa                                                                                               >  2:2849416/1‑100 (MQ=255)
        cAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAcc                                                                                             >  2:1168066/1‑100 (MQ=255)
        cAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAcc                                                                                             <  2:2870544/100‑1 (MQ=255)
         aaCGAGAGGGGGCAGGCGCATAAATGTTTCCTTGCTGGTTTTTTTCTAAGTGTGGATGACAGCAAGGAAGAGCGGGGAAAGGGTCCGAAGTTTCGGAcct                                                                                            <  1:1030902/100‑2 (MQ=255)
         aaCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcct                                                                                            <  1:232539/100‑2 (MQ=255)
         aaCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcct                                                                                            <  1:2640247/100‑2 (MQ=255)
         aaCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcct                                                                                            <  2:662811/100‑2 (MQ=255)
         aaCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAccc                                                                                            >  2:395486/1‑100 (MQ=255)
          aCGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcctt                                                                                           <  2:2886402/100‑3 (MQ=255)
           cGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttt                                                                                          <  2:2652641/100‑4 (MQ=255)
           cGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGa                                                                                          >  1:1735195/1‑100 (MQ=255)
           cGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGa                                                                                          >  2:2695132/1‑100 (MQ=255)
             agagGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttc                                                                                        >  1:1612376/1‑95 (MQ=255)
              gagGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATgg                                                                                       >  1:4341770/1‑100 (MQ=255)
              gagGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATAACAGCCAGGAAAAGCGTAGAAAAGGTCCGAAAATTAGGAACCGa                                                                                          >  1:3823492/1‑97 (MQ=255)
               agGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttctc                                                                                      <  1:3163439/100‑8 (MQ=255)
                gggAGCAGGCGCATAAATGTTTCCTTGCTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttctc                                                                                      <  2:3789448/99‑8 (MQ=38)
                gggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttctcc                                                                                     <  1:2676356/100‑9 (MQ=38)
                gggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGaa                                                                                     >  2:4262854/1‑100 (MQ=255)
                gggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGaa                                                                                     >  2:974970/1‑100 (MQ=255)
                 ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAt                                                                                    >  2:3566907/1‑100 (MQ=255)
                  gAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGACAAGCGGAGAAAAAGACCGAAAATTCGGACCCGATGGAAtt                                                                                   >  1:1340101/1‑100 (MQ=255)
                    gCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAtttt                                                                                 <  2:1096099/100‑1 (MQ=255)
                      aGGCGGATAAATGTTTCCTTACGGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGTCCCGATGGAATTTTAc                                                                               <  2:3087179/100‑1 (MQ=255)
                      aGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAATGCGGAGAAAATGTCCGAAAATTCGGACCGGATGGAAttttgc                                                                               <  2:2981305/100‑3 (MQ=255)
                      aGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAc                                                                               <  2:1634940/100‑1 (MQ=255)
                       ggCGCATAAATGTTTCCTTACTGGTTTTTTTTTTTTTTTTGATGACAGCAAGGAAAAGCTGAGTAAAGGTCGGAAAATTCGGACCCGATGGAATTTTACt                                                                              >  1:2460951/1‑100 (MQ=255)
                       ggCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACt                                                                              >  1:939192/1‑100 (MQ=255)
                       ggCGCATAAATGTTTCCTTACTGGTTTTTGTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTAc                                                                               <  1:2838059/99‑1 (MQ=255)
                        gcgcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTg                                                                             <  1:2426719/100‑1 (MQ=255)
                        gcgcATAAATGTTTCCTGACTGGTTTTTTTCTAAGTGTAGAAGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGATCCGATGGAATTTTACTg                                                                             <  1:1819275/100‑1 (MQ=255)
                         cgcATAAATGTTTCCTTACTGGTTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGGCCCGATGGAATTTTACTg                                                                             <  1:3210344/100‑1 (MQ=255)
                         cgcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGt                                                                            <  2:2334073/100‑1 (MQ=255)
                          gcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATCTCGGCGATCGCCGGCTCATTAAAAGTCAGAAAATTCGGACCCGATGGAATTTTACTGTa                                                                           >  1:1561038/1‑100 (MQ=255)
                           cATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTa                                                                           >  1:1757732/1‑99 (MQ=255)
                             tAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTagag                                                                        <  2:2019728/100‑1 (MQ=255)
                              aaaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc                                                                       <  2:4412737/100‑1 (MQ=255)
                              aaaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc                                                                       <  1:2011156/100‑1 (MQ=255)
                                  ttttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTa                                                                   <  1:879531/99‑1 (MQ=255)
                                       cTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAAt                                                              >  2:488029/1‑100 (MQ=255)
                                       cTTACTGGTTTTTTTCTAAGTGTAGATCACAGCTAGGAAAGAAGGAGTAAAGGTCAGAAAATTCGGCCTCGATGGGATCTTACAGTATAGCGTTaaaga                                                               >  1:3992965/1‑96 (MQ=255)
                                        ttACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATc                                                             <  1:2582678/100‑1 (MQ=255)
                                         tACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCt                                                            >  2:1024955/1‑100 (MQ=255)
                                           cTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTg                                                           >  1:1658103/1‑99 (MQ=255)
                                           cTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                          <  2:822862/100‑1 (MQ=255)
                                            ttccgatctttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGc                                                          <  2:2453132/92‑1 (MQ=255)
                                              gTTTTTTGCTAAGGGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAg                                                       <  1:3935893/100‑1 (MQ=255)
                                                tttttttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGATCCGAAAATTCGGACACGATGGAATTTTACTGTAGAGCGTTAAGAATATGc                                                          >  1:2521982/1‑96 (MQ=255)
                                                    tttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcggc                                                  >  2:3543585/1‑100 (MQ=255)
                                                      tCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTACCACGCGGCGt                                                >  2:1583924/1‑100 (MQ=255)
                                                       cTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                               <  2:3443311/100‑1 (MQ=255)
                                                       cTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTc                                               <  2:967140/100‑1 (MQ=255)
                                                        tAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCg                                              <  1:3795491/100‑1 (MQ=255)
                                                            tgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGActc                                          <  2:943260/100‑1 (MQ=255)
                                                             gtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGACGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGActct                                         <  2:2115643/100‑1 (MQ=255)
                                                              tAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACtctc                                        >  2:1693597/1‑100 (MQ=255)
                                                              tAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACtctc                                        <  2:468184/100‑1 (MQ=255)
                                                               aGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCg                                       <  2:193612/100‑1 (MQ=255)
                                                                 aTGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGACTCTGCCAGGCCGCGTCGACTCTCGc                                      >  2:4338440/1‑99 (MQ=255)
                                                                  tGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTATCGCTg                                    >  1:338340/1‑100 (MQ=255)
                                                                    aCAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGa                                  >  1:1221529/1‑100 (MQ=255)
                                                                     cAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAt                                 >  1:420097/1‑100 (MQ=255)
                                                                       gCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACCGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAATg                               >  1:3726339/1‑100 (MQ=255)
                                                                        cAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTgg                              >  1:3417036/1‑100 (MQ=255)
                                                                        cAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTgg                              >  1:335716/1‑100 (MQ=255)
                                                                         aaGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTggg                             <  1:3595739/100‑1 (MQ=255)
                                                                          aGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGa                            <  1:407701/100‑1 (MQ=255)
                                                                          aGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTGACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGa                            <  2:4104541/100‑1 (MQ=255)
                                                                           ggAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGAt                           >  2:1376925/1‑100 (MQ=255)
                                                                           ggAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGAt                           >  2:4361341/1‑100 (MQ=255)
                                                                                aGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttt                      >  2:4232535/1‑100 (MQ=255)
                                                                                   ggAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttg                   <  1:1142193/100‑1 (MQ=255)
                                                                                   ggAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttg                   <  1:2860384/100‑1 (MQ=255)
                                                                                    tccatcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGACAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgt                  >  2:3241386/8‑100 (MQ=255)
                                                                                    gagaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgt                  <  2:1757739/100‑1 (MQ=255)
                                                                                     agaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtt                 >  2:2554706/1‑100 (MQ=255)
                                                                                     agaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgt                  <  2:2717032/99‑1 (MQ=255)
                                                                                       atcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTTTCGCTGGATTGGGATAGTTtttgttt                >  2:4139196/5‑99 (MQ=255)
                                                                                       aaaaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg               <  2:878140/100‑1 (MQ=255)
                                                                                       aaaaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgttt                <  1:4318767/99‑1 (MQ=255)
                                                                                        aaaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg               <  1:3750131/99‑1 (MQ=255)
                                                                                        aaaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGc              >  2:2497447/1‑100 (MQ=255)
                                                                                         aaGGTCCGAAAATTCGGGCCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGcc             <  2:3598091/100‑1 (MQ=255)
                                                                                         aaGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGcc             <  1:3459446/100‑1 (MQ=255)
                                                                                          gggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCa            >  1:4282307/2‑100 (MQ=255)
                                                                                          gggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCa            >  2:4172198/2‑100 (MQ=255)
                                                                                           ggTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAg           <  2:1235498/100‑1 (MQ=255)
                                                                                             tCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCa         >  2:906651/1‑100 (MQ=255)
                                                                                               cGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATc       >  1:647916/1‑100 (MQ=255)
                                                                                               cGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATc       >  1:3443218/1‑100 (MQ=255)
                                                                                                gAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCa      >  1:522467/1‑100 (MQ=255)
                                                                                                gAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCa      >  1:2239431/1‑100 (MQ=255)
                                                                                                 aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCAc     <  2:2947674/100‑1 (MQ=255)
                                                                                                  aaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACg    <  1:2990274/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
CGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGAT  >  NC_000913/4377693‑4377891

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