Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 3 195.4 13365061 92.0% 12295856 72.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_012967 3,483,588 T→G L635R (CTG→CGG)  malT → HTH‑type transcriptional regulator MalT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0129673,483,5880TG99.0% 818.7 / ‑3.7 191L635R (CTG→CGG) malTHTH‑type transcriptional regulator MalT
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (76/113);  total (78/113)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01

Only 100 of 202 total aligned reads displayed.

TGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGA  >  NC_012967/3483519‑3483656
                                                                     |                                                                    
 ggCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGa                                                                   <  1:5071642/72‑1 (MQ=255)
  gCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGaa                                                                  >  1:3436614/1‑72 (MQ=255)
    aaTGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGaaaa                                                                <  2:3559794/72‑1 (MQ=255)
    aaTGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGaaaa                                                                >  2:3602593/1‑72 (MQ=255)
    aaTGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGaaaa                                                                >  2:2828789/1‑72 (MQ=255)
    aaTGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGaaaa                                                                <  1:6642516/72‑1 (MQ=255)
     aTGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAAt                                                               >  2:2983181/1‑72 (MQ=255)
     aTGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAAt                                                               >  1:1882945/1‑72 (MQ=255)
     aTGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAAt                                                               <  1:6369791/72‑1 (MQ=255)
      tGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATc                                                              <  2:5625322/72‑1 (MQ=255)
        ttGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTNAACCGTCGGGAAAATctg                                                            >  1:2562351/1‑72 (MQ=255)
         tGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATctgc                                                           <  2:5349547/72‑1 (MQ=255)
          taGTCAATGCTGGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATctgct                                                          <  2:3841895/71‑1 (MQ=255)
          gAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATctgct                                                          <  1:5368621/72‑1 (MQ=255)
           aGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATctgctg                                                         <  1:3522426/72‑1 (MQ=255)
              cAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTgggg                                                      <  2:2056258/72‑1 (MQ=255)
               aaTGCTCGCTGGCCCTTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGa                                                     <  2:5765892/72‑1 (MQ=255)
               aaTGCTCGCTGGCCCGTGGTGATCTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGGAAATCTGCTGGGGa                                                     >  1:4088509/1‑72 (MQ=255)
                 tGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGaaa                                                   >  1:3620288/1‑71 (MQ=255)
                 tGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAAt                                                   <  2:1423206/72‑1 (MQ=255)
                 tGCTCGCTGGCCCGTGGTGATTTAGATAACGCACGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAAt                                                   <  1:959016/72‑1 (MQ=255)
                 ggCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAAt                                                   <  1:2612984/71‑1 (MQ=255)
                   cTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATgg                                                 <  1:6446214/72‑1 (MQ=255)
                    tCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGc                                                >  1:3102210/1‑72 (MQ=255)
                    tCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGc                                                <  2:6539388/72‑1 (MQ=255)
                    tCGCTGGCCCGTGGTGATTTAGATAACGCACGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGc                                                <  1:5255731/72‑1 (MQ=255)
                      gCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCaa                                              <  1:2834306/72‑1 (MQ=255)
                      gCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCaa                                              <  2:5956734/72‑1 (MQ=255)
                       cTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCaaa                                             <  2:5049582/72‑1 (MQ=255)
                       cTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCaaa                                             <  1:989837/72‑1 (MQ=255)
                       cTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCaaa                                             >  1:3738041/1‑72 (MQ=255)
                         ggCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAata                                           >  2:6380857/1‑72 (MQ=255)
                          gCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAatat                                          <  1:4467673/72‑1 (MQ=255)
                          gCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGCGAATGGCAAatat                                          >  2:5150731/1‑72 (MQ=255)
                             cGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATcac                                       <  1:5310218/72‑1 (MQ=255)
                              gtggtgATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATcaca                                      >  1:679600/1‑72 (MQ=255)
                              gtggtgATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATcaca                                      >  1:1975598/1‑72 (MQ=255)
                               tggtgATTTAGATAACGCCCGTAGCTAGCGGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAg                                     <  1:3937907/72‑1 (MQ=255)
                               tggtgATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAg                                     <  2:1880473/72‑1 (MQ=255)
                               tggtgATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAg                                     >  2:5976622/1‑72 (MQ=255)
                                ggtgATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGc                                    >  1:4927707/1‑72 (MQ=255)
                                ggtgATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGc                                    <  1:1720861/72‑1 (MQ=255)
                                ggtgATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGc                                    <  2:1971580/72‑1 (MQ=255)
                                 gtgATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCCCAACg                                   >  1:2438103/1‑72 (MQ=255)
                                 gtgATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCg                                   <  1:6014602/72‑1 (MQ=255)
                                   gATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGAc                                 <  2:1213986/72‑1 (MQ=255)
                                     tttAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTg                               <  2:1320915/72‑1 (MQ=255)
                                     tttAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTg                               <  1:226429/72‑1 (MQ=255)
                                      tgaGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTgg                              <  2:3150104/70‑1 (MQ=255)
                                       tAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGa                             >  2:5031434/1‑72 (MQ=255)
                                        aGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGAt                            >  2:6033339/1‑72 (MQ=255)
                                        aGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGAt                            >  1:4561272/1‑72 (MQ=255)
                                         gATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGAtc                           <  2:3310229/72‑1 (MQ=255)
                                         caTAACGCCCGTAGCAAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGAtc                           <  2:2439730/71‑1 (MQ=255)
                                          aTAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGAtct                          <  1:4813058/72‑1 (MQ=255)
                                           tAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGAtctc                         >  2:6443697/1‑72 (MQ=255)
                                           tAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGAtctc                         >  1:4832419/1‑72 (MQ=255)
                                            aaCGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATctct                        <  2:5271509/72‑1 (MQ=255)
                                            aaCGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATctct                        >  1:1430329/1‑72 (MQ=255)
                                             actcccGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATTGCAAATATCACAGCGACTGGATCTCTa                       <  2:274960/69‑1 (MQ=255)
                                              cGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTaa                      <  1:1142318/72‑1 (MQ=255)
                                              cGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTaa                      <  2:4960655/72‑1 (MQ=255)
                                               gCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAAc                     <  1:6381085/72‑1 (MQ=255)
                                               gCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAAc                     >  1:518101/1‑72 (MQ=255)
                                                cccGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACg                    <  1:893932/72‑1 (MQ=255)
                                                cccGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACg                    >  2:4476102/1‑72 (MQ=255)
                                                cccGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACg                    >  2:2283230/1‑72 (MQ=255)
                                                cccGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACg                    <  1:5874533/72‑1 (MQ=255)
                                                cccGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACg                    >  1:5657304/1‑72 (MQ=255)
                                                cccGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACg                    <  1:1531835/72‑1 (MQ=255)
                                                cccGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACg                    <  2:3530497/72‑1 (MQ=255)
                                                 ccGTANCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGc                   <  2:591404/72‑1 (MQ=255)
                                                 ccGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGc                   <  1:1209221/72‑1 (MQ=255)
                                                 ccGTACCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGc                   <  2:798574/72‑1 (MQ=255)
                                                   gTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCca                 >  1:5122500/1‑72 (MQ=255)
                                                    tAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCcaa                >  1:79234/1‑72 (MQ=255)
                                                    tAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATACATAACGCcaa                >  1:4686139/1‑72 (MQ=255)
                                                     aGCCAGCTGAACCGTCTGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCcaac               >  1:6142180/1‑72 (MQ=255)
                                                     aGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCcaac               >  1:2416676/1‑72 (MQ=255)
                                                     aGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCcaac               >  1:3206727/1‑72 (MQ=255)
                                                      gCCAGCTGAACCGTCGGGAAAATCTGCTGGTGAATGGCAAATATCACAGCGACTGGATCTCTAACGCcaaca              >  1:2071992/1‑72 (MQ=255)
                                                      gCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCcaaca              <  2:4908870/72‑1 (MQ=255)
                                                      gCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCcaaca              <  2:1261266/72‑1 (MQ=255)
                                                      gCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCcaaca              >  1:5912469/1‑72 (MQ=255)
                                                        cAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCNACaaa            >  2:2527223/1‑72 (MQ=255)
                                                        cAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACaaa            >  1:137057/1‑72 (MQ=255)
                                                         aGCTGAACCGTCGGGAAANTCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAg           <  1:2205676/72‑1 (MQ=255)
                                                         aGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAg           <  2:3889982/72‑1 (MQ=255)
                                                          gCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGt          <  2:2716121/72‑1 (MQ=255)
                                                          gCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGt          <  1:1078108/72‑1 (MQ=255)
                                                          gCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGt          <  2:4760285/72‑1 (MQ=255)
                                                          gCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGt          <  2:2204627/72‑1 (MQ=255)
                                                           cTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTc         <  2:3456291/72‑1 (MQ=255)
                                                            tGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTcc        >  1:5450883/1‑72 (MQ=255)
                                                             gAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCg       >  2:1708055/1‑72 (MQ=255)
                                                               aCCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCgtg     >  1:1964169/1‑70 (MQ=255)
                                                               aCCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCggg     <  2:83264/72‑1 (MQ=255)
                                                                ccGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGt    >  2:1454026/1‑72 (MQ=255)
                                                                ccGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGt    <  1:3837571/72‑1 (MQ=255)
                                                                  gTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGa  <  2:1530256/72‑1 (MQ=255)
                                                                     |                                                                    
TGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGA  >  NC_012967/3483519‑3483656

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTTCAGTGCCTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCTGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGG  >  NC_012967/3483509‑3483664
                                                                                |                                                                            
acgtAGTGCCTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCG                                                                              >  SRR333729.4128458/5‑80 (MQ=60)
     acgtCTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAA                                                                         >  SRR333729.1373036/5‑80 (MQ=60)
     acgtCTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAAACGTCGGGGAA                                                                         >  SRR333729.849869/5‑80 (MQ=60)
       acgTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAAT                                                                       >  SRR333729.5412036/4‑80 (MQ=60)
        acgtGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATC                                                                      >  SRR333729.3745319/5‑80 (MQ=60)
         cggcaatggtggagccatggctgGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTGGCCAGCTGAACCGTCGGGAAAAcgt                                                                     <  SRR333729.885564/57‑4 (MQ=60)
         acgtCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATac                                                                     >  SRR333729.1578034/5‑78 (MQ=60)
         acgtCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCT                                                                     >  SRR333729.4003527/5‑80 (MQ=60)
         CTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAAcgt                                                                     <  SRR333729.5861072/80‑4 (MQ=60)
         CTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCCGAACCGTCGGGAAAAcgt                                                                     <  SRR333729.4210662/80‑4 (MQ=60)
          acgtAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTG                                                                    >  SRR333729.6949908/5‑80 (MQ=60)
          TGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAAacgt                                                                    <  SRR333729.7221544/80‑5 (MQ=60)
           acgtATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGC                                                                   >  SRR333729.2056615/5‑80 (MQ=60)
           GGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATacgt                                                                   <  SRR333729.1269094/80‑5 (MQ=60)
           GGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATacgt                                                                   <  SRR333729.6927725/80‑5 (MQ=60)
             acgTGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTG                                                                 >  SRR333729.3782843/4‑80 (MQ=60)
              acgtTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTNAACCGTCGGGAAAATCTGCTGG                                                                >  SRR333729.2796649/5‑80 (MQ=60)
               ACGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGGAAATTTGCTGGG                                                               >  SRR333729.6489003/1‑80 (MQ=60)
               ACGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGG                                                               >  SRR333729.5579686/1‑80 (MQ=60)
                TGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTacgt                                                              <  SRR333729.5844784/80‑5 (MQ=60)
                 GTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGacgt                                                             <  SRR333729.3838022/80‑5 (MQ=60)
                   TGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGacgt                                                           <  SRR333729.3277868/80‑5 (MQ=60)
                   TGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGacgt                                                           <  SRR333729.5082268/80‑5 (MQ=60)
                     acgtAATGCTCGCTGGCCCGTGGTGATCTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGGAAATCTGCTGGGGAATGG                                                         >  SRR333729.4454858/5‑80 (MQ=60)
                      GTCAGTGCTCGCTGTCCCGTGGTGATTTAGATAACG‑CCCTTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAAacgt                                                        <  SRR333729.2232590/80‑5 (MQ=60)
                       tacagGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATacgt                                                       <  SRR333729.2851846/75‑5 (MQ=60)
                       acgtTGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAAAGGCA                                                       >  SRR333729.3943819/5‑80 (MQ=60)
                       TCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATacgt                                                       <  SRR333729.5792507/80‑5 (MQ=60)
                       TCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CACGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATacgt                                                       <  SRR333729.1054127/80‑5 (MQ=60)
                        CAATGCTCTCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGacgt                                                      <  SRR333729.6079751/80‑5 (MQ=60)
                         AATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGacgt                                                     <  SRR333729.7010233/80‑5 (MQ=60)
                          acgtTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAAT                                                    >  SRR333729.3380926/5‑80 (MQ=60)
                          ATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CACGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAcgt                                                    <  SRR333729.5723028/80‑4 (MQ=60)
                           TGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAcgt                                                   <  SRR333729.5180236/80‑4 (MQ=60)
                            GCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAAcgt                                                  <  SRR333729.3091276/80‑4 (MQ=60)
                             acgtCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATC                                                 >  SRR333729.4071413/5‑80 (MQ=60)
                             GTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAAacgt                                                 <  SRR333729.1088169/80‑5 (MQ=60)
                             CTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAAacgt                                                 <  SRR333729.2447375/80‑5 (MQ=60)
                              TCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATAcgt                                                <  SRR333729.554841/80‑4 (MQ=60)
                               acgTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCACATATCAC                                               >  SRR333729.1426554/4‑80 (MQ=60)
                                GCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATacgt                                              <  SRR333729.4869334/80‑5 (MQ=60)
                                 acgtCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAG                                             >  SRR333729.6227202/5‑80 (MQ=60)
                                  acgtCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGC                                            >  SRR333729.2963722/5‑80 (MQ=60)
                                   acgtCGTGGTGATTTAGATAACG‑CCCGTAGCCCGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCG                                           >  SRR333729.3873151/5‑80 (MQ=60)
                                   GGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAcgt                                           <  SRR333729.5781971/80‑4 (MQ=60)
                                   GGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAcgt                                           <  SRR333729.2146743/80‑4 (MQ=60)
                                    acgtGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGA                                          >  SRR333729.744326/5‑80 (MQ=60)
                                    acgtGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGA                                          >  SRR333729.2156296/5‑80 (MQ=60)
                                    GCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAacgt                                          <  SRR333729.7245871/80‑5 (MQ=60)
                                    GCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAacgt                                          <  SRR333729.865729/80‑5 (MQ=60)
                                     CCCGTGGTGATTTAGATAACG‑CCCGTAGCTAGCGGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGacgt                                         <  SRR333729.4289279/80‑5 (MQ=60)
                                      CCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCacgt                                        <  SRR333729.1881947/80‑5 (MQ=60)
                                      ACGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACG                                        >  SRR333729.5369062/1‑80 (MQ=60)
                                       acgtGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCccaacgacag                                       >  SRR333729.2660964/5‑70 (MQ=60)
                                       CGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACgt                                       <  SRR333729.1546218/80‑3 (MQ=60)
                                       CGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACgt                                       <  SRR333729.6545313/80‑3 (MQ=60)
                                        acgtTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGG                                      >  SRR333729.1961686/5‑80 (MQ=60)
                                        GTGGTGATTTAGTTAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGAacgt                                      <  SRR333729.5696365/80‑5 (MQ=60)
                                          GGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTacgt                                    <  SRR333729.1723987/80‑5 (MQ=60)
                                           GTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGNGGAATGGCAAATATCACAGCGACTGacgt                                   <  SRR333729.806776/80‑5 (MQ=60)
                                           GTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGacgt                                   <  SRR333729.248905/80‑5 (MQ=60)
                                            acgTTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACGGGAGCT                                  >  SRR333729.109162/4‑80 (MQ=60)
                                              acgTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCT                                >  SRR333729.4970646/4‑80 (MQ=60)
                                               acgtGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTA                               >  SRR333729.5779265/5‑80 (MQ=60)
                                                TTAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTacgt                              <  SRR333729.5244472/80‑5 (MQ=60)
                                                 acgtTAACG‑CCCGTAGCCAGCTGAACCGTCGGGAANATCTGCNGGGGAATGGCAAATATCACAGCGACTGGATCTCTAAC                             >  SRR333729.419223/5‑80 (MQ=60)
                                                 acgtTAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAAC                             >  SRR333729.5265480/5‑80 (MQ=60)
                                                 TAGATAACG‑CCCGTAGCCAGCTGAACCGTCGGGANAATCTGCTGGGGANTGGCAAATATCACAGCGACTGGATCTCacgt                             <  SRR333729.362441/80‑5 (MQ=60)
                                                  acgTAACG‑CCCGTAGCCCGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGCCTGGATtctaaagg                            >  SRR333729.3126489/4‑72 (MQ=60)
                                                  acgTAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATNGCAAATATCACAGCGACTGGATCTCTAACG                            >  SRR333729.498348/4‑80 (MQ=60)
                                                  acgTAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACG                            >  SRR333729.1565299/4‑80 (MQ=60)
                                                  acgTAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACG                            >  SRR333729.6742956/4‑80 (MQ=60)
                                                   acgtACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGC                           >  SRR333729.4674736/5‑80 (MQ=60)
                                                   GATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGT                           <  SRR333729.5324297/80‑1 (MQ=60)
                                                    ATCACG‑CCCGTATCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAAacgt                          <  SRR333729.2677456/80‑5 (MQ=60)
                                                    ATAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAAacgt                          <  SRR333729.1254138/80‑5 (MQ=60)
                                                     acgtG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATTGCAAATATCACAGCGACTGGATCTCTAACGCCA                         >  SRR333729.3494652/5‑80 (MQ=60)
                                                     acgtG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA                         >  SRR333729.566275/5‑80 (MQ=60)
                                                     TAACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACacgt                         <  SRR333729.6939876/80‑5 (MQ=60)
                                                      aacccgccgtagcggGCTGACCCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACacgt                         <  SRR333729.2639477/65‑5 (MQ=60)
                                                      AACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGacgt                        <  SRR333729.1675840/80‑5 (MQ=60)
                                                      AACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGacgt                        <  SRR333729.6393810/80‑5 (MQ=60)
                                                      AACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGacgt                        <  SRR333729.982526/80‑5 (MQ=60)
                                                       ACGTCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAA                        >  SRR333729.6157259/1‑80 (MQ=60)
                                                       ACGTCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAA                        >  SRR333729.6935606/1‑80 (MQ=60)
                                                       ACG‑TCCGCAGCCACCTGAACCGTCGGGAAAATCTGCCGGGGGATGGCAAATATCACAACGACTGGATCTCTAACGCCAAC                       >  SRR333729.4124875/1‑80 (MQ=56)
                                                       ACG‑CCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCacgt                       <  SRR333729.1326471/80‑5 (MQ=60)
                                                          acgtGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAA                     >  SRR333729.2131400/5‑80 (MQ=60)
                                                          acgtGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAA                     >  SRR333729.5579432/5‑80 (MQ=60)
                                                           acgtTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAA                    >  SRR333729.87963/5‑80 (MQ=60)
                                                           acgtTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATACATAACGCCAACAAA                    >  SRR333729.5105727/5‑80 (MQ=60)
                                                            ACGTAGCCAGCTGAACCGTCTGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAG                   >  SRR333729.6682927/1‑80 (MQ=60)
                                                            ACGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAG                   >  SRR333729.2637263/1‑80 (MQ=60)
                                                            ACGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAG                   >  SRR333729.3494885/1‑80 (MQ=60)
                                                             acgtGCCAGCTGAACCGTCGGGAAAATCTGCTGGTGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGT                  >  SRR333729.2259739/5‑80 (MQ=60)
                                                             acgtGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGGAAATATCACCGCGACTGGATCTCTAACGCCAACAAAGT                  >  SRR333729.2501329/5‑80 (MQ=60)
                                                             acgtGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGT                  >  SRR333729.6434898/5‑80 (MQ=60)
                                                             CGTAGCCAGCTGAACCGTCGGGNAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAACGT                  <  SRR333729.2799121/80‑1 (MQ=60)
                                                               acgtCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCC                >  SRR333729.4416734/5‑80 (MQ=60)
                                                               acgtCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCC                >  SRR333729.151516/5‑80 (MQ=60)
                                                                agcaggaTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCAAACAAAGacgt               <  SRR333729.3024169/73‑5 (MQ=60)
                                                                acgtAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCG               >  SRR333729.6598742/5‑80 (MQ=60)
                                                                AGCCAGCTGAACCGTCGGGAAANTCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGacgt               <  SRR333729.2404635/80‑5 (MQ=60)
                                                                 GCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTacgt              <  SRR333729.6415571/80‑5 (MQ=60)
                                                                 GCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTacgt              <  SRR333729.1184487/80‑5 (MQ=60)
                                                                   acgtTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGT            >  SRR333729.5933433/5‑80 (MQ=60)
                                                                   CAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCacgt            <  SRR333729.1249690/80‑5 (MQ=60)
                                                                   CAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCacgt            <  SRR333729.942961/80‑5 (MQ=60)
                                                                    AGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGacgt           <  SRR333729.5246412/80‑5 (MQ=60)
                                                                      acgtACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGTGTGAT         >  SRR333729.2144044/5‑80 (MQ=60)
                                                                       TGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTacgt        <  SRR333729.7129628/80‑5 (MQ=60)
                                                                       TGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTacgt        <  SRR333729.4179677/80‑5 (MQ=60)
                                                                             GTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTAacgt  <  SRR333729.3251448/80‑5 (MQ=60)
                                                                                |                                                                            
CTTCAGTGCCTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACG‑CCCGTAGCCAGCTGAACCGTCTGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGG  >  NC_012967/3483509‑3483664

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 2 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑