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 F121 I0 R1 856 273.6 9081716 94.1% 8545894 150.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,543,959 G→T 85.0% A625E (GCG→GAG)  malP ← maltodextrin phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,543,9590GT85.0% 542.8 / 90.0 241A625E (GCG→GAG) malPmaltodextrin phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (19/17);  new base T (99/101);  total (123/118)
Fisher's exact test for biased strand distribution p-value = 8.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

ATATTTTCTTCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGTAGCCCGG  >  CP009273/3543812‑3544104
                                                                                                                                                   |                                                                                                                                                 
atatTTTCTTCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGACGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATccgcctcc                                                                                                                                                 <  2:1729285/150‑1 (MQ=255)
 tatTTTCTTCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATccgccgccg                                                                                                                                                >  2:1359244/1‑150 (MQ=255)
     ttcttcACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGAt                                                                                                                                            >  2:3785275/1‑150 (MQ=255)
      tcttcACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATc                                                                                                                                           >  1:3890139/1‑150 (MQ=255)
       cttcACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCa                                                                                                                                          >  2:2705739/1‑150 (MQ=255)
       cttcACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCa                                                                                                                                          <  2:4261645/150‑1 (MQ=255)
         tcACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTTTTCGGTGATATCACCCACCGGGCTCAGt                                                                                                                                        >  1:8173/1‑150 (MQ=255)
         tcACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGt                                                                                                                                        >  1:4056538/1‑150 (MQ=255)
           aCCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGttt                                                                                                                                      <  1:3836352/150‑1 (MQ=255)
            ccGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTAGGAGATATCCGCCACCCGGATCAGtttt                                                                                                                                     >  1:3126846/1‑150 (MQ=255)
             cGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGttttt                                                                                                                                    <  1:1509665/150‑1 (MQ=255)
             cGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGttttt                                                                                                                                    <  1:3139254/150‑1 (MQ=255)
              gACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTc                                                                                                                                   >  1:4033027/1‑150 (MQ=255)
              gACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTc                                                                                                                                   >  2:3919012/1‑150 (MQ=255)
              gACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCGTCTTTACCTGCAGTCGAAATTTTTTAGGAGATATCCGCCTCCGGGATCAGttttta                                                                                                                                   >  1:3458319/1‑149 (MQ=255)
               aCTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTcc                                                                                                                                  <  2:2067141/150‑1 (MQ=255)
               aCTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTcc                                                                                                                                  <  2:492854/150‑1 (MQ=255)
                cTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTccg                                                                                                                                 <  2:1252729/150‑1 (MQ=255)
                cTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTccg                                                                                                                                 <  2:666884/150‑1 (MQ=255)
                 tttCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTAGGAGATATCCGCCACCGGGCTCAGTTTTTacgc                                                                                                                                >  1:4327451/1‑150 (MQ=255)
                     tcGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGCTCAGTTTTTccgccgcc                                                                                                                            >  1:4380058/1‑150 (MQ=255)
                     tcGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTccgccgcc                                                                                                                            >  1:4379975/1‑150 (MQ=255)
                      cGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTccgccgccg                                                                                                                           >  1:1354066/1‑150 (MQ=255)
                      cGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTccgccgccg                                                                                                                           >  1:247876/1‑150 (MQ=255)
                       ggCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGa                                                                                                                          <  1:279361/150‑1 (MQ=255)
                          gATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAAc                                                                                                                       >  2:4367566/1‑150 (MQ=255)
                          gATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAAc                                                                                                                       >  2:288915/1‑150 (MQ=255)
                                aaCGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCaata                                                                                                                 <  1:1171908/150‑1 (MQ=255)
                                aaCGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCGTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTGTTTCCGCCGCCGAAACGCacta                                                                                                                 >  2:3932980/1‑150 (MQ=255)
                                  cGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCaataat                                                                                                               <  2:3737144/150‑1 (MQ=255)
                                   gTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATc                                                                                                              <  2:3830444/150‑1 (MQ=255)
                                   gTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATc                                                                                                              <  2:2133653/150‑1 (MQ=255)
                                     tCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCg                                                                                                            >  2:1765998/1‑150 (MQ=255)
                                      cGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCgg                                                                                                           <  2:522103/150‑1 (MQ=255)
                                          ccATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAgg                                                                                                       <  2:4134496/150‑1 (MQ=255)
                                              ccAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAaca                                                                                                   >  2:2239044/1‑150 (MQ=255)
                                              ccAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAaca                                                                                                   >  2:2239394/1‑150 (MQ=255)
                                                 gCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAAcacca                                                                                                <  2:219107/150‑1 (MQ=255)
                                                  cGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAAcaccac                                                                                               >  1:3361164/1‑150 (MQ=255)
                                                  cGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAAcaccac                                                                                               <  2:4086349/150‑1 (MQ=255)
                                                   gTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACaccacc                                                                                              >  2:4416857/1‑150 (MQ=255)
                                                   gTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACaccacc                                                                                              <  2:2750653/150‑1 (MQ=255)
                                                   gTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACaccacc                                                                                              >  1:1099794/1‑150 (MQ=255)
                                                   gTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACaccacc                                                                                              >  2:292264/1‑150 (MQ=255)
                                                     aCCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGCTCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCtt                                                                                            >  2:3135553/1‑150 (MQ=255)
                                                     aCCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCtt                                                                                            <  2:3205061/150‑1 (MQ=255)
                                                     aCCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCtt                                                                                            <  2:3205223/150‑1 (MQ=255)
                                                       cGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCa                                                                                          <  1:2811478/150‑1 (MQ=255)
                                                           aGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACtt                                                                                      <  2:1631388/150‑1 (MQ=255)
                                                              aaGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATc                                                                                   >  2:4463010/1‑150 (MQ=255)
                                                              aaGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATc                                                                                   <  1:4497891/150‑1 (MQ=255)
                                                               aGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCg                                                                                  <  1:2536445/150‑1 (MQ=255)
                                                               aGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCg                                                                                  <  1:2571224/150‑1 (MQ=255)
                                                                  gcACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGCTCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCATCAACTTATCGCCa                                                                               >  1:946920/1‑150 (MQ=255)
                                                                  gcACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATACGGCAGGAACACCACCTTCAACTTATCGCCa                                                                               >  1:3694806/1‑150 (MQ=255)
                                                                  gcACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTACGCCGCCGAAACGCAATAATACGGCAGGAACACCACCTTCAACTTATAGCCa                                                                               >  1:951362/1‑150 (MQ=255)
                                                                     ccATTGAGCGCCAGTGTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCCAAACGCAATAATCCGCCAGGAACCCCACCTTCAACATATAGCCAAcc                                                                            >  2:2489647/1‑150 (MQ=255)
                                                                        ttGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGc                                                                         <  2:1909467/150‑1 (MQ=255)
                                                                         tGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTTTTGGGAGATATACGCCCCCGGGATCAGTTTTTACGCCGCCGAAAAGCCAAAAAACGGCAGGAACACCACCATCAACATAAAGCCAAcccccg                                                                        >  1:3784375/1‑146 (MQ=255)
                                                                         tGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCg                                                                        >  2:1580035/1‑150 (MQ=255)
                                                                         tGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCg                                                                        >  2:3329619/1‑150 (MQ=255)
                                                                                cAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCgtt                                                                 <  1:1593660/150‑1 (MQ=255)
                                                                                  gTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCgttgt                                                               >  1:1394481/1‑150 (MQ=255)
                                                                                  gTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCgttgt                                                               <  2:1320764/150‑1 (MQ=255)
                                                                                  gTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCATCAACTTATCGCCAACCAGCGGATCgttgt                                                               >  1:1770701/1‑150 (MQ=255)
                                                                                    ttCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGttgttg                                                             >  1:2852029/1‑150 (MQ=255)
                                                                                    ttCATATTGCCGGTACCGGAAGCGTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGttgttg                                                             <  2:3439284/150‑1 (MQ=255)
                                                                                      cATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGAt                                                           <  2:4033027/150‑1 (MQ=255)
                                                                                       atatTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATc                                                          <  1:128818/150‑1 (MQ=255)
                                                                                       atatTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATc                                                          <  2:2208722/150‑1 (MQ=255)
                                                                                       atatTGCCGGTACCGGAAGCTTCGTTACCATCAGTCGAAAATTGTTCGGAGATATACGCCTCCGGGGTCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCATGAACACCACCTTCAACTTATCGCCAAACAGCGGCTCGTTGTTGATc                                                          >  2:3409187/1‑150 (MQ=255)
                                                                                         atTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAAATTGTTCGGAGATATACGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCATCAAAATATAGCCAACCAGCGGATCGTTGTTGATCAc                                                        >  1:3585647/1‑150 (MQ=255)
                                                                                             ccGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCa                                                    <  2:4056538/150‑1 (MQ=255)
                                                                                                gTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGcc                                                 <  1:3785275/150‑1 (MQ=255)
                                                                                                  aCCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTACGCCGCCGAAACGCAAAAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCAc                                               >  1:1024833/1‑150 (MQ=255)
                                                                                                  aCCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCAc                                               <  1:198614/150‑1 (MQ=255)
                                                                                                   ccGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACt                                              <  1:697804/150‑1 (MQ=255)
                                                                                                    cGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACtt                                             >  2:1842533/1‑150 (MQ=255)
                                                                                                    cGGAAGCTTCTGTACCTGCAGGCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACtt                                             <  2:433590/150‑1 (MQ=255)
                                                                                                     ggAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATAACGTCAGCCACttt                                            >  1:4316846/1‑150 (MQ=255)
                                                                                                          cttcttTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGa                                       <  1:1974079/150‑1 (MQ=255)
                                                                                                             cttTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCg                                    <  1:559034/150‑1 (MQ=255)
                                                                                                                tACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCaa                                 >  1:2715949/1‑150 (MQ=255)
                                                                                                                tACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCaa                                 <  2:3477805/150‑1 (MQ=255)
                                                                                                                   cTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGa                              <  1:2903255/150‑1 (MQ=255)
                                                                                                                          cGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATAtt                       <  2:4441986/150‑1 (MQ=255)
                                                                                                                           gAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATAttc                      <  2:4361272/150‑1 (MQ=255)
                                                                                                                              aTTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATAttcttc                   <  2:2530702/150‑1 (MQ=255)
                                                                                                                              aTTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCTCGTCAGCCACTTTGTTGATCGCAAAGATAATAttcttc                   <  2:4496935/150‑1 (MQ=255)
                                                                                                                              aTTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATAttcttc                   >  1:2216662/1‑150 (MQ=255)
                                                                                                                               tttGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCg                  <  1:2489647/150‑1 (MQ=255)
                                                                                                                                ttGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTAGc                 <  1:3322564/150‑1 (MQ=255)
                                                                                                                                   ttCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAg              <  2:1979964/150‑1 (MQ=255)
                                                                                                                                     cGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGgt            <  2:3126846/150‑1 (MQ=255)
                                                                                                                                     cGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGgt            <  2:861400/150‑1 (MQ=255)
                                                                                                                                      ggAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGgta           >  1:2780235/1‑150 (MQ=255)
                                                                                                                                        agaTATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGgtagt         >  1:722934/1‑150 (MQ=255)
                                                                                                                                        agaTATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGgtagt         <  1:7934/150‑1 (MQ=255)
                                                                                                                                           tatCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGTAGc      >  1:162209/1‑150 (MQ=255)
                                                                                                                                               cgcctccgGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGTAGCCCgg  >  2:1789406/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATATTTTCTTCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGTAGCCCGG  >  CP009273/3543812‑3544104

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTTTCTTCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGTAGCCCGGTGCC  >  CP009273/3543815‑3544108
                                                                                                                                                |                                                                                                                                                     
TTTTCTTCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTAACTGCAGTCGAAAATTGTTCTGAGATATACGCCTCCGGG                                                                                                                                                  >  SRR19973080.797351/1‑150 (MQ=60)
 TTTCTTCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGA                                                                                                                                                 >  SRR19973080.1120993/1‑150 (MQ=60)
   TCTTCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATC                                                                                                                                               >  SRR19973080.3911580/1‑150 (MQ=60)
      TCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTTTTCGGTGATATCACCCACCGGGCTCAGT                                                                                                                                            >  SRR19973080.8268/1‑150 (MQ=60)
      TCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGT                                                                                                                                            >  SRR19973080.4078709/1‑150 (MQ=60)
      TCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTAGGAGATATACGCCTCCGGGATCAGT                                                                                                                                            >  SRR19973080.4077455/1‑150 (MQ=60)
      TCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCAGCAGTCGAAATTTGTTCGGAGATATACGCCTCCGGGCTCAGT                                                                                                                                            >  SRR19973080.3189935/1‑150 (MQ=60)
        ACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTT                                                                                                                                          <  SRR19973080.3857424/150‑1 (MQ=60)
         CCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTAGGAGATATCCGCCACCCGGATCAGTTTT                                                                                                                                         >  SRR19973080.3144054/1‑150 (MQ=60)
         CCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTATTTACCAGCAGTCGAAATTTTTTAGGAGATATCCCCCTCCGGGCTCAGTTTT                                                                                                                                         >  SRR19973080.866644/1‑150 (MQ=60)
          CGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTT                                                                                                                                        <  SRR19973080.1518295/150‑1 (MQ=60)
          CGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTT                                                                                                                                        <  SRR19973080.3156484/150‑1 (MQ=60)
          CGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTT                                                                                                                                        <  SRR19973080.3156533/150‑1 (MQ=60)
           GACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTC                                                                                                                                       >  SRR19973080.4055115/1‑150 (MQ=60)
           GACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTC                                                                                                                                       >  SRR19973080.3358694/1‑150 (MQ=60)
           GACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCGTCTTTACCTGCAGTCGAAATTTTTTAGGAGATATCCGCCTCCGGGATCAGTTTTTA                                                                                                                                       >  SRR19973080.3477264/1‑150 (MQ=60)
             CTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCG                                                                                                                                     <  SRR19973080.2836397/150‑1 (MQ=60)
              TTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATACGCCTCCGGGATCAGTTTTTCCGC                                                                                                                                    >  SRR19973080.4350952/1‑150 (MQ=60)
              TTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTAGGAGATATCCGCCACCGGGCTCAGTTTTTACGC                                                                                                                                    >  SRR19973080.4351088/1‑150 (MQ=60)
               TTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCC                                                                                                                                   >  SRR19973080.2982153/1‑150 (MQ=60)
                  TCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGCTCAGTTTTTCCGCCGCC                                                                                                                                >  SRR19973080.4403931/1‑150 (MQ=60)
                  TCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCC                                                                                                                                >  SRR19973080.4403848/1‑150 (MQ=60)
                  TCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCC                                                                                                                                >  SRR19973080.4403875/1‑150 (MQ=60)
                   CGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCG                                                                                                                               >  SRR19973080.1363211/1‑150 (MQ=60)
                   CGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCG                                                                                                                               >  SRR19973080.249673/1‑150 (MQ=60)
                   CGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCG                                                                                                                               >  SRR19973080.249581/1‑150 (MQ=60)
                   CGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCG                                                                                                                               >  SRR19973080.1362184/1‑150 (MQ=60)
                    GGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGA                                                                                                                              <  SRR19973080.281189/150‑1 (MQ=60)
                    GGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGA                                                                                                                              <  SRR19973080.280384/150‑1 (MQ=60)
                    GGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGA                                                                                                                              <  SRR19973080.280308/150‑1 (MQ=60)
                       GATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTAGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAAC                                                                                                                           >  SRR19973080.2817416/1‑150 (MQ=60)
                             AACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATA                                                                                                                     <  SRR19973080.1178885/150‑1 (MQ=60)
                                GTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATC                                                                                                                  >  SRR19973080.3952258/1‑150 (MQ=60)
                                GTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATC                                                                                                                  >  SRR19973080.2589656/1‑150 (MQ=60)
                                GTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCGTCTTTACCTGCCGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTACGCCGCCGCAACGCCAGAATC                                                                                                                  >  SRR19973080.3862536/1‑150 (MQ=60)
                                       CCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATACGCCGCCGGGCTCAGTTTTTCCGCCGCCGAAACGCAAAAAACCCGCAGG                                                                                                           >  SRR19973080.3107799/1‑150 (MQ=60)
                                       CCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTAACTGCAGTCGAAATTTGTTCGGAGATATACGCCGCCGGGATCAGTTTTTACGCCGCCGAAACCCCAAAATCCGGCAGG                                                                                                           >  SRR19973080.2309997/1‑150 (MQ=60)
                                            CAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACAC                                                                                                      <  SRR19973080.289905/150‑1 (MQ=60)
                                               CGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCAC                                                                                                   >  SRR19973080.3379403/1‑150 (MQ=60)
                                               CGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAAAAATACGGCAGGAACACCAC                                                                                                   >  SRR19973080.97965/1‑150 (MQ=60)
                                                GTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACC                                                                                                  >  SRR19973080.1106123/1‑150 (MQ=60)
                                                  ACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTT                                                                                                >  SRR19973080.380517/1‑150 (MQ=60)
                                                    CGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCA                                                                                              <  SRR19973080.2826962/150‑1 (MQ=60)
                                                    CGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCA                                                                                              >  SRR19973080.606095/1‑150 (MQ=60)
                                                         GTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTA                                                                                         >  SRR19973080.168053/1‑150 (MQ=60)
                                                         GTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCATCAACTTA                                                                                         >  SRR19973080.168043/1‑150 (MQ=60)
                                                           AAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATC                                                                                       <  SRR19973080.4522190/150‑1 (MQ=60)
                                                           AAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATACGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCATCAACatata                                                                                       >  SRR19973080.1593372/1‑145 (MQ=60)
                                                           AAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTAGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATACGGCAGGACCACCACCATCAACTTATC                                                                                       >  SRR19973080.4084683/1‑150 (MQ=60)
                                                            AGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCG                                                                                      <  SRR19973080.2550363/150‑1 (MQ=60)
                                                            AGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCG                                                                                      <  SRR19973080.878857/150‑1 (MQ=60)
                                                            AGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCG                                                                                      <  SRR19973080.2585262/150‑1 (MQ=60)
                                                            AGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCG                                                                                      <  SRR19973080.2586483/150‑1 (MQ=60)
                                                               GCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGCTCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCATCAACTTATCGCCA                                                                                   >  SRR19973080.952683/1‑150 (MQ=60)
                                                               GCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATACGGCAGGAACACCACCTTCAACTTATCGCCA                                                                                   >  SRR19973080.3715111/1‑150 (MQ=60)
                                                               GCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTACGCCGCCGAAACGCAATAATACGGCAGGAACACCACCTTCAACTTATAGCCA                                                                                   >  SRR19973080.957153/1‑150 (MQ=60)
                                                               GCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCAGCAGTCGAAAATTGTTCGGAGATATACGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCGTCAACATATAGCCA                                                                                   >  SRR19973080.952797/1‑150 (MQ=60)
                                                                    ATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTAATTAACTGCAGTAGAAAATTGTTAGGAGATATACCCCACCCGGATCAGTTTTTCCGCCGCCGAAAAGCAATAAAACCGCAGGAACACCACCataaacatatagccaaccac                                                                              >  SRR19973080.3850712/1‑130 (MQ=59)
                                                                      TGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTTTTGGGAGATATACGCCCCCGGGATCAGTTTTTACGCCGCCGAAAAGCCAAAAAACGGCAGGAACACCACCATCAACataaagccaacccccg                                                                            >  SRR19973080.3804948/1‑134 (MQ=60)
                                                                      TGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTATTTACCAGCAGTCGCAATTTTTTAGGAGATATACGCCGCCGGGATCAGTTTTTACGCCGCCCAAACGCAAAAATACCGCAGGAACACCACCATCAACATATAGCCAAACACCG                                                                            >  SRR19973080.3805496/1‑150 (MQ=60)
                                                                      TGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTATTTAACAGCAGTCGAAATTTTTTATGAGATATACGCCGCCCGGATCAGTTTTTACGCCGCCGCAACGCAAAAAAACGGCAGGAAAACCACCATCAACATATCGCCAACCACCG                                                                            >  SRR19973080.3804909/1‑150 (MQ=60)
                                                                           GCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCG                                                                       <  SRR19973080.4403183/150‑1 (MQ=60)
                                                                           GCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCG                                                                       <  SRR19973080.4401582/150‑1 (MQ=60)
                                                                            CCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGT                                                                      <  SRR19973080.3188412/150‑1 (MQ=60)
                                                                             CAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTT                                                                     <  SRR19973080.1602976/150‑1 (MQ=60)
                                                                             CAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTT                                                                     <  SRR19973080.1602982/150‑1 (MQ=60)
                                                                              AGTTTCATATTGCCGGTACCGGAAGCTTATTTAACTGCAGTAGAAAATTGTTAGGCTATATCCGCCGACGGGCTAAATTTTTACCCCGCCGAAAAGCAAAAATCCGGCAGGAACACAACCTTCAAaaaatagccaaccagccgatcgtta                                                                    >  SRR19973080.1275627/1‑125 (MQ=35)
                                                                               GTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGT                                                                   >  SRR19973080.1402762/1‑150 (MQ=60)
                                                                               GTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCATCAACTTATCGCCAACCAGCGGATCGTTGT                                                                   >  SRR19973080.1780862/1‑150 (MQ=60)
                                                                                 TTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTG                                                                 >  SRR19973080.2867674/1‑150 (MQ=60)
                                                                                    ATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATC                                                              <  SRR19973080.129585/150‑1 (MQ=60)
                                                                                    ATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATC                                                              <  SRR19973080.129909/150‑1 (MQ=60)
                                                                                     TATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATACGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAAAAATCCCGCAGGAACCCCACCATCAACTTATCGCCAACCAGCGGATCGTTGTTGATCA                                                             >  SRR19973080.578711/1‑150 (MQ=60)
                                                                                     TATTGCCGGTACCGGAAGCTTATTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAAAGCAAAAATCAGGCAGGAAAACCACCATAAACTTATCGCCAACCAGCGGATCGTTGTTGATCA                                                             >  SRR19973080.2855132/1‑150 (MQ=60)
                                                                                      ATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAAATTGTTCGGAGATATACGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCATCAAAATATAGCCAACCAGCGGATCGTTGTTGATCAC                                                            >  SRR19973080.3605472/1‑150 (MQ=60)
                                                                                            GGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTAGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGC                                                      >  SRR19973080.2567258/1‑150 (MQ=60)
                                                                                             GTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCC                                                     <  SRR19973080.4329162/150‑1 (MQ=60)
                                                                                             GTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCC                                                     <  SRR19973080.3805857/150‑1 (MQ=60)
                                                                                             GTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCC                                                     <  SRR19973080.3806304/150‑1 (MQ=60)
                                                                                               ACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTACGCCGCCGAAACGCAAAAATCCGGCAGGAACACCACCATCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCAC                                                   >  SRR19973080.1030924/1‑150 (MQ=60)
                                                                                               ACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCAC                                                   <  SRR19973080.1367382/150‑1 (MQ=60)
                                                                                               ACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCAC                                                   <  SRR19973080.200098/150‑1 (MQ=60)
                                                                                                CCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACT                                                  <  SRR19973080.701892/150‑1 (MQ=60)
                                                                                                CCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACT                                                  <  SRR19973080.703604/150‑1 (MQ=60)
                                                                                                  GGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATAACGTCAGCCACTTT                                                >  SRR19973080.4340426/1‑150 (MQ=60)
                                                                                                     AGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTT                                             <  SRR19973080.3526621/150‑1 (MQ=60)
                                                                                                     AGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTT                                             >  SRR19973080.4408679/1‑150 (MQ=60)
                                                                                                     AGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTT                                             <  SRR19973080.1275572/150‑1 (MQ=60)
                                                                                                      GCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTG                                            >  SRR19973080.2545133/1‑150 (MQ=60)
                                                                                                      GCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTG                                            <  SRR19973080.2382204/150‑1 (MQ=60)
                                                                                                       CTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGA                                           <  SRR19973080.1985600/150‑1 (MQ=60)
                                                                                                       CTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGA                                           <  SRR19973080.1985594/150‑1 (MQ=60)
                                                                                                       CTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGA                                           <  SRR19973080.1985295/150‑1 (MQ=60)
                                                                                                         TCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATC                                         <  SRR19973080.1457859/150‑1 (MQ=60)
                                                                                                          CTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCG                                        <  SRR19973080.562556/150‑1 (MQ=60)
                                                                                                            TTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCA                                      >  SRR19973080.812611/1‑150 (MQ=60)
                                                                                                            TTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCA                                      >  SRR19973080.339990/1‑150 (MQ=60)
                                                                                                             TACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAA                                     >  SRR19973080.2730844/1‑150 (MQ=60)
                                                                                                             TACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAA                                     >  SRR19973080.451615/1‑150 (MQ=60)
                                                                                                                CTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGA                                  <  SRR19973080.2919100/150‑1 (MQ=60)
                                                                                                                 TGCAGTCGAAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGAT                                 <  SRR19973080.2531539/150‑1 (MQ=60)
                                                                                                                         AAATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCT                         >  SRR19973080.1036256/1‑150 (MQ=60)
                                                                                                                          AATTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTT                        >  SRR19973080.1926480/1‑150 (MQ=60)
                                                                                                                           ATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTC                       >  SRR19973080.2228908/1‑150 (MQ=60)
                                                                                                                            TTTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCG                      <  SRR19973080.2503393/150‑1 (MQ=60)
                                                                                                                             TTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGC                     <  SRR19973080.2754613/150‑1 (MQ=60)
                                                                                                                             TTGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTAGC                     <  SRR19973080.3340696/150‑1 (MQ=60)
                                                                                                                              TGTTCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCC                    <  SRR19973080.3659623/150‑1 (MQ=60)
                                                                                                                                 TCGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGG                 <  SRR19973080.3550529/150‑1 (MQ=60)
                                                                                                                                  CGGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGT                >  SRR19973080.1563504/1‑150 (MQ=60)
                                                                                                                                   GGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTA               >  SRR19973080.2795575/1‑150 (MQ=60)
                                                                                                                                   GGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTA               <  SRR19973080.394793/150‑1 (MQ=60)
                                                                                                                                   GGAGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTA               >  SRR19973080.521774/1‑150 (MQ=60)
                                                                                                                                     AGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGT             >  SRR19973080.727827/1‑150 (MQ=60)
                                                                                                                                     AGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGT             >  SRR19973080.727088/1‑150 (MQ=60)
                                                                                                                                     AGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGT             <  SRR19973080.8027/150‑1 (MQ=60)
                                                                                                                                     AGATATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGT             >  SRR19973080.2132357/1‑150 (MQ=60)
                                                                                                                                        TATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGTAGC          >  SRR19973080.163503/1‑150 (MQ=60)
                                                                                                                                         ATCCGCCTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGTAGCC         <  SRR19973080.134038/150‑1 (MQ=60)
                                                                                                                                               CTCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGTAGCCCGGTGC   >  SRR19973080.659026/1‑150 (MQ=60)
                                                                                                                                                TCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGTAGCCCGGTGCC  >  SRR19973080.2977693/1‑150 (MQ=60)
                                                                                                                                                |                                                                                                                                                     
TTTTCTTCACCGACTTTCTCGGCGATTTCAACGTTCGCCCCATCCAGCGTACCGACAGTAAGCGCACCATTGAGCGCCAGTTTCATATTGCCGGTACCGGAAGCTTCTTTACCTGCAGTCGAAATTTGTTCGGAGATATCCGCCGCCGGGATCAGTTTTTCCGCCGCCGAAACGCAATAATCCGGCAGGAACACCACCTTCAACTTATCGCCAACCAGCGGATCGTTGTTGATCACGTCAGCCACTTTGTTGATCGCAAAGATAATATTCTTCGCCAGGTAGTAGCCCGGTGCC  >  CP009273/3543815‑3544108

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