Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 343 268.7 10873516 97.0% 10547310 150.6

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,030,927 G→A intergenic (+69/‑715) RPPX_17905 → / → RPPX_17915 23S rRNA methyltransferase/sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,030,9270GA98.1% 545.2 / NA 160intergenic (+69/‑715)RPPX_17905/RPPX_1791523S rRNA methyltransferase/sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (67/90);  minor base C (1/2);  total (68/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

CGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCA‑GGATGGGGCCTTCGGCGCGGTCCT  >  CP009974/4030779‑4031048
                                                                                                                                                    |                                                                                                                          
cGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATa                                                                                                                          >  2:388280/1‑151 (MQ=255)
 gATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAg                                                                                                                         <  1:738031/151‑1 (MQ=255)
  aTTGCCTGCAAGCAGCTGTACCACGACCGTAAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGt                                                                                                                        >  1:2294955/1‑151 (MQ=255)
   ttGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGtt                                                                                                                       <  1:4076652/151‑1 (MQ=255)
   ttGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGtt                                                                                                                       <  1:4471390/151‑1 (MQ=255)
   ttGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGtt                                                                                                                       >  2:4431200/1‑151 (MQ=255)
    tGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTc                                                                                                                      <  2:2728388/151‑1 (MQ=255)
     gCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTCCTAAGCGCGCGATAGTCCa                                                                                                                     >  1:4111059/1‑151 (MQ=255)
      ccTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAg                                                                                                                    <  1:4292481/151‑1 (MQ=255)
      ccTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAg                                                                                                                    <  2:232834/151‑1 (MQ=255)
         gCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGATCGACTTGAAGCCCCGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCg                                                                                                                 <  2:5201630/151‑1 (MQ=255)
           aaGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGAc                                                                                                               <  1:138483/151‑1 (MQ=255)
             gcagcTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGc                                                                                                             <  1:2048089/151‑1 (MQ=255)
                 cTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGAt                                                                                                         <  1:2294149/151‑1 (MQ=255)
                  tGTACCACGACCGTGAAGAAGTAACGTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGAtt                                                                                                        <  1:4006951/151‑1 (MQ=255)
                  tGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGAtt                                                                                                        <  2:843027/151‑1 (MQ=255)
                     aCCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCg                                                                                                     >  1:1985109/1‑151 (MQ=255)
                     aCCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCg                                                                                                     >  1:2551744/1‑151 (MQ=255)
                        aCGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACt                                                                                                  <  2:638520/151‑1 (MQ=255)
                          gACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTc                                                                                                <  2:2445131/151‑1 (MQ=255)
                          gACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTc                                                                                                >  2:620332/1‑151 (MQ=255)
                           aCCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCa                                                                                               >  2:1978442/1‑151 (MQ=255)
                            ccGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCag                                                                                              >  1:1942533/1‑151 (MQ=255)
                             cGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGACGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCagg                                                                                             >  2:4984587/1‑151 (MQ=255)
                              gTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCagga                                                                                            <  1:1021588/151‑1 (MQ=255)
                                gaagaaGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCaggagg                                                                                          <  1:5099177/151‑1 (MQ=255)
                                 aagaagTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGc                                                                                         <  1:5149727/151‑1 (MQ=255)
                                    aagTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGa                                                                                      <  1:5420715/151‑1 (MQ=255)
                                     agTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAg                                                                                     >  1:3559732/1‑151 (MQ=255)
                                       tAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTg                                                                                   <  1:619216/151‑1 (MQ=255)
                                         aCCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCa                                                                                 <  1:952883/151‑1 (MQ=255)
                                          ccTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGCTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAg                                                                                <  1:4079945/151‑1 (MQ=255)
                                          ccTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAg                                                                                <  1:1112304/151‑1 (MQ=255)
                                            tGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAggg                                                                              <  1:4206801/151‑1 (MQ=255)
                                             gCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGc                                                                             <  1:757528/151‑1 (MQ=255)
                                              ccACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCt                                                                            <  2:4054371/151‑1 (MQ=255)
                                                aCCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTg                                                                          >  2:5342047/1‑151 (MQ=255)
                                                aCCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTg                                                                          >  1:3905956/1‑151 (MQ=255)
                                                     cggcggCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAg                                                                     >  1:2030559/1‑151 (MQ=255)
                                                      ggcggcTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACACAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAgg                                                                    >  1:3681190/1‑151 (MQ=255)
                                                       gcggcTCGACTTGAAGCCATGCTGAAGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGAAGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAggg                                                                   <  1:3694702/151‑1 (MQ=255)
                                                        cggcTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAg                                                                                                           >  1:1124645/1‑110 (MQ=255)
                                                        cggcTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAg                                                                                                           >  1:1125301/1‑110 (MQ=255)
                                                        cggcTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAg                                                                                                           <  2:1125301/110‑1 (MQ=255)
                                                         ggcTCGACTTGAAGCCCCGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAg                                                                 <  2:1864818/151‑1 (MQ=255)
                                                           cTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGt                                                               <  1:2556878/151‑1 (MQ=255)
                                                           cTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGt                                                               <  1:2896224/151‑1 (MQ=255)
                                                           cTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGt                                                               <  1:3793700/151‑1 (MQ=255)
                                                            tCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTg                                                              >  1:3292229/1‑151 (MQ=255)
                                                             cGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGa                                                             <  2:5331885/151‑1 (MQ=255)
                                                             cGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGa                                                             >  2:1636595/1‑151 (MQ=255)
                                                              gACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGAc                                                            <  1:4966107/151‑1 (MQ=255)
                                                               aCTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAggcgccg                                                           >  2:4092809/1‑146 (MQ=255)
                                                               aCTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACg                                                           >  2:4758038/1‑151 (MQ=255)
                                                               aCTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACg                                                           <  2:281635/151‑1 (MQ=255)
                                                                cTTGAAGCCACGCTGCGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACgg                                                          <  1:4033439/151‑1 (MQ=255)
                                                                  tGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGc                                                        <  1:2691560/151‑1 (MQ=255)
                                                                   gAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCa                                                       >  2:2838263/1‑151 (MQ=255)
                                                                    aaGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCAt                                                      <  1:4705310/151‑1 (MQ=255)
                                                                        cACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGAt                                                  >  1:3454556/1‑151 (MQ=255)
                                                                          cGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGTTTGCGAATTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGc                                                >  1:1814768/1‑151 (MQ=255)
                                                                           gCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGcc                                               >  2:3019652/1‑151 (MQ=255)
                                                                            cTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGccc                                              >  2:247064/1‑151 (MQ=255)
                                                                             tGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCg                                             <  1:3218101/151‑1 (MQ=255)
                                                                              gAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGt                                            <  2:3952338/151‑1 (MQ=255)
                                                                              gAGGCCAGCACGGCGGTCTGCGTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGAGAGTTCAGGCGGCGCAGGTTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGt                                            <  2:2633545/151‑1 (MQ=14)
                                                                               aGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTc                                           <  1:5342047/151‑1 (MQ=255)
                                                                               aGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTc                                           >  2:3996750/1‑151 (MQ=255)
                                                                               aGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATCCATGCCCGTc                                           <  2:1544430/151‑1 (MQ=255)
                                                                                ggCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCa                                          <  1:3768977/151‑1 (MQ=255)
                                                                                  ccaccccgtcgtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTGTTGCTTCTCAGCGCGCGCTAGTTCGGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAg                                        <  1:4086501/145‑1 (MQ=11)
                                                                                  ccAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAg                                        <  1:3601170/151‑1 (MQ=255)
                                                                                    aGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAgcg                                      <  1:1201751/151‑1 (MQ=255)
                                                                                    aGCACGTCGTCCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAgcg                                      <  1:5008321/151‑1 (MQ=255)
                                                                                     gCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAgcgc                                     >  1:4238876/1‑151 (MQ=255)
                                                                                     gCACGTCGTACTGCTTCTGCTTCTGCTTCTGATTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGCCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAgcgc                                     <  2:4259249/151‑1 (MQ=255)
                                                                                      cACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCt                                    >  2:4417489/1‑151 (MQ=255)
                                                                                       aCGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTg                                   >  1:2054707/1‑151 (MQ=255)
                                                                                        cgtcgtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTgg                                  <  1:1240874/151‑2 (MQ=255)
                                                                                        cgtcgtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTgg                                  <  2:5117976/151‑2 (MQ=255)
                                                                                           cgtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTgggcc                               >  2:4236483/1‑147 (MQ=255)
                                                                                            gtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGccc                             >  1:3106084/1‑151 (MQ=255)
                                                                                            gtTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGccc                             >  2:270408/1‑151 (MQ=16)
                                                                                             ttctgcttctgcttctgcttctgcttctgcttctgcTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGcccc                            <  1:852286/151‑1 (MQ=255)
                                                                                              tctgcttctgcttctgcttctgcttctgcttctgctTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCa                           <  2:5059353/151‑1 (MQ=16)
                                                                                               ttgcttcggattctgcttctgcttctgcttctgcttTTGCTTCTCAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑g                         <  1:2073309/150‑1 (MQ=11)
                                                                                                tgcttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGCCGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑gg                        >  2:2142805/1‑151 (MQ=255)
                                                                                                 gcttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGa                       >  1:1763211/1‑151 (MQ=16)
                                                                                                 gctgctgcttctgcttctgcttctgcttctgcttTTGCTTCTCAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGACGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGa                       <  2:3481242/151‑1 (MQ=11)
                                                                                                   ttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATg                     >  2:1052902/1‑151 (MQ=16)
                                                                                                   ttctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATg                     <  2:3049066/151‑1 (MQ=16)
                                                                                                     ctgcttctgcttctgcttctgcttctgcttTTGCTTCTAGGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATggg                   <  1:5295473/151‑1 (MQ=255)
                                                                                                     ctgcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATggg                   >  1:323683/1‑151 (MQ=255)
                                                                                                     ctgcttctgcttcggcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATggg                   <  1:3858899/151‑1 (MQ=255)
                                                                                                       gcttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGc                 <  2:1763211/151‑1 (MQ=16)
                                                                                                        cttctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGcc                >  2:317539/1‑151 (MQ=255)
                                                                                                         ttcggcttctgcttcggcttctgcttTTGCGTCAACGCGCGCGATGGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGCGCGCAGGGCTGGCGAAGGGGGGTGACGGGCTTTGATGCCCGTCACGCGCTGAGGACCCATGGATGGGGc                 <  2:2858868/151‑1 (MQ=11)
                                                                                                          tctgcttctgcttctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCtg              >  2:3273285/1‑150 (MQ=16)
                                                                                                               ttctgctgcggcttcggcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCtgcagcg         <  2:559177/151‑7 (MQ=1)
                                                                                                                      tctgcttctgcttTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCA‑GGATGGGGCCTGCAGCGCGGTCCt  >  1:399776/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                          
CGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCA‑GGATGGGGCCTTCGGCGCGGTCCT  >  CP009974/4030779‑4031048

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTCGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAG  >  CP009974/4030777‑4031056
                                                                                                                                                      |                                                                                                                                 
GTCGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGA                                                                                                                                   <  SRR10918230.213563/151‑1 (MQ=60)
GTCGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGA                                                                                                                                   >  SRR10918230.1110177/1‑151 (MQ=60)
GTCGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGA                                                                                                                                   <  SRR10918230.213769/151‑1 (MQ=60)
GTCGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGA                                                                                                                                   >  SRR10918230.4522296/1‑151 (MQ=60)
 TCGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGAT                                                                                                                                  <  SRR10918230.3118154/151‑1 (MQ=60)
   GATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAG                                                                                                                                <  SRR10918230.740774/151‑1 (MQ=60)
    ATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGT                                                                                                                               >  SRR10918230.3719514/1‑151 (MQ=60)
    ATTGCCTGCAAGCAGCTGTACCACGACCGTAAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGT                                                                                                                               >  SRR10918230.2304703/1‑151 (MQ=60)
     TTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTT                                                                                                                              <  SRR10918230.4491606/151‑1 (MQ=60)
     TTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTT                                                                                                                              <  SRR10918230.4101531/151‑1 (MQ=60)
     TTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTT                                                                                                                              <  SRR10918230.4094819/151‑1 (MQ=60)
       GCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTCCTAAGCGCGCGATAGTCCA                                                                                                                            >  SRR10918230.4129383/1‑151 (MQ=60)
        CCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAG                                                                                                                           <  SRR10918230.4311791/151‑1 (MQ=60)
             AAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGAC                                                                                                                      <  SRR10918230.3165554/151‑1 (MQ=60)
             AAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGAC                                                                                                                      <  SRR10918230.138971/151‑1 (MQ=60)
              AGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACG                                                                                                                     <  SRR10918230.3873369/151‑1 (MQ=60)
               GCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGC                                                                                                                    <  SRR10918230.2056730/151‑1 (MQ=60)
                   attgggtgagtttctgtcccgtcctatggataaggatatcacctgcatggctgccctcggcctagtgtccctggaagaacgctcttgctctttctctttctttggatttTGCTTTTGCTTCTAAGCGCGCTGTAGTTCAGGCGATGCAGAT                                                                                                                <  SRR10918230.1473609/42‑1 (MQ=60)
                   CTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGAT                                                                                                                <  SRR10918230.2303894/151‑1 (MQ=60)
                    TGTACCACGACCGTGAAGAAGTAACGTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATT                                                                                                               <  SRR10918230.4024666/151‑1 (MQ=60)
                     GTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTG                                                                                                              <  SRR10918230.1884966/151‑1 (MQ=60)
                       ACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCG                                                                                                            >  SRR10918230.1993489/1‑151 (MQ=60)
                       ACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCG                                                                                                            >  SRR10918230.2562791/1‑151 (MQ=60)
                              CCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAG                                                                                                     >  SRR10918230.1950756/1‑151 (MQ=60)
                              CCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAG                                                                                                     >  SRR10918230.543714/1‑151 (MQ=60)
                                GTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGA                                                                                                   <  SRR10918230.1025807/151‑1 (MQ=60)
                                  GAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGG                                                                                                 <  SRR10918230.4512365/151‑1 (MQ=60)
                                  GAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGG                                                                                                 <  SRR10918230.5122565/151‑1 (MQ=60)
                                   catctggataaggatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGC                                                                                                >  SRR10918230.4723140/94‑151 (MQ=12)
                                   catctggataaggatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGC                                                                                                >  SRR10918230.634030/94‑151 (MQ=12)
                                   AAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGC                                                                                                <  SRR10918230.5173402/151‑1 (MQ=60)
                                   AAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGC                                                                                                <  SRR10918230.1452161/151‑1 (MQ=60)
                                    atctggataaggatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCC                                                                                               <  SRR10918230.1438992/59‑1 (MQ=12)
                                    atctggataaggatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCC                                                                                               >  SRR10918230.3048030/93‑151 (MQ=12)
                                     GAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCG                                                                                              >  SRR10918230.4403400/1‑151 (MQ=60)
                                      ggtcatctggataaggatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGA                                                                                             >  SRR10918230.4170952/91‑151 (MQ=12)
                                      ggtcatctggataaggagagcgccggcatggctgacctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGA                                                                                             <  SRR10918230.1278604/61‑1 (MQ=12)
                                      AAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGA                                                                                             >  SRR10918230.2169158/1‑151 (MQ=60)
                                      AAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGA                                                                                             <  SRR10918230.5446268/151‑1 (MQ=60)
                                       tggataaggatatcgtctgcatggctgccctcggtctggtgtacctggaagaatgctgttgcgcttgctcttgctcttgcttttgatttTGCTTTTGCTTCTAAGCGCGAGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAG                                                                                            <  SRR10918230.1050318/62‑1 (MQ=11)
                                       AGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAG                                                                                            >  SRR10918230.3575127/1‑151 (MQ=60)
                                         TAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTG                                                                                          <  SRR10918230.621428/151‑1 (MQ=60)
                                          atctggataaggatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGC                                                                                         >  SRR10918230.185645/87‑151 (MQ=3)
                                          atctggataaggatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGC                                                                                         <  SRR10918230.4232061/65‑1 (MQ=3)
                                          ataaggatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGC                                                                                         <  SRR10918230.4243535/65‑1 (MQ=3)
                                           taaggatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCA                                                                                        >  SRR10918230.959201/86‑151 (MQ=0)
                                           ACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCA                                                                                        <  SRR10918230.956719/151‑1 (MQ=60)
                                           ACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCA                                                                                        <  SRR10918230.2832127/151‑1 (MQ=60)
                                            aaggatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAG                                                                                       <  SRR10918230.5283276/67‑1 (MQ=0)
                                            CCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGCTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAG                                                                                       <  SRR10918230.4098132/151‑1 (MQ=60)
                                            CCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAG                                                                                       <  SRR10918230.1120739/151‑1 (MQ=60)
                                            CCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAG                                                                                       <  SRR10918230.1117037/151‑1 (MQ=60)
                                              ggataaggatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGG                                                                                     >  SRR10918230.4548987/83‑151 (MQ=0)
                                              TGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGG                                                                                     <  SRR10918230.4225640/151‑1 (MQ=60)
                                               gatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGC                                                                                    >  SRR10918230.4321681/82‑151 (MQ=0)
                                               GCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGC                                                                                    <  SRR10918230.760376/151‑1 (MQ=60)
                                                CCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCT                                                                                   >  SRR10918230.865287/1‑151 (MQ=60)
                                                  atcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTG                                                                                 >  SRR10918230.833395/79‑151 (MQ=0)
                                                  aaggatatcgccggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTG                                                                                 >  SRR10918230.1880996/79‑151 (MQ=0)
                                                  ACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTG                                                                                 >  SRR10918230.3923049/1‑151 (MQ=60)
                                                    CTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCG                                                                               <  SRR10918230.2617364/151‑1 (MQ=60)
                                                      GCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAA                                                                             <  SRR10918230.4022173/151‑1 (MQ=60)
                                                       cggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAG                                                                            >  SRR10918230.2619978/74‑151 (MQ=0)
                                                       cggcatggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAG                                                                            >  SRR10918230.2620171/74‑151 (MQ=0)
                                                       CGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAG                                                                            >  SRR10918230.2041822/1‑151 (MQ=60)
                                                       CGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAG                                                                            >  SRR10918230.2039137/1‑151 (MQ=60)
                                                        ggcgtggctgcccttgcgctggtgtacctggaagggcgctcttgctcttgcttttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGG                                                                           <  SRR10918230.920542/79‑1 (MQ=0)
                                                        GGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACACAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGG                                                                           >  SRR10918230.3697166/1‑151 (MQ=60)
                                                         GCGGCTCGACTTGAAGCCATGCTGAAGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGAAGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGG                                                                          <  SRR10918230.3710747/151‑1 (MQ=60)
                                                          CGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGAgatcggaagagcacacgtctgaactccagtcacggacatc                                                                         >  SRR10918230.1129441/1‑111 (MQ=60)
                                                          CGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGAgatcggaagagcacacgtctgaactccagtcacggacatc                                                                         >  SRR10918230.1130104/1‑111 (MQ=60)
                                                             ggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGT                                                                      <  SRR10918230.771139/84‑1 (MQ=0)
                                                             ggctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGG                                                                      >  SRR10918230.2957403/68‑151 (MQ=0)
                                                             CTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGT                                                                      <  SRR10918230.2908830/151‑1 (MQ=60)
                                                             CTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGT                                                                      <  SRR10918230.3810204/151‑1 (MQ=60)
                                                             CTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGT                                                                      <  SRR10918230.2567947/151‑1 (MQ=60)
                                                             CTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGT                                                                      >  SRR10918230.751375/1‑151 (MQ=60)
                                                              TCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTG                                                                     >  SRR10918230.3306534/1‑151 (MQ=60)
                                                               CGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGA                                                                    <  SRR10918230.2430324/151‑1 (MQ=60)
                                                               CGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGA                                                                    >  SRR10918230.3580030/1‑151 (MQ=60)
                                                                tgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGCGGTGAC                                                                   >  SRR10918230.547919/65‑151 (MQ=0)
                                                                GACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGAC                                                                   <  SRR10918230.4988741/151‑1 (MQ=60)
                                                                 tccctcggcctggtgtacctggaagtacgctcttgctcttgctcttgctcttgcttttgatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTGAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACG                                                                  <  SRR10918230.1315243/88‑1 (MQ=0)
                                                                  ccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGGGACGG                                                                 >  SRR10918230.1261572/63‑151 (MQ=0)
                                                                  CTTGAAGCCACGCTGCGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGG                                                                 <  SRR10918230.4051329/151‑1 (MQ=60)
                                                                    gctgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGC                                                               >  SRR10918230.1448468/61‑151 (MQ=0)
                                                                    TGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGC                                                               <  SRR10918230.32802/151‑1 (MQ=60)
                                                                    TGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGC                                                               <  SRR10918230.2703338/151‑1 (MQ=60)
                                                                      tgccctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCAT                                                             >  SRR10918230.4407980/59‑151 (MQ=0)
                                                                      AAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCAT                                                             <  SRR10918230.4726654/151‑1 (MQ=60)
                                                                       ggactggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATG                                                            <  SRR10918230.5411702/94‑1 (MQ=0)
                                                                          ctggtgtacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGAT                                                         >  SRR10918230.27541/55‑151 (MQ=0)
                                                                          ctcggcctggtgtacctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGAT                                                         >  SRR10918230.657199/55‑151 (MQ=0)
                                                                          CACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGAT                                                         >  SRR10918230.3469461/1‑151 (MQ=60)
                                                                            CGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGTTTGCGAATTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGC                                                       >  SRR10918230.1822575/1‑151 (MQ=60)
                                                                              gcctggtgtacctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCC                                                     >  SRR10918230.765242/51‑151 (MQ=0)
                                                                               cctggtgtacctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGagatcggaagagca                                                    >  SRR10918230.4763418/50‑137 (MQ=0)
                                                                               TGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCG                                                    <  SRR10918230.3232090/151‑1 (MQ=60)
                                                                                tacctggaagaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGT                                                   >  SRR10918230.1361738/49‑151 (MQ=0)
                                                                                 AGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTC                                                  <  SRR10918230.5367041/151‑1 (MQ=60)
                                                                                  GGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCA                                                 <  SRR10918230.3785362/151‑1 (MQ=59)
                                                                                  GGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCA                                                 >  SRR10918230.4090531/1‑151 (MQ=59)
                                                                                    CCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAG                                               <  SRR10918230.3616764/151‑1 (MQ=54)
                                                                                    CCACCCCGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTGTTGCTTCTCAGCGCGCGCTAGTTCGGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAG                                               <  SRR10918230.4104714/151‑1 (MQ=27)
                                                                                      tacctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCG                                             >  SRR10918230.3848127/43‑151 (MQ=0)
                                                                                      AGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCG                                             >  SRR10918230.5003167/1‑151 (MQ=50)
                                                                                      AGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCG                                             <  SRR10918230.1206968/151‑1 (MQ=50)
                                                                                      AGCACGTCGTCCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCG                                             <  SRR10918230.5031238/151‑1 (MQ=43)
                                                                                       GCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGC                                            >  SRR10918230.4257888/1‑151 (MQ=48)
                                                                                        cctggaagaacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCT                                           <  SRR10918230.3862921/111‑1 (MQ=0)
                                                                                        agaacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGGTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCT                                           >  SRR10918230.2454900/41‑151 (MQ=0)
                                                                                         ACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG                                          >  SRR10918230.2063376/1‑151 (MQ=44)
                                                                                         ACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG                                          >  SRR10918230.2070378/1‑151 (MQ=44)
                                                                                         ACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG                                          >  SRR10918230.5266875/1‑151 (MQ=44)
                                                                                          aacgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGG                                         >  SRR10918230.3660964/39‑151 (MQ=0)
                                                                                          CGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGG                                         <  SRR10918230.1246326/151‑1 (MQ=40)
                                                                                           acgctcttgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGG                                        >  SRR10918230.640075/38‑151 (MQ=0)
                                                                                              GTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCC                                    >  SRR10918230.3697619/1‑151 (MQ=22)
                                                                                              GTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCC                                    >  SRR10918230.3119618/1‑151 (MQ=22)
                                                                                               TTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCC                                   <  SRR10918230.855610/151‑1 (MQ=18)
                                                                                               TTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCC                                   <  SRR10918230.699962/151‑1 (MQ=18)
                                                                                                aacgctcttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATCCATGCCCGTCctctccctccgcaag                                   >  SRR10918230.3225841/33‑136 (MQ=0)
                                                                                                 ttgcttcggaTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTCAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAG                                 <  SRR10918230.2082066/141‑1 (MQ=16)
                                                                                                   GCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGA                               >  SRR10918230.1770844/1‑151 (MQ=16)
                                                                                                     TTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATG                             <  SRR10918230.5148601/151‑1 (MQ=16)
                                                                                                      cttgctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATCCATGCCCGTCAtctccagatcggaagagcac                             >  SRR10918230.4976042/27‑131 (MQ=0)
                                                                                                       CTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAGGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGG                           <  SRR10918230.5320124/151‑1 (MQ=16)
                                                                                                       CTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGG                           >  SRR10918230.324772/1‑151 (MQ=16)
                                                                                                       CTGCTTCTGCTTCGGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGG                           <  SRR10918230.3875742/151‑1 (MQ=16)
                                                                                                         gctcttgctcttgctttggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGAGCAGGGCTTGCGGAGGGAGGTGACGGGCATctattcccctcatctccctccgcaaagccgggccg                          >  SRR10918230.497905/24‑116 (MQ=11)
                                                                                                           TTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGCCT                       >  SRR10918230.1217039/1‑151 (MQ=16)
                                                                                                            gactggagttcagacgtgtgctcttccgatctTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAcctgcctccgcaaggcctgcgttcgg                       <  SRR10918230.556074/119‑27 (MQ=25)
                                                                                                            TCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGCCTG                      >  SRR10918230.1187106/1‑151 (MQ=11)
                                                                                                              TGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGCGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGCCTgca                    <  SRR10918230.4111624/151‑4 (MQ=7)
                                                                                                                CTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGCCTgcagc                  >  SRR10918230.3352682/1‑146 (MQ=3)
                                                                                                                 TTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTG‑GGCCCCAGGATGGGGCCTGCAGCG                 <  SRR10918230.3956859/151‑1 (MQ=1)
                                                                                                                   CTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATCCATGCCCGTCAcctgcctccgcaaggcctgcgttcggcctccag                >  SRR10918230.740119/1‑118 (MQ=52)
                                                                                                                          ggatttTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCagatcggaagagcacacgtctgaactccagtcacgga         >  SRR10918230.4258695/7‑114 (MQ=0)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCTGTAGTTCAGGCGATGCAGATTGCTACT                                                                                                         <  SRR10918230.1648525/49‑1 (MQ=60)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTTAGGCGATG                                                                                                                     >  SRR10918230.4567052/1‑37 (MQ=23)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGN                                                                                                                 <  SRR10918230.5214858/41‑1 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCGA                                                                                             >  SRR10918230.1434277/1‑61 (MQ=12)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGAGGCCG                                                                                              <  SRR10918230.2638085/60‑1 (MQ=12)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGGA                                                                                                   >  SRR10918230.78795/1‑55 (MQ=13)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTCAGG                                                                                                    >  SRR10918230.727345/1‑54 (MQ=13)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACTTC                                                                                                       <  SRR10918230.1932859/51‑1 (MQ=13)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGACT                                                                                                         >  SRR10918230.1102147/1‑49 (MQ=25)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGAC                                                                                                          >  SRR10918230.3954103/1‑48 (MQ=25)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGAC                                                                                                          >  SRR10918230.1291511/1‑48 (MQ=25)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGA                                                                                                           >  SRR10918230.1203864/1‑47 (MQ=25)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATTGCGA                                                                                                           >  SRR10918230.4606780/1‑47 (MQ=25)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATT                                                                                                               >  SRR10918230.3170194/1‑43 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGATT                                                                                                               <  SRR10918230.1458809/43‑1 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGAT                                                                                                                >  SRR10918230.4707171/1‑42 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGAT                                                                                                                >  SRR10918230.3859113/1‑42 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGA                                                                                                                 >  SRR10918230.4520820/1‑41 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGA                                                                                                                 <  SRR10918230.3428201/41‑1 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGA                                                                                                                 >  SRR10918230.2750084/1‑41 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAGA                                                                                                                 >  SRR10918230.763189/1‑41 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAG                                                                                                                  <  SRR10918230.202209/40‑1 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAG                                                                                                                  >  SRR10918230.4738542/1‑40 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAG                                                                                                                  >  SRR10918230.3060009/1‑40 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCAG                                                                                                                  >  SRR10918230.844068/1‑40 (MQ=24)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGCA                                                                                                                   >  SRR10918230.5185140/1‑39 (MQ=25)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGC                                                                                                                    >  SRR10918230.3377003/1‑38 (MQ=25)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATGC                                                                                                                    <  SRR10918230.5132275/38‑1 (MQ=25)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATG                                                                                                                     >  SRR10918230.1029384/1‑37 (MQ=26)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGATG                                                                                                                     >  SRR10918230.667352/1‑37 (MQ=26)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGAT                                                                                                                      <  SRR10918230.2986707/36‑1 (MQ=26)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGAT                                                                                                                      >  SRR10918230.3066113/1‑36 (MQ=26)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGA                                                                                                                       <  SRR10918230.3881379/35‑1 (MQ=27)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCG                                                                                                                        >  SRR10918230.501422/1‑34 (MQ=27)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGC                                                                                                                         >  SRR10918230.2759134/1‑33 (MQ=27)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGG                                                                                                                          <  SRR10918230.1798517/32‑1 (MQ=27)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAG                                                                                                                           >  SRR10918230.153003/1‑31 (MQ=27)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGTAGTTCAG                                                                                                                           <  SRR10918230.902657/31‑1 (MQ=27)
                                                                                                                                TGCTTTTGCTTCTAAGCGCGCGGCAGTTCAGGCGATGCAGATTGCGACTT                                                                                                        >  SRR10918230.3198108/1‑50 (MQ=60)
                                                                                                                                 GCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGA                                                                                                   >  SRR10918230.1194678/1‑54 (MQ=11)
                                                                                                                                 GCTTTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATCCATGCCCGTCAcctgcctccgcaaggcctgcgttcggcctccagaagtcgcgaagatc  >  SRR10918230.3166402/1‑104 (MQ=22)
                                                                                                                                                      |                                                                                                                                 
GTCGATTGCCTGCAAGCAGCTGTACCACGACCGTGAAGAAGTAACCTGCCACCTGCGGCGGCTCGACTTGAAGCCACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTCTGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAG  >  CP009974/4030777‑4031056

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