Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,158,718 T→C Q686R (CAG→CGG)  fhuE ← outer membrane receptor for ferric iron uptake

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,158,7180TC99.3% 924.1 / ‑5.1 272Q686R (CAG→CGG) fhuEouter membrane receptor for ferric iron uptake
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (146/124);  total (147/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

CCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGA  >  NC_000913/1158628‑1158807
                                                                                          |                                                                                         
ccGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCgg                                                                                         <  1:2868640/93‑1 (MQ=255)
ccGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCgg                                                                                         >  2:3534178/1‑93 (MQ=255)
 cGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGt                                                                                        <  1:4738522/93‑1 (MQ=255)
 cGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGt                                                                                        <  1:1891469/93‑1 (MQ=255)
 cGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGt                                                                                        <  2:7081172/93‑1 (MQ=255)
 cGTAGACGATAGAACCTTCCACGTTGGTATAGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGt                                                                                        >  1:7474954/1‑93 (MQ=255)
  gTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTa                                                                                       <  2:6582987/93‑1 (MQ=255)
  gTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTa                                                                                       <  1:3495980/93‑1 (MQ=255)
  gTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCCAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTa                                                                                       <  2:7239237/93‑1 (MQ=255)
   tAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAg                                                                                      <  1:2988548/93‑1 (MQ=255)
   tAGACGATAGAACCTTCCACGTTGGGATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGGTTTTCGTCACCCGGTAg                                                                                      >  1:8100215/1‑93 (MQ=255)
     gACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCg                                                                                    >  2:5977821/1‑93 (MQ=255)
     gACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCg                                                                                    <  2:91030/93‑1 (MQ=255)
        gATAGACCCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGt                                                                                 <  2:2534877/93‑1 (MQ=255)
        gATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGCAGCGGGt                                                                                 >  2:331050/1‑93 (MQ=255)
         aTAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTg                                                                                <  2:5954267/93‑1 (MQ=255)
           aGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGaa                                                                              >  2:4356865/1‑93 (MQ=255)
             aaCCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAaga                                                                            <  1:8263498/93‑1 (MQ=255)
             aaCCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAaga                                                                            <  1:1049201/93‑1 (MQ=255)
               ccTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAgaga                                                                          <  1:6632156/93‑1 (MQ=255)
               ccTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAgaga                                                                          >  1:2535563/1‑93 (MQ=255)
               ccTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAgaga                                                                          <  2:2799135/93‑1 (MQ=255)
                cTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGAt                                                                         <  1:3256613/93‑1 (MQ=255)
                cTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGAt                                                                         >  1:5507534/1‑93 (MQ=255)
                 ttCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCCCCCGGTAGCGGGCGGAGAGATc                                                                        >  2:1739584/1‑93 (MQ=255)
                 ttCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATc                                                                        <  2:6223663/93‑1 (MQ=255)
                 ttCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATc                                                                        >  2:1416326/1‑93 (MQ=255)
                    cACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATccac                                                                     <  1:7373970/93‑1 (MQ=255)
                    cACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATccac                                                                     >  1:8231545/1‑93 (MQ=255)
                      cGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACTAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATccacca                                                                   <  2:4286576/93‑1 (MQ=255)
                       gTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAg                                                                  <  1:5062400/93‑1 (MQ=255)
                       gTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAg                                                                  <  2:163360/93‑1 (MQ=255)
                       gTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAg                                                                  >  1:998518/1‑93 (MQ=255)
                       gTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAg                                                                  >  1:8010780/1‑93 (MQ=255)
                       gTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAg                                                                  >  1:7625912/1‑93 (MQ=255)
                        ttGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAgc                                                                 >  2:8644597/1‑93 (MQ=255)
                         tGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAgcg                                                                <  1:6987696/93‑1 (MQ=255)
                          ggTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAgcgc                                                               <  1:5801654/93‑1 (MQ=255)
                           gtgtCGTAGCTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCCCCCGGTAGCGGGTGAAGAGATCCACCAGcgcg                                                              <  2:6298935/90‑1 (MQ=255)
                           gTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGcgcg                                                              >  2:4471887/1‑93 (MQ=255)
                           gTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGcgcg                                                              >  1:4626314/1‑93 (MQ=255)
                            tATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGt                                                             <  1:1001362/93‑1 (MQ=255)
                              tCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAg                                                           >  1:7856024/1‑93 (MQ=255)
                               cGTAGGTTTTGTCGAACAGGTTATTGCCGTTCCCCTGTAACGAGCAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGc                                                          <  1:2793884/93‑1 (MQ=255)
                                ttAGGGTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCt                                                         <  1:2463290/92‑1 (MQ=255)
                                gTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCt                                                         >  1:8476148/1‑93 (MQ=255)
                                gTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCt                                                         >  2:1902256/1‑93 (MQ=255)
                                gTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCt                                                         >  1:4056853/1‑93 (MQ=255)
                                 tAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTTACCCGGTAGCGGGTGAAGAGAGCCACCCGCGCGTAGCTa                                                        >  1:7164907/1‑93 (MQ=255)
                                 tAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTa                                                        >  2:8242778/1‑93 (MQ=255)
                                 tAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTa                                                        >  1:7746518/1‑93 (MQ=255)
                                 tAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTa                                                        >  2:6972245/1‑93 (MQ=255)
                                 tAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTa                                                        >  2:1840806/1‑93 (MQ=255)
                                  aGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTAc                                                       >  1:7224837/1‑93 (MQ=255)
                                  aGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTAc                                                       >  2:4709407/1‑93 (MQ=255)
                                  aGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTAc                                                       >  1:1351734/1‑93 (MQ=255)
                                    gTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCt                                                     <  2:7888618/93‑1 (MQ=255)
                                       ttGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGc                                                  >  2:1534354/1‑93 (MQ=255)
                                          tCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCg                                               <  1:1234281/93‑1 (MQ=255)
                                          tCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTAACTTGCTCg                                               >  2:879204/1‑93 (MQ=255)
                                            gAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGgc                                             <  1:2192890/93‑1 (MQ=255)
                                             aaCAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGgcg                                            >  2:1433617/1‑93 (MQ=255)
                                             aaCAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGgcg                                            >  2:3078090/1‑93 (MQ=255)
                                               cAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGcgcg                                          >  2:5730626/1‑93 (MQ=255)
                                                 ggTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGa                                        <  2:6482853/93‑1 (MQ=255)
                                                 ggTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGa                                        >  1:5776412/1‑93 (MQ=255)
                                                 ggTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGa                                        >  2:6183885/1‑93 (MQ=255)
                                                  gTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGaa                                       <  1:4903450/93‑1 (MQ=255)
                                                  gTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGaa                                       <  1:4422781/93‑1 (MQ=255)
                                                    tATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAgg                                     <  2:8715268/93‑1 (MQ=255)
                                                     aTTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGt                                    >  1:2293128/1‑93 (MQ=255)
                                                       tGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGc                                  >  1:1608003/1‑93 (MQ=255)
                                                       tGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGc                                  >  2:3785123/1‑93 (MQ=255)
                                                        gACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGcc                                 >  1:6815760/1‑93 (MQ=255)
                                                        gACGTCCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCCCCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGcc                                 <  2:8534382/93‑1 (MQ=255)
                                                            ttCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATAc                             <  1:4850827/93‑1 (MQ=255)
                                                             tCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACg                            >  2:3156774/1‑93 (MQ=255)
                                                             tCCCCTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACg                            <  1:5380203/93‑1 (MQ=255)
                                                              ccccTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACgg                           >  2:5614872/1‑93 (MQ=255)
                                                               cccTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGc                          <  1:1806865/93‑1 (MQ=255)
                                                               cccTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGc                          >  1:8794541/1‑93 (MQ=255)
                                                                 cTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGt                        >  1:3562169/1‑93 (MQ=255)
                                                                 cTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTACCTACCTGCCTCGGCGCGGGAGGTGCCATACGGCgc                        >  2:1160600/1‑92 (MQ=255)
                                                                 cTGTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGCGAAGAGATCCACCAGCGCGCAGCTACCTTGCTCGGTGCGGCAGGTGCCATACGGCga                        >  2:4946997/1‑92 (MQ=255)
                                                                   gTAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCa                      >  2:2280689/1‑93 (MQ=255)
                                                                    tAACGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCAc                     <  2:768522/93‑1 (MQ=255)
                                                                      aCGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACgg                   <  1:3921896/93‑1 (MQ=255)
                                                                      aCGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACgg                   <  1:5716011/93‑1 (MQ=255)
                                                                       cGAGAAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgt                  <  2:6762269/93‑1 (MQ=255)
                                                                        gagaAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgtg                 >  2:5355030/1‑93 (MQ=255)
                                                                        gagaAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgtg                 >  2:4546247/1‑93 (MQ=255)
                                                                         agaAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgtgt                <  1:4144591/93‑1 (MQ=255)
                                                                         agaAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgtgt                <  1:3695854/93‑1 (MQ=255)
                                                                         agaAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgtgt                <  2:2148576/93‑1 (MQ=255)
                                                                         agaAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgtgt                <  1:6011144/93‑1 (MQ=255)
                                                                         agaAGTTTTTCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGgtgt                >  1:7260080/1‑93 (MQ=255)
                                                                                 ttCGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATac        >  1:6111088/1‑93 (MQ=255)
                                                                                   cGTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATacac      >  2:7364415/1‑93 (MQ=255)
                                                                                    gTCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACAcg     >  2:1340101/1‑93 (MQ=255)
                                                                                     tCACCCGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACAcgc    <  1:2690977/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCGTAGACGATAGAACCTTCCACGTTGGTATCGTAGGTTTTGTCGAACAGGTTATTGACGTTCCCCTGTAACGAGAAGTTTTTCGTCACCTGGTAGCGGGTGAAGAGATCCACCAGCGCGTAGCTACCTTGCTCGGCGCGGAAGGTGCCATACGGCGTCACGGTGTCGGTATACACGCGA  >  NC_000913/1158628‑1158807

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