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 2,738,529 G→T A248S (GCC→TCT)  RPPX_12165 → benzene 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099742,738,5290GT98.8% 395.4 / ‑4.3 160A248S (GCC→TCT) RPPX_12165benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (78/80);  total (78/82)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

AGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCCGCGCAGACTG  >  CP009974/2738393‑2738674
                                                                                                                                        |                                                                                                                                                 
 ggtaaCTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAgg                                                                                                                                    >  2:324255/4‑151 (MQ=0)
 ggtaaCTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAgg                                                                                                                                    >  2:4119968/4‑151 (MQ=0)
          aaGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGtt                                                                                                                           >  2:2054432/1‑151 (MQ=0)
                      gagaaCGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCg                                                                                                               >  2:2858554/4‑151 (MQ=0)
                           cGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCAt                                                                                                          <  1:4204967/151‑1 (MQ=2)
                               gcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGac                                                                                                      <  1:5318288/151‑1 (MQ=2)
                                   aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                                                                                  >  2:668703/1‑151 (MQ=2)
                                   aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                                                                                  >  2:2644653/1‑151 (MQ=2)
                                   aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                                                                                  >  2:2646067/1‑151 (MQ=2)
                                   aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                                                                                  <  2:4236934/151‑1 (MQ=11)
                                   aCGGCTATCATGTTAGCGCGGTGCACTGGAACTACGCCGCTAATCGCAACCAGCGCAAGCAGCGCGAGGCAGGCGACTAAGTTAAATCCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTTGTTGTTTCTATTCCTTCGAGCATGGACTCct                                                                                                  <  2:3041600/151‑1 (MQ=2)
                                    cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg                                                                                                 >  1:2506923/1‑151 (MQ=11)
                                    cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg                                                                                                 <  2:2365483/151‑1 (MQ=2)
                                    cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGGTTCTATTCCTTCCACCATGGACACctg                                                                                                 >  2:4326824/1‑151 (MQ=2)
                                     ggCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctgc                                                                                                >  2:493223/1‑151 (MQ=11)
                                      gCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctgct                                                                                               >  1:4515264/1‑151 (MQ=11)
                                       cTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTCCTATTCCTTCGACCATGGACACctgctg                                                                                              >  2:3228592/1‑151 (MQ=11)
                                          tcatgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctc                                                                                           <  2:4031951/144‑1 (MQ=11)
                                          tcatgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctc                                                                                           <  1:4858875/144‑1 (MQ=11)
                                          tcatgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctc                                                                                           <  1:3610142/144‑1 (MQ=2)
                                            atgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTg                                                                                         >  2:4730553/6‑151 (MQ=11)
                                            atgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTg                                                                                         >  1:220555/6‑151 (MQ=11)
                                             tgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTgg                                                                                        <  1:3256087/147‑1 (MQ=11)
                                               ttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGAc                                                                                      <  2:2715610/149‑1 (MQ=11)
                                                taGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACt                                                                                     <  2:1772422/150‑1 (MQ=14)
                                                   cTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTATGCAGGGCGGTGGTTTCTATTCCTTCGACCATGGACACCTCCTGCTCTGGCCTCgg                                                                                  >  2:3366426/1‑151 (MQ=14)
                                                   cTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCgg                                                                                  >  2:1838182/1‑151 (MQ=14)
                                                     ccGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTg                                                                                <  1:2099453/151‑1 (MQ=14)
                                                     ccGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTg                                                                                <  2:1923450/151‑1 (MQ=14)
                                                      cGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTgg                                                                               >  1:4538964/1‑151 (MQ=14)
                                                      cGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTgg                                                                               <  1:104673/151‑1 (MQ=14)
                                                      cGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTgg                                                                               <  1:105033/151‑1 (MQ=14)
                                                       gTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTggg                                                                              <  1:809449/151‑1 (MQ=14)
                                                           aCTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGaa                                                                          <  2:5313792/151‑1 (MQ=14)
                                                            cTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTTCTCTGGACTCGGTGGGCGAAc                                                                         >  1:1760739/1‑151 (MQ=17)
                                                             tGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcc                                                                        >  2:321968/1‑151 (MQ=14)
                                                               gAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcccc                                                                      <  2:1385394/151‑1 (MQ=17)
                                                                  cTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAg                                                                   <  2:219614/151‑1 (MQ=17)
                                                                   tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                                  <  2:1714847/151‑1 (MQ=17)
                                                                   tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                                  >  2:3905787/1‑151 (MQ=17)
                                                                    aCGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAGCCCCGAGGa                                                                 >  2:1298712/1‑151 (MQ=17)
                                                                     cgccgcCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAt                                                                <  1:4949585/151‑1 (MQ=17)
                                                                         gCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGcc                                                            <  1:890926/151‑1 (MQ=17)
                                                                          ccACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGccc                                                           >  1:786268/1‑151 (MQ=17)
                                                                             ctcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGc                                                        >  1:3522604/3‑151 (MQ=18)
                                                                             ctcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGc                                                        >  1:1429328/3‑151 (MQ=18)
                                                                              tcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCt                                                       <  2:368254/150‑1 (MQ=18)
                                                                                aGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTa                                                     <  1:4031362/151‑1 (MQ=18)
                                                                                 gAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTAt                                                    <  2:3773630/151‑1 (MQ=18)
                                                                                  aaCCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATg                                                   >  1:425207/1‑151 (MQ=18)
                                                                                    ccAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAg                                                 >  1:4701486/1‑151 (MQ=18)
                                                                                     cAGCGCCAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGGGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAgc                                                <  2:3863099/151‑1 (MQ=255)
                                                                                      aGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAgcg                                               <  2:2573089/151‑1 (MQ=18)
                                                                                      aGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAgcg                                               <  2:2401391/151‑1 (MQ=18)
                                                                                        cgcAAGCAGCGCGAGGCAGGCGACGAAATCAACACCAGGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                                             <  2:2752582/151‑1 (MQ=255)
                                                                                        cgcAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                                             <  2:834715/151‑1 (MQ=255)
                                                                                        cgcAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                                             <  2:4735549/151‑1 (MQ=255)
                                                                                        cgcAAGCAGCGCGAGCCAGGCGACGACATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAGCGCCCAGCTTATGAGCGcc                                             <  2:3538825/151‑1 (MQ=255)
                                                                                         gcAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCg                                            >  1:4393400/1‑151 (MQ=18)
                                                                                          cAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGt                                           >  2:4508740/1‑151 (MQ=255)
                                                                                          cAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGt                                           >  2:4510787/1‑151 (MQ=255)
                                                                                          cAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGt                                           >  1:4221109/1‑151 (MQ=255)
                                                                                           aaGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTg                                          >  1:2552307/1‑151 (MQ=255)
                                                                                           aaGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTg                                          >  1:772281/1‑151 (MQ=255)
                                                                                           aaGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTg                                          <  1:592820/151‑1 (MQ=255)
                                                                                              cagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAcc                                       <  1:5143558/151‑1 (MQ=18)
                                                                                                gcgcGAGGCAGGCGACGCAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCCCTGGACTCGGTGGGCGACCTCCGAGGATCGCCCCGCTTATTAGCGCCGTGACCac                                     >  2:2698408/1‑150 (MQ=255)
                                                                                                gcgcGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCaa                                     <  1:22896/151‑1 (MQ=255)
                                                                                                gcgcGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCaa                                     >  1:5195923/1‑151 (MQ=255)
                                                                                                gcgcGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCaa                                     >  1:5196064/1‑151 (MQ=255)
                                                                                                 cgcgAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAAt                                    >  2:2961693/1‑151 (MQ=18)
                                                                                                  gcgAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAAtt                                   >  2:1685961/1‑151 (MQ=255)
                                                                                                   cgAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTg                                  >  2:2850142/1‑151 (MQ=255)
                                                                                                     aggcaggcGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGc                                <  1:144461/151‑1 (MQ=255)
                                                                                                      ggcgggcgACGAAATCAAAACCAGGAGCGCCGGCTCGTGGGCTAAGCAGGGTGGTGGGTTCTAGTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCt                               <  2:5143999/151‑1 (MQ=255)
                                                                                                      ggcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCt                               >  1:3288601/8‑151 (MQ=255)
                                                                                                      ggcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCt                               <  1:3818708/144‑1 (MQ=21)
                                                                                                       gcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTc                              >  1:1306264/7‑151 (MQ=255)
                                                                                                       gcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTc                              <  2:200642/145‑1 (MQ=255)
                                                                                                       gcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTc                              >  2:3271331/7‑151 (MQ=255)
                                                                                                        caggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCa                             <  1:727060/146‑1 (MQ=255)
                                                                                                           gcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGa                          >  1:4142511/3‑151 (MQ=255)
                                                                                                             gACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAtt                        <  2:1133734/151‑1 (MQ=255)
                                                                                                             gACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAtt                        <  1:1408691/151‑1 (MQ=255)
                                                                                                              aCGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAttt                       >  1:4149203/1‑151 (MQ=255)
                                                                                                                gaaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                                             >  1:4462052/4‑127 (MQ=21)
                                                                                                                gaaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                                             <  2:4462052/124‑1 (MQ=255)
                                                                                                                gaaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                                             >  1:4462108/4‑127 (MQ=21)
                                                                                                                gaaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCg                     <  1:4362098/148‑1 (MQ=255)
                                                                                                                gaaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCg                     >  1:2637477/4‑151 (MQ=255)
                                                                                                                  aaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGc                   <  2:3967119/150‑1 (MQ=25)
                                                                                                                  aaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGc                   >  2:5044641/2‑151 (MQ=255)
                                                                                                                   aTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGcc                  >  1:202987/1‑151 (MQ=255)
                                                                                                                     cAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAg                <  1:1232330/151‑1 (MQ=255)
                                                                                                                      aaaaCCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAgg               <  1:2927645/148‑1 (MQ=25)
                                                                                                                      aaaaCCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAgg               <  1:2925308/148‑1 (MQ=25)
                                                                                                                          ccATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCcg           >  1:3870207/1‑151 (MQ=255)
                                                                                                                           cATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCcgc          <  2:1804417/151‑1 (MQ=255)
                                                                                                                           cATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCcgc          <  2:1935188/151‑1 (MQ=255)
                                                                                                                           cATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCcgc          <  1:3434643/151‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                 
AGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCCGCGCAGACTG  >  CP009974/2738393‑2738674

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCCGC  >  CP009974/2738379‑2738666
                                                                                                                                                      |                                                                                                                                         
TCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCG                                                                                                                                           >  SRR10918230.3891025/1‑151 (MQ=0)
TCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCG                                                                                                                                           <  SRR10918230.2816143/151‑1 (MQ=0)
TCGTACATCTACGAAGGCAACTAGAAGCTCACCGCCCAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCG                                                                                                                                           <  SRR10918230.1060081/151‑1 (MQ=0)
     CATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGG                                                                                                                                      >  SRR10918230.4252741/1‑151 (MQ=0)
             AAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCA                                                                                                                              <  SRR10918230.1375586/151‑1 (MQ=0)
             AAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCA                                                                                                                              >  SRR10918230.502404/1‑151 (MQ=0)
             AAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCA                                                                                                                              <  SRR10918230.4805898/151‑1 (MQ=0)
              AGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAG                                                                                                                             <  SRR10918230.1577503/151‑1 (MQ=0)
              AGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAG                                                                                                                             <  SRR10918230.325859/151‑1 (MQ=0)
                GTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGG                                                                                                                           <  SRR10918230.789566/151‑1 (MQ=0)
                GTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGG                                                                                                                           >  SRR10918230.4623686/1‑151 (MQ=0)
                     TGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCG                                                                                                                      >  SRR10918230.4469737/1‑151 (MQ=0)
                           CTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCT                                                                                                                <  SRR10918230.388504/151‑1 (MQ=0)
                            TGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTA                                                                                                               >  SRR10918230.4848879/1‑151 (MQ=0)
                             GACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTAT                                                                                                              >  SRR10918230.1167569/1‑151 (MQ=0)
                               CCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTC                                                                                                            >  SRR10918230.201305/1‑151 (MQ=0)
                                CGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCC                                                                                                           >  SRR10918230.2250580/1‑151 (MQ=0)
                                 GCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGCGGTTTCTATTCCT                                                                                                          >  SRR10918230.2695884/1‑151 (MQ=0)
                                 GCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGCGGATTCTATTCCT                                                                                                          >  SRR10918230.2693964/1‑151 (MQ=0)
                                 GCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCT                                                                                                          >  SRR10918230.5040984/1‑151 (MQ=0)
                                  CCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTT                                                                                                         <  SRR10918230.5002979/151‑1 (MQ=0)
                                   CGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTC                                                                                                        <  SRR10918230.2166137/151‑1 (MQ=0)
                                     AAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGA                                                                                                      >  SRR10918230.20154/1‑151 (MQ=0)
                                      GAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGAC                                                                                                     <  SRR10918230.4786773/151‑1 (MQ=0)
                                        ACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCA                                                                                                   <  SRR10918230.177090/151‑1 (MQ=0)
                                         CGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCAT                                                                                                  >  SRR10918230.2654119/1‑151 (MQ=0)
                                         CGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCAT                                                                                                  <  SRR10918230.4223794/151‑1 (MQ=0)
                                          GGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATG                                                                                                 <  SRR10918230.3511779/151‑1 (MQ=0)
                                           GCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGG                                                                                                >  SRR10918230.5250218/1‑151 (MQ=1)
                                             GCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGAC                                                                                              >  SRR10918230.5039585/1‑151 (MQ=4)
                                             GCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGAC                                                                                              <  SRR10918230.5343107/151‑1 (MQ=4)
                                              CCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACA                                                                                             <  SRR10918230.746282/151‑1 (MQ=5)
                                                 ACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCT                                                                                          >  SRR10918230.2607789/1‑151 (MQ=12)
                                                  CGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTCAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTG                                                                                         >  SRR10918230.2517728/1‑151 (MQ=14)
                                                   GGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGC                                                                                        <  SRR10918230.708753/151‑1 (MQ=14)
                                                    GCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCT                                                                                       <  SRR10918230.1749834/151‑1 (MQ=14)
                                                    GCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCT                                                                                       >  SRR10918230.4535676/1‑151 (MQ=14)
                                                     ctatcatgttAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTG                                                                                      <  SRR10918230.238250/141‑1 (MQ=15)
                                                     ctatcatgttAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTG                                                                                      <  SRR10918230.2367765/141‑1 (MQ=15)
                                                     ctatcatgttAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTG                                                                                      <  SRR10918230.3865058/141‑1 (MQ=15)
                                                       atcatgttAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCT                                                                                    <  SRR10918230.2037743/143‑1 (MQ=15)
                                                        tcatgttAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTC                                                                                   <  SRR10918230.3625773/144‑1 (MQ=17)
                                                        tcatgttAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTC                                                                                   <  SRR10918230.4880980/144‑1 (MQ=17)
                                                         CATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCT                                                                                  <  SRR10918230.1421798/151‑1 (MQ=17)
                                                          atgttAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTG                                                                                 >  SRR10918230.221300/6‑151 (MQ=20)
                                                           tgttAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGG                                                                                <  SRR10918230.3270233/147‑1 (MQ=21)
                                                                   CCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTG                                                                        <  SRR10918230.2108329/151‑1 (MQ=32)
                                                                    CGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGG                                                                       <  SRR10918230.105386/151‑1 (MQ=33)
                                                                    CGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGG                                                                       >  SRR10918230.4559480/1‑151 (MQ=33)
                                                                    CGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGG                                                                       <  SRR10918230.105024/151‑1 (MQ=33)
                                                                     GTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGG                                                                      <  SRR10918230.812567/151‑1 (MQ=35)
                                                                         ACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAA                                                                  >  SRR10918230.4758246/1‑151 (MQ=41)
                                                                         ACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAA                                                                  >  SRR10918230.2173658/1‑151 (MQ=41)
                                                                          CTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTTCTCTGGACTCGGTGGGCGAAC                                                                 >  SRR10918230.1768360/1‑151 (MQ=39)
                                                                            GGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCC                                                               <  SRR10918230.796879/151‑1 (MQ=45)
                                                                                 TACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGG                                                          <  SRR10918230.4536127/151‑1 (MQ=53)
                                                                                 TACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGG                                                          <  SRR10918230.968426/151‑1 (MQ=53)
                                                                                   CGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAT                                                        <  SRR10918230.4972108/151‑1 (MQ=56)
                                                                                      CGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTTCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGtttgc                                                     >  SRR10918230.409737/1‑146 (MQ=52)
                                                                                       GCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCC                                                    <  SRR10918230.894451/151‑1 (MQ=60)
                                                                                        CCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCC                                                   >  SRR10918230.516597/1‑151 (MQ=60)
                                                                                        CCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCC                                                   >  SRR10918230.789268/1‑151 (MQ=60)
                                                                                           CTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGC                                                >  SRR10918230.3537817/1‑151 (MQ=60)
                                                                                           CTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGC                                                >  SRR10918230.1435715/1‑151 (MQ=60)
                                                                                              AGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTA                                             <  SRR10918230.4049241/151‑1 (MQ=60)
                                                                                               GAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTAT                                            <  SRR10918230.3016939/151‑1 (MQ=60)
                                                                                                AACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATG                                           >  SRR10918230.426677/1‑151 (MQ=60)
                                                                                                  CCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAG                                         >  SRR10918230.4722814/1‑151 (MQ=60)
                                                                                                       GCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCG                                    >  SRR10918230.4413178/1‑151 (MQ=60)
                                                                                                        CAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGT                                   >  SRR10918230.4240038/1‑151 (MQ=60)
                                                                                                         AAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTG                                  >  SRR10918230.775209/1‑151 (MQ=60)
                                                                                                         AAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTG                                  >  SRR10918230.2563356/1‑151 (MQ=60)
                                                                                                         AAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTG                                  <  SRR10918230.594933/151‑1 (MQ=60)
                                                                                                          AGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGA                                 >  SRR10918230.4335998/1‑151 (MQ=60)
                                                                                                           GCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAC                                <  SRR10918230.20654/151‑1 (MQ=60)
                                                                                                            CAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACC                               <  SRR10918230.4515279/151‑1 (MQ=60)
                                                                                                            CAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACC                               <  SRR10918230.5167199/151‑1 (MQ=60)
                                                                                                              GCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAA                             <  SRR10918230.22963/151‑1 (MQ=60)
                                                                                                              GCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAA                             >  SRR10918230.5219878/1‑151 (MQ=60)
                                                                                                              GCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAA                             >  SRR10918230.5220021/1‑151 (MQ=60)
                                                                                                                   AGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGC                        <  SRR10918230.144823/151‑1 (MQ=60)
                                                                                                                   AGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGC                        <  SRR10918230.144965/151‑1 (MQ=60)
                                                                                                                    GGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCT                       >  SRR10918230.3302890/1‑151 (MQ=60)
                                                                                                                    GGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCT                       <  SRR10918230.3835360/151‑1 (MQ=60)
                                                                                                                     GCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTC                      >  SRR10918230.1312157/1‑151 (MQ=60)
                                                                                                                         GCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGA                  >  SRR10918230.4160991/1‑151 (MQ=60)
                                                                                                                           GACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATT                <  SRR10918230.1415012/151‑1 (MQ=60)
                                                                                                                           GACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATT                >  SRR10918230.3795293/1‑151 (MQ=60)
                                                                                                                            ACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTT               >  SRR10918230.4167727/1‑151 (MQ=60)
                                                                                                                             CGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTC              >  SRR10918230.1352049/1‑151 (MQ=60)
                                                                                                                              GAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCagatcggaagagcacacgtctgaa             >  SRR10918230.4482276/1‑127 (MQ=60)
                                                                                                                              GAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCagatcggaagagcacacgtctgaa             >  SRR10918230.4482220/1‑127 (MQ=60)
                                                                                                                              GAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCagatcggaagagcacacgtctgaa             >  SRR10918230.4482203/1‑127 (MQ=60)
                                                                                                                              GAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCG             >  SRR10918230.1213137/1‑151 (MQ=60)
                                                                                                                              GAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCG             >  SRR10918230.2648965/1‑151 (MQ=60)
                                                                                                                              GAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCG             <  SRR10918230.3365186/151‑1 (MQ=60)
                                                                                                                              GAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCG             <  SRR10918230.4381743/151‑1 (MQ=60)
                                                                                                                               AAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGG            <  SRR10918230.4215211/151‑1 (MQ=60)
                                                                                                                                 ATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCC          >  SRR10918230.203656/1‑151 (MQ=60)
                                                                                                                                   CAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAG        <  SRR10918230.1237728/151‑1 (MQ=60)
                                                                                                                                    AAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGG       <  SRR10918230.2938021/151‑1 (MQ=60)
                                                                                                                                    AAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGG       <  SRR10918230.2940366/151‑1 (MQ=60)
                                                                                                                                       ACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCC    >  SRR10918230.158976/1‑151 (MQ=60)
                                                                                                                                       ACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCC    >  SRR10918230.2458199/1‑151 (MQ=60)
                                                                                                                                       ACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCC    <  SRR10918230.723621/151‑1 (MQ=60)
                                                                                                                                        CCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCCG   >  SRR10918230.3887109/1‑151 (MQ=60)
                                                                                                                                         CATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCCGC  <  SRR10918230.3447488/151‑1 (MQ=60)
                                                                                                                                         CATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCCGC  <  SRR10918230.3449473/151‑1 (MQ=60)
                                                                                                                                                      |                                                                                                                                         
TCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGGCCCGC  >  CP009974/2738379‑2738666

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