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
A2 F1 I1 R1 2 240.8 15455571 92.9% 14358225 72.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_012967 1,003,919 G→T N52K (AAC→AAA ompF ← porin OmpF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0129671,003,9190GT100.0% 564.2 / NA 164N52K (AAC→AAAompFporin OmpF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (68/96);  total (68/96)

Only 100 of 174 total aligned reads displayed.

GAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCA  >  NC_012967/1003851‑1003988
                                                                    |                                                                     
gAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTtttt                                                                    >  2:2148655/1‑72 (MQ=255)
gAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTtttt                                                                    >  2:7522381/1‑72 (MQ=255)
 aaTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACtttttt                                                                   >  1:1280467/1‑72 (MQ=255)
 aaTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACtttttt                                                                   <  2:6904819/72‑1 (MQ=255)
  aTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTc                                                                  <  1:2396828/72‑1 (MQ=255)
     gATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCAcc                                                               <  1:811599/72‑1 (MQ=255)
       tttGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGt                                                             <  2:5422402/72‑1 (MQ=255)
       tttGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGt                                                             <  1:3785911/72‑1 (MQ=255)
        ttGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGtt                                                            <  2:1282981/72‑1 (MQ=255)
        ttGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGtt                                                            <  1:5539319/72‑1 (MQ=255)
         tGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTa                                                           <  1:1514124/72‑1 (MQ=255)
          gAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTAc                                                          <  2:4694168/72‑1 (MQ=255)
          gAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTAc                                                          >  1:247723/1‑72 (MQ=255)
           aGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACAGTAACTTTTTTCACCGTTAcc                                                         <  1:5554017/72‑1 (MQ=255)
            gTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTAccc                                                        <  2:5042196/72‑1 (MQ=255)
            gTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTAccc                                                        >  1:3971716/1‑72 (MQ=255)
              ttCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCtt                                                      >  2:5047173/1‑72 (MQ=255)
              ttCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCtt                                                      >  2:7586194/1‑72 (MQ=255)
               tCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCtcg                                                     >  1:42634/1‑70 (MQ=255)
               tCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTg                                                     <  2:5596461/72‑1 (MQ=255)
               tCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTg                                                     <  1:3331416/72‑1 (MQ=255)
               tCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTg                                                     >  2:2113700/1‑72 (MQ=255)
                ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTgg                                                    >  2:2523953/1‑72 (MQ=255)
                ccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTgg                                                    >  2:5329164/1‑72 (MQ=255)
                  ccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGgag                                                  <  1:151468/72‑1 (MQ=255)
                   cTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGgaga                                                 >  1:6112330/1‑72 (MQ=255)
                    tttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGaa                                                >  2:4890297/1‑72 (MQ=255)
                    tttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGaa                                                <  2:7437271/72‑1 (MQ=255)
                     ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGaaa                                               <  2:4586313/72‑1 (MQ=255)
                       aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAaata                                             >  2:4741099/1‑72 (MQ=255)
                       aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAaata                                             <  1:4984595/72‑1 (MQ=255)
                       aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAaata                                             >  2:5784185/1‑72 (MQ=255)
                       aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAaata                                             >  2:1968246/1‑72 (MQ=255)
                        aaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGagaaagaa                                            <  1:6438020/72‑1 (MQ=255)
                        aaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAaataa                                            <  2:3835048/72‑1 (MQ=255)
                        aaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAaataa                                            >  1:168100/1‑72 (MQ=255)
                        aaaCCAAGACGGGCATAGGTCATGTCGACATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAaataa                                            >  2:4101996/1‑72 (MQ=255)
                         aaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAaataat                                           >  1:5254698/1‑72 (MQ=255)
                           ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGc                                         >  2:2734116/1‑72 (MQ=255)
                             aaGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAg                                       >  1:1847470/1‑72 (MQ=255)
                              aGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGa                                      >  2:1789782/1‑72 (MQ=255)
                              aGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGa                                      <  2:5201901/72‑1 (MQ=255)
                               gACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGAc                                     >  2:6214245/1‑72 (MQ=255)
                               gACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGAc                                     >  2:6474829/1‑72 (MQ=255)
                                aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGAcc                                    >  2:6663771/1‑72 (MQ=255)
                                aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGAcc                                    <  2:6684218/72‑1 (MQ=255)
                                aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGAcc                                    >  1:6171750/1‑72 (MQ=255)
                                 cGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCg                                   >  2:5945126/1‑72 (MQ=255)
                                 cGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCg                                   >  2:6363216/1‑72 (MQ=255)
                                  gggCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGa                                  >  1:3032805/1‑72 (MQ=255)
                                  gggCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGa                                  <  2:4312920/72‑1 (MQ=255)
                                   ggCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGAc                                 <  1:2159285/72‑1 (MQ=255)
                                    gCATAGTTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACa                                <  1:667079/72‑1 (MQ=255)
                                     cATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAg                               <  1:7038552/72‑1 (MQ=255)
                                     cATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAg                               <  1:3640543/72‑1 (MQ=255)
                                     cATAGGTCATGTCGCCATTGCCACCGTAACTNTTTTCANCGTTACCCTTGGAGAAATAATGCAGACCGACAg                               >  1:443099/1‑72 (MQ=255)
                                      aTAGGTCATGTCGCCATTGCCACNGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGc                              >  1:4950902/1‑72 (MQ=255)
                                      aTAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGc                              <  2:5948550/72‑1 (MQ=255)
                                       tAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCt                             <  2:5547879/72‑1 (MQ=255)
                                         gctCATCTCGCCAGTGCCACCGTAACTTTTTTCGCCGTTACCCTGGGAGAAATAATGCAGGCCGACAGCttt                           <  2:7177466/70‑1 (MQ=255)
                                          gTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTa                          >  1:3437094/1‑72 (MQ=255)
                                           tCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAc                         >  2:6140874/1‑72 (MQ=255)
                                           tCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAc                         <  1:3873829/72‑1 (MQ=255)
                                            cATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                        >  2:2923948/1‑72 (MQ=255)
                                            cATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                        >  2:3643890/1‑72 (MQ=255)
                                            cATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTAcc                        >  2:2534175/1‑72 (MQ=255)
                                               gTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTa                     <  2:3174423/72‑1 (MQ=255)
                                                tCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTAc                    >  2:3264074/1‑72 (MQ=255)
                                                tCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTAc                    >  1:7637688/1‑72 (MQ=255)
                                                tCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTAc                    <  2:5432834/72‑1 (MQ=255)
                                                 cGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACa                   >  1:142721/1‑72 (MQ=255)
                                                  gCCGTTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAg                  <  1:582974/72‑1 (MQ=255)
                                                  gCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAg                  >  2:7402841/1‑72 (MQ=255)
                                                  gCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAg                  <  2:2003491/72‑1 (MQ=255)
                                                   cctttGCCACCGTAACTTNTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGa                 <  1:2996664/69‑1 (MQ=255)
                                                   ccATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGa                 <  1:5734519/72‑1 (MQ=255)
                                                    cATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGAt                <  2:1425483/72‑1 (MQ=255)
                                                    cATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGAt                >  2:6588041/1‑72 (MQ=255)
                                                      ttGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCt              <  1:3077628/72‑1 (MQ=255)
                                                      ttGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCt              <  1:1598990/72‑1 (MQ=255)
                                                      ttGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCt              <  2:5706702/72‑1 (MQ=255)
                                                       tGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTa             >  2:1559097/1‑72 (MQ=255)
                                                       tGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTa             <  2:4836626/72‑1 (MQ=255)
                                                       tGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTa             <  2:3410077/72‑1 (MQ=255)
                                                       tGCCACCGTAACTTTTTTCACCGNTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTa             >  1:4969302/1‑72 (MQ=255)
                                                         ccACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACt           <  1:1794384/72‑1 (MQ=255)
                                                         ccACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACt           >  1:2557621/1‑72 (MQ=255)
                                                         ccACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACt           >  1:2605782/1‑72 (MQ=255)
                                                           aCCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACttt         <  1:1171314/72‑1 (MQ=255)
                                                             cGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTttgt       <  1:7001363/72‑1 (MQ=255)
                                                             cGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTttgt       >  2:5350501/1‑72 (MQ=255)
                                                              gTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTttgtt      <  1:2627122/72‑1 (MQ=255)
                                                              gTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTttgtt      <  1:92697/72‑1 (MQ=255)
                                                              gTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTttgtt      <  1:4154034/72‑1 (MQ=255)
                                                              gTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTttgtt      <  2:5291106/72‑1 (MQ=255)
                                                               tAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTttgttg     >  1:3741166/1‑72 (MQ=255)
                                                               tAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTttgttg     <  2:1852479/72‑1 (MQ=255)
                                                               tAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTttgtcg     <  1:430275/72‑3 (MQ=255)
                                                                aaCTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGc    >  1:4340551/1‑72 (MQ=255)
                                                                 aCTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGcc   <  2:3962639/72‑1 (MQ=255)
                                                                    |                                                                     
GAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCA  >  NC_012967/1003851‑1003988

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTAT  >  NC_012967/1003844‑1003998
                                                                           |                                                                               
gtctTCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTT                                                                             >  SRR333730.7977276/5‑80 (MQ=60)
 gtcTCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTC                                                                            >  SRR333730.3363863/4‑80 (MQ=60)
    gtctAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACC                                                                         >  SRR333730.1419511/5‑80 (MQ=60)
     CGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCAgac                                                                        <  SRR333730.2634900/80‑4 (MQ=60)
        AATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCagac                                                                     <  SRR333730.896221/80‑5 (MQ=60)
        AATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCagac                                                                     <  SRR333730.2693590/80‑5 (MQ=60)
         gtctATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTAC                                                                    >  SRR333730.2610256/5‑80 (MQ=60)
         ATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGAATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGagac                                                                    <  SRR333730.1667859/80‑5 (MQ=60)
          gtctTTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACC                                                                   >  SRR333730.2867065/5‑80 (MQ=60)
          TTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTagac                                                                   <  SRR333730.4159515/80‑5 (MQ=60)
           TGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTAgac                                                                  <  SRR333730.4597141/80‑4 (MQ=60)
           TGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTAgac                                                                  <  SRR333730.6076376/80‑4 (MQ=60)
            gtcTTGAGTTTCCCCTTTAAAACCAAGACNGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCT                                                                 >  SRR333730.941652/4‑80 (MQ=60)
            GATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTAagac                                                                 <  SRR333730.1672542/80‑5 (MQ=60)
             gtcTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTT                                                                >  SRR333730.274385/4‑80 (MQ=60)
              gtctAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTG                                                               >  SRR333730.5271216/5‑80 (MQ=60)
              TTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACAGTAACTTTTTTCACCGTTACCagac                                                               <  SRR333730.6092207/80‑5 (MQ=60)
               gtctGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGG                                                              >  SRR333730.4361668/5‑80 (MQ=60)
                gtctTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGA                                                             >  SRR333730.1917043/5‑80 (MQ=60)
                  gtcTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTCGGAGA                                                           >  SRR333730.47915/4‑80 (MQ=60)
                  AGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGagac                                                           <  SRR333730.4205609/80‑5 (MQ=60)
                  AGTCTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGagac                                                           <  SRR333730.3663101/80‑5 (MQ=60)
                   ccTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAC                                                          <  SRR333730.5413148/78‑1 (MQ=60)
                    gtctCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAAGTTTCTTCACCGTTACCCTTGGAGAAA                                                         >  SRR333730.2850029/5‑80 (MQ=60)
                    TTTCCCCTTGAAGACCAAGCCATGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTATCCTTGGAagac                                                         <  SRR333730.5944785/80‑5 (MQ=48)
                     TTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAgac                                                        <  SRR333730.168645/80‑4 (MQ=60)
                      gtctCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATA                                                       >  SRR333730.6697708/5‑80 (MQ=60)
                        gtCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAAT                                                     >  SRR333730.1554532/3‑80 (MQ=60)
                         CCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAgac                                                    <  SRR333730.2183157/80‑4 (MQ=60)
                          gtcTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGC                                                   >  SRR333730.7296328/4‑80 (MQ=60)
                          CTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAAgac                                                   <  SRR333730.5471825/80‑4 (MQ=60)
                           gtctAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCA                                                  >  SRR333730.1761715/5‑80 (MQ=60)
                           gtctAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCA                                                  >  SRR333730.186875/5‑80 (MQ=60)
                           TTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAAgaaagac                                                  <  SRR333730.7050404/80‑8 (MQ=60)
                            gtctAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAG                                                 >  SRR333730.5765888/5‑80 (MQ=60)
                            TTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAAT‑‑AGAC                                               <  SRR333730.8366961/80‑1 (MQ=60)
                             TAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGAGAC                                               <  SRR333730.3681269/80‑1 (MQ=60)
                              AAAACCAGGCCGGGCATAGGTCCTGTAGCCCTTGCCACAGTAACTTTTTTAACCGTTACCCTTGGAGAAATAATGCAGAC                                               <  SRR333730.705753/80‑1 (MQ=39)
                                gtctAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCG                                             >  SRR333730.2038717/5‑80 (MQ=60)
                                 ACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGAagac                                            <  SRR333730.6568400/80‑5 (MQ=60)
                                   gtctACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACA                                          >  SRR333730.6762111/5‑80 (MQ=60)
                                    gtctCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAG                                         >  SRR333730.616877/5‑80 (MQ=60)
                                     gtctGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAag                                        >  SRR333730.3339265/5‑78 (MQ=60)
                                      GACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGac                                       <  SRR333730.2378427/80‑3 (MQ=60)
                                       ACGGGCATAGTTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAagac                                      <  SRR333730.734504/80‑5 (MQ=60)
                                       ACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACagacaagac                                      <  SRR333730.6858126/80‑10 (MQ=60)
                                        gtctCATAGGTCATGTCGCCATTGCCACCGTAACTNTTTTCANCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTT                                     >  SRR333730.489228/5‑80 (MQ=60)
                                        CGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGagac                                     <  SRR333730.3999193/80‑5 (MQ=60)
                                        CGGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGagac                                     <  SRR333730.7704458/80‑5 (MQ=60)
                                         gtctATAGGTCATGTCGCCATTGCCACNGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTA                                    >  SRR333730.5434702/5‑80 (MQ=60)
                                         GGGCATAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCagac                                    <  SRR333730.4459122/80‑5 (MQ=60)
                                             gtctGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGT                                >  SRR333730.3777721/5‑80 (MQ=60)
                                              TAGGTCATGTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACagac                               <  SRR333730.4255098/80‑5 (MQ=60)
                                                   gtctTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGAT                          >  SRR333730.8351056/5‑80 (MQ=60)
                                                    gtcTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATC                         >  SRR333730.158947/4‑80 (MQ=60)
                                                    ATCTCNCCATTCCCACCGTAACTTTTTTCACNGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAagac                         <  SRR333730.842112/80‑5 (MQ=60)
                                                     TGTCGCCGTTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGAgac                        <  SRR333730.642758/80‑4 (MQ=60)
                                                      tccgcctTTGCCACCGTAACTTNTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGAagac                       <  SRR333730.3299707/73‑5 (MQ=60)
                                                      GTCTCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTA                       >  SRR333730.2985933/1‑80 (MQ=60)
                                                      GTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGAagac                       <  SRR333730.3364580/80‑5 (MQ=60)
                                                      GTCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGAagac                       <  SRR333730.6288582/80‑5 (MQ=60)
                                                       TCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATagac                      <  SRR333730.4111744/80‑5 (MQ=60)
                                                       TCGCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATagac                      <  SRR333730.4781029/80‑5 (MQ=60)
                                                        gtctATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACT                     >  SRR333730.4753924/5‑80 (MQ=60)
                                                        cccaattcCCACCGTAGCTTTTTTCACCGTTACCCTTGAAGAAAAAATGCAGACCGACAGCTTTACCGTACAGATCagac                     <  SRR333730.7370365/72‑5 (MQ=60)
                                                         GCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTAgac                    <  SRR333730.3388190/80‑4 (MQ=60)
                                                         GCCATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTAgac                    <  SRR333730.1764893/80‑4 (MQ=60)
                                                          gtcTTGCCACCGTAACTTTTTTCACCGNTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTT                   >  SRR333730.5454989/4‑80 (MQ=60)
                                                           CATTGCCGCCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACagac                  <  SRR333730.642609/80‑5 (MQ=60)
                                                            gtctCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGT                 >  SRR333730.2863674/5‑80 (MQ=60)
                                                            gtctCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGT                 >  SRR333730.5471281/5‑80 (MQ=60)
                                                            gtctCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGT                 >  SRR333730.2810146/5‑80 (MQ=60)
                                                            ATTGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTagac                 <  SRR333730.1979731/80‑5 (MQ=60)
                                                            ATGGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTagac                 <  SRR333730.3493707/80‑5 (MQ=60)
                                                              TGCCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTagac               <  SRR333730.1299247/80‑5 (MQ=60)
                                                               gtctCCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGC              >  SRR333730.5472800/5‑80 (MQ=60)
                                                                CCACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTagac             <  SRR333730.7663800/80‑5 (MQ=60)
                                                                 CACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTagac            <  SRR333730.2887623/80‑5 (MQ=60)
                                                                 CACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTagac            <  SRR333730.4561343/80‑5 (MQ=60)
                                                                 CACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTagac            <  SRR333730.103563/80‑5 (MQ=60)
                                                                  gtctTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCAT           >  SRR333730.4110338/5‑80 (MQ=60)
                                                                  ACCGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTcgagac           <  SRR333730.475083/80‑7 (MQ=60)
                                                                   gtcTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATC          >  SRR333730.4763748/4‑80 (MQ=60)
                                                                    CGTAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCAgac         <  SRR333730.4008464/80‑4 (MQ=60)
                                                                     GTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCAagac        <  SRR333730.183927/80‑5 (MQ=60)
                                                                      TAACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATagac       <  SRR333730.4585571/80‑5 (MQ=60)
                                                                       AACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCagac      <  SRR333730.5274648/80‑5 (MQ=60)
                                                                       AACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCagac      <  SRR333730.4874024/80‑5 (MQ=60)
                                                                        ACTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTagac     <  SRR333730.3794229/80‑5 (MQ=60)
                                                                        ACTTTTTTCACCGTTACCCTTGGAGAAATAAGGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTagac     <  SRR333730.2421354/80‑5 (MQ=60)
                                                                         CTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCaatctcagac    <  SRR333730.5872706/80‑11 (MQ=60)
                                                                         CTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTagac    <  SRR333730.5218277/80‑5 (MQ=60)
                                                                         CTTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCAATCTTagac    <  SRR333730.5872707/80‑5 (MQ=60)
                                                                          TTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTagac   <  SRR333730.3389731/80‑5 (MQ=60)
                                                                          TTTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTagac   <  SRR333730.1882833/80‑5 (MQ=60)
                                                                           TTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGagac  <  SRR333730.7755749/80‑5 (MQ=60)
                                                                           TTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGagac  <  SRR333730.5556658/80‑5 (MQ=60)
                                                                           GTCTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTAT  >  SRR333730.2093474/1‑80 (MQ=60)
                                                                           |                                                                               
CAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAGAAATAATGCAGACCGACAGCTTTACCGTACAGATCTACTTTGTTGCCATCTTTGTTAT  >  NC_012967/1003844‑1003998

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