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
CTTCAGTGCCTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGG  >  NC_012967/3483509‑3483664
                                                                                |                                                                            
acgtAGTGCCTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCG                                                                              >  SRR333729.4128458/5‑80 (MQ=60)
     acgtCTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAA                                                                         >  SRR333729.1373036/5‑80 (MQ=60)
     acgtCTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAAACGTCGGGGAA                                                                         >  SRR333729.849869/5‑80 (MQ=60)
       acgTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAAT                                                                       >  SRR333729.5412036/4‑80 (MQ=60)
        acgtGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATC                                                                      >  SRR333729.3745319/5‑80 (MQ=60)
         cggcaatggtggagccatggctgGCTGGCCCGTGGTGATTTAGATAACGCCCGTGGCCAGCTGAACCGTCGGGAAAAcgt                                                                     <  SRR333729.885564/57‑4 (MQ=60)
         acgtCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATac                                                                     >  SRR333729.1578034/5‑78 (MQ=60)
         acgtCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCT                                                                     >  SRR333729.4003527/5‑80 (MQ=60)
         CTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAAcgt                                                                     <  SRR333729.5861072/80‑4 (MQ=60)
         CTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCCGAACCGTCGGGAAAAcgt                                                                     <  SRR333729.4210662/80‑4 (MQ=60)
          acgtAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTG                                                                    >  SRR333729.6949908/5‑80 (MQ=60)
          TGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAAacgt                                                                    <  SRR333729.7221544/80‑5 (MQ=60)
           acgtATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGC                                                                   >  SRR333729.2056615/5‑80 (MQ=60)
           GGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATacgt                                                                   <  SRR333729.1269094/80‑5 (MQ=60)
           GGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATacgt                                                                   <  SRR333729.6927725/80‑5 (MQ=60)
             acgTGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTG                                                                 >  SRR333729.3782843/4‑80 (MQ=60)
              acgtTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTNAACCGTCGGGAAAATCTGCTGG                                                                >  SRR333729.2796649/5‑80 (MQ=60)
               ACGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGGAAATTTGCTGGG                                                               >  SRR333729.6489003/1‑80 (MQ=60)
               ACGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGG                                                               >  SRR333729.5579686/1‑80 (MQ=60)
                TGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTacgt                                                              <  SRR333729.5844784/80‑5 (MQ=60)
                 GTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGacgt                                                             <  SRR333729.3838022/80‑5 (MQ=60)
                   TGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGacgt                                                           <  SRR333729.3277868/80‑5 (MQ=60)
                   TGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGacgt                                                           <  SRR333729.5082268/80‑5 (MQ=60)
                     acgtAATGCTCGCTGGCCCGTGGTGATCTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGGAAATCTGCTGGGGAATGG                                                         >  SRR333729.4454858/5‑80 (MQ=60)
                      GTCAGTGCTCGCTGTCCCGTGGTGATTTAGATAACGCCCTTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAAacgt                                                        <  SRR333729.2232590/80‑5 (MQ=60)
                       tacagGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATacgt                                                       <  SRR333729.2851846/75‑5 (MQ=60)
                       acgtTGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAAAGGCA                                                       >  SRR333729.3943819/5‑80 (MQ=60)
                       TCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATacgt                                                       <  SRR333729.5792507/80‑5 (MQ=60)
                       TCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCACGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATacgt                                                       <  SRR333729.1054127/80‑5 (MQ=60)
                        CAATGCTCTCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGacgt                                                      <  SRR333729.6079751/80‑5 (MQ=60)
                         AATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGacgt                                                     <  SRR333729.7010233/80‑5 (MQ=60)
                          acgtTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAAT                                                    >  SRR333729.3380926/5‑80 (MQ=60)
                          ATGCTCGCTGGCCCGTGGTGATTTAGATAACGCACGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAcgt                                                    <  SRR333729.5723028/80‑4 (MQ=60)
                           TGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAcgt                                                   <  SRR333729.5180236/80‑4 (MQ=60)
                            GCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAAcgt                                                  <  SRR333729.3091276/80‑4 (MQ=60)
                             acgtCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATC                                                 >  SRR333729.4071413/5‑80 (MQ=60)
                             GTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAAacgt                                                 <  SRR333729.1088169/80‑5 (MQ=60)
                             CTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAAacgt                                                 <  SRR333729.2447375/80‑5 (MQ=60)
                              TCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATAcgt                                                <  SRR333729.554841/80‑4 (MQ=60)
                               acgTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCACATATCAC                                               >  SRR333729.1426554/4‑80 (MQ=60)
                                GCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATacgt                                              <  SRR333729.4869334/80‑5 (MQ=60)
                                 acgtCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAG                                             >  SRR333729.6227202/5‑80 (MQ=60)
                                  acgtCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGC                                            >  SRR333729.2963722/5‑80 (MQ=60)
                                   acgtCGTGGTGATTTAGATAACGCCCGTAGCCCGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCG                                           >  SRR333729.3873151/5‑80 (MQ=60)
                                   GGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAcgt                                           <  SRR333729.5781971/80‑4 (MQ=60)
                                   GGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAcgt                                           <  SRR333729.2146743/80‑4 (MQ=60)
                                    acgtGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGA                                          >  SRR333729.744326/5‑80 (MQ=60)
                                    acgtGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGA                                          >  SRR333729.2156296/5‑80 (MQ=60)
                                    GCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAacgt                                          <  SRR333729.7245871/80‑5 (MQ=60)
                                    GCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAacgt                                          <  SRR333729.865729/80‑5 (MQ=60)
                                     CCCGTGGTGATTTAGATAACGCCCGTAGCTAGCGGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGacgt                                         <  SRR333729.4289279/80‑5 (MQ=60)
                                      CCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCacgt                                        <  SRR333729.1881947/80‑5 (MQ=60)
                                      ACGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACG                                        >  SRR333729.5369062/1‑80 (MQ=60)
                                       acgtGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCccaacgacag                                       >  SRR333729.2660964/5‑70 (MQ=60)
                                       CGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACgt                                       <  SRR333729.1546218/80‑3 (MQ=60)
                                       CGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACgt                                       <  SRR333729.6545313/80‑3 (MQ=60)
                                        acgtTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGG                                      >  SRR333729.1961686/5‑80 (MQ=60)
                                        GTGGTGATTTAGTTAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGAacgt                                      <  SRR333729.5696365/80‑5 (MQ=60)
                                          GGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTacgt                                    <  SRR333729.1723987/80‑5 (MQ=60)
                                           GTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGNGGAATGGCAAATATCACAGCGACTGacgt                                   <  SRR333729.806776/80‑5 (MQ=60)
                                           GTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGacgt                                   <  SRR333729.248905/80‑5 (MQ=60)
                                            acgTTTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACGGGAGCT                                  >  SRR333729.109162/4‑80 (MQ=60)
                                              acgTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCT                                >  SRR333729.4970646/4‑80 (MQ=60)
                                               acgtGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTA                               >  SRR333729.5779265/5‑80 (MQ=60)
                                                TTAGATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTacgt                              <  SRR333729.5244472/80‑5 (MQ=60)
                                                 acgtTAACGCCCGTAGCCAGCTGAACCGTCGGGAANATCTGCNGGGGAATGGCAAATATCACAGCGACTGGATCTCTAAC                             >  SRR333729.419223/5‑80 (MQ=60)
                                                 acgtTAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAAC                             >  SRR333729.5265480/5‑80 (MQ=60)
                                                 TAGATAACGCCCGTAGCCAGCTGAACCGTCGGGANAATCTGCTGGGGANTGGCAAATATCACAGCGACTGGATCTCacgt                             <  SRR333729.362441/80‑5 (MQ=60)
                                                  acgTAACGCCCGTAGCCCGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGCCTGGATtctaaagg                            >  SRR333729.3126489/4‑72 (MQ=60)
                                                  acgTAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATNGCAAATATCACAGCGACTGGATCTCTAACG                            >  SRR333729.498348/4‑80 (MQ=60)
                                                  acgTAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACG                            >  SRR333729.1565299/4‑80 (MQ=60)
                                                  acgTAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACG                            >  SRR333729.6742956/4‑80 (MQ=60)
                                                   acgtACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGC                           >  SRR333729.4674736/5‑80 (MQ=60)
                                                   GATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGT                           <  SRR333729.5324297/80‑1 (MQ=60)
                                                    ATCACGCCCGTATCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAAacgt                          <  SRR333729.2677456/80‑5 (MQ=60)
                                                    ATAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAAacgt                          <  SRR333729.1254138/80‑5 (MQ=60)
                                                     acgtGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATTGCAAATATCACAGCGACTGGATCTCTAACGCCA                         >  SRR333729.3494652/5‑80 (MQ=60)
                                                     acgtGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA                         >  SRR333729.566275/5‑80 (MQ=60)
                                                     TAACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACacgt                         <  SRR333729.6939876/80‑5 (MQ=60)
                                                      aacccgccgtagcggGCTGACCCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACacgt                         <  SRR333729.2639477/65‑5 (MQ=60)
                                                      AACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGacgt                        <  SRR333729.1675840/80‑5 (MQ=60)
                                                      AACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGacgt                        <  SRR333729.6393810/80‑5 (MQ=60)
                                                      AACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGacgt                        <  SRR333729.982526/80‑5 (MQ=60)
                                                       ACGTCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAA                        >  SRR333729.6157259/1‑80 (MQ=60)
                                                       ACGTCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAA                        >  SRR333729.6935606/1‑80 (MQ=60)
                                                       ACGTCCGCAGCCACCTGAACCGTCGGGAAAATCTGCCGGGGGATGGCAAATATCACAACGACTGGATCTCTAACGCCAAC                       >  SRR333729.4124875/1‑80 (MQ=56)
                                                       ACGCCCGTAGCCAGCTGAACCGTCGGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCacgt                       <  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)
                                                                                |                                                                            
CTTCAGTGCCTGGCAATGTTGAGTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAATCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCAACAAAGTCCGGGTGATTTACTGG  >  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: