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,513 G→A K242K (AAG→AAA RPPX_12165 → benzene 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099742,738,5130GA99.4% 349.2 / ‑5.0 157K242K (AAG→AAARPPX_12165benzene 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (80/75);  total (80/77)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 157 total aligned reads displayed.

CGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAG  >  CP009974/2738369‑2738660
                                                                                                                                                |                                                                                                                                                   
tggagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATg                                                                                                                                               <  1:2678076/144‑1 (MQ=1)
  gagctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAg                                                                                                                                             >  1:423224/6‑151 (MQ=1)
   agctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgc                                                                                                                                            >  2:4290040/5‑150 (MQ=1)
    gctaGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAgcg                                                                                                                                           >  1:3967119/4‑149 (MQ=1)
                   accggaGGCAACTGGAAGCTCCCCGCCGAAACCGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGAAACTACGCTGCCACCCAGAACCAGCGAAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCa                                                                                                                            <  2:1457271/146‑1 (MQ=1)
                         ggtaaCTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAgg                                                                                                                      >  2:1779345/4‑151 (MQ=0)
                                  aaGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGtt                                                                                                             >  2:2054432/1‑151 (MQ=0)
                                     ctgaCCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCt                                                                                                          <  2:2876788/148‑1 (MQ=0)
                                                   cGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCAt                                                                                            <  1:4204967/151‑1 (MQ=2)
                                                    ggCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATg                                                                                           <  1:3496694/151‑1 (MQ=2)
                                                     gcgcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATgg                                                                                          >  1:5226025/1‑151 (MQ=2)
                                                     gcgcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATgg                                                                                          >  2:4214548/1‑151 (MQ=11)
                                                     gcgcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATgg                                                                                          <  2:3370886/151‑1 (MQ=2)
                                                       gcCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGac                                                                                        >  2:3806900/1‑151 (MQ=2)
                                                        ccGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGaca                                                                                       <  1:743505/151‑1 (MQ=2)
                                                           aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                                                                    >  2:2646067/1‑151 (MQ=2)
                                                           aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                                                                    >  1:2596495/1‑151 (MQ=2)
                                                           aCGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACct                                                                                    >  2:668703/1‑151 (MQ=2)
                                                           aCGGCTATCATGTTAGCGCGGTGCACTGGAACTACGCCGCTAATCGCAACCAGCGCAAGCAGCGCGAGGCAGGCGACTAAGTTAAATCCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTTGTTGTTTCTATTCCTTCGAGCATGGACTCct                                                                                    <  2:3041600/151‑1 (MQ=2)
                                                            cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctg                                                                                   >  1:2506923/1‑151 (MQ=11)
                                                            cGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGGTTCTATTCCTTCCACCATGGACACctg                                                                                   >  2:4326824/1‑151 (MQ=2)
                                                             ggCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACctgc                                                                                  >  2:493223/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)
                                                                   catgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGctct                                                                            <  1:1415458/145‑1 (MQ=2)
                                                                    atgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTg                                                                           >  2:4730553/6‑151 (MQ=11)
                                                                     tgttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTgg                                                                          <  1:3256087/147‑1 (MQ=11)
                                                                       ttaGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGAc                                                                        <  2:2715610/149‑1 (MQ=11)
                                                                           cTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCgg                                                                    >  2:1838182/1‑151 (MQ=14)
                                                                             ccGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTg                                                                  <  2:1923450/151‑1 (MQ=14)
                                                                             ccGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTg                                                                  <  1:2099453/151‑1 (MQ=14)
                                                                              cGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTgg                                                                 <  1:104673/151‑1 (MQ=14)
                                                                              cGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTgg                                                                 <  1:105033/151‑1 (MQ=14)
                                                                              cGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTgg                                                                 >  1:4538964/1‑151 (MQ=14)
                                                                                   aCTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGaa                                                            <  2:5313792/151‑1 (MQ=14)
                                                                                   aCTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGaa                                                            >  1:4736720/1‑151 (MQ=14)
                                                                                   aCTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGaa                                                            >  1:2164522/1‑151 (MQ=14)
                                                                                     tGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcc                                                          >  2:2817590/1‑151 (MQ=14)
                                                                                      ggAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAccc                                                         <  1:793839/151‑1 (MQ=17)
                                                                                       gAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAAcccc                                                        <  2:1385394/151‑1 (MQ=17)
                                                                                          cTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAg                                                     <  2:219614/151‑1 (MQ=17)
                                                                                          cTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAg                                                     >  2:3220652/1‑151 (MQ=17)
                                                                                           tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                    <  1:4515714/151‑1 (MQ=17)
                                                                                           tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                    >  2:3268979/1‑151 (MQ=17)
                                                                                           tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                    <  2:1714847/151‑1 (MQ=17)
                                                                                           tACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAgg                                                    <  1:964545/151‑1 (MQ=17)
                                                                                            aCGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGa                                                   >  2:31300/1‑151 (MQ=17)
                                                                                             cgccgcCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGAt                                                  <  1:4949585/151‑1 (MQ=17)
                                                                                                cgCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTTCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAggtttgc                                               >  1:408337/2‑146 (MQ=17)
                                                                                                cgCCACTCAGAACCAACGCAAGCAGCGCGAGGCAGGCGACGAAAGCAAAACCAGGAGCGCCGGCTCTTAGGCTAAGCAGGGTGGTGGTGTCTATTCCTTCGACCAGGGACGCCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGc                                               <  2:3620336/150‑1 (MQ=255)
                                                                                                 gCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGcc                                              <  1:890926/151‑1 (MQ=17)
                                                                                                  ccACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGccc                                             >  1:786268/1‑151 (MQ=17)
                                                                                                  ccACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGccc                                             >  1:514803/1‑151 (MQ=17)
                                                                                                    actcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAg                                           >  2:1165595/4‑151 (MQ=17)
                                                                                                     ctcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGc                                          >  1:3522604/3‑151 (MQ=18)
                                                                                                     ctcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGc                                          >  1:1429328/3‑151 (MQ=18)
                                                                                                      tcAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCt                                         <  2:368254/150‑1 (MQ=18)
                                                                                                             cAGCGCCAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGGGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAgc                                  <  2:3863099/151‑1 (MQ=255)
                                                                                                              aGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAgcg                                 <  2:2401391/151‑1 (MQ=18)
                                                                                                                cgcAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                               <  2:834715/151‑1 (MQ=255)
                                                                                                                 gcAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCg                              >  1:4393400/1‑151 (MQ=18)
                                                                                                                  cAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGt                             >  2:4510787/1‑151 (MQ=255)
                                                                                                                   aaGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTg                            <  1:592820/151‑1 (MQ=255)
                                                                                                                    agcagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGa                           >  2:1731813/1‑151 (MQ=255)
                                                                                                                    agcagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGa                           >  1:4316570/1‑151 (MQ=255)
                                                                                                                      cagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAcc                         <  1:5143558/151‑1 (MQ=18)
                                                                                                                      cagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGAcc                         <  1:4494949/151‑1 (MQ=18)
                                                                                                                      cagcGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGCATCGCCCAGCTTATGAGCGCCGTGAcc                         >  2:2738838/1‑151 (MQ=255)
                                                                                                                        gcgcGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCaa                       >  1:5196064/1‑151 (MQ=255)
                                                                                                                        gcgcGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCaa                       <  1:22896/151‑1 (MQ=255)
                                                                                                                         cgcgAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCTGCTTATGAGCGCCGTGACCAAt                      <  2:4831024/151‑1 (MQ=255)
                                                                                                                           cgAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTg                    >  2:2850142/1‑151 (MQ=255)
                                                                                                                             aggcaggcGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGc                  <  1:144461/151‑1 (MQ=255)
                                                                                                                              ggcgggcgACGAAATCAAAACCAGGAGCGCCGGCTCGTGGGCTAAGCAGGGTGGTGGGTTCTAGTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCt                 <  2:5143999/151‑1 (MQ=255)
                                                                                                                              ggcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCt                 <  2:2036425/144‑1 (MQ=21)
                                                                                                                              ggcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCt                 >  1:3288601/8‑151 (MQ=255)
                                                                                                                              ggcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCt                 <  1:3818708/144‑1 (MQ=21)
                                                                                                                               gcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTc                >  1:1306264/7‑151 (MQ=255)
                                                                                                                               gcaggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTc                >  2:3271331/7‑151 (MQ=255)
                                                                                                                                caggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCa               <  1:727060/146‑1 (MQ=255)
                                                                                                                                  ggcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTGCGATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAg             <  2:2685210/148‑1 (MQ=255)
                                                                                                                                   gcgACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGa            >  1:4142511/3‑151 (MQ=255)
                                                                                                                                     gACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAtt          >  1:3778863/1‑151 (MQ=21)
                                                                                                                                     gACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAtt          >  2:2507410/1‑151 (MQ=25)
                                                                                                                                     gACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGAtt          <  1:1408691/151‑1 (MQ=255)
                                                                                                                                        gaaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                               >  1:4462035/4‑127 (MQ=21)
                                                                                                                                        gaaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                               <  2:4462052/124‑1 (MQ=255)
                                                                                                                                        gaaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                               <  2:4462108/124‑1 (MQ=255)
                                                                                                                                        gaaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGcc                               >  1:4462108/4‑127 (MQ=21)
                                                                                                                                        gaaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCg       >  1:2637477/4‑151 (MQ=255)
                                                                                                                                        gaaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCg       <  2:2642600/148‑1 (MQ=255)
                                                                                                                                        gaaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGAGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTCGGCGATCCCCGAGGACCGCGCAGCTTATGAGCGCCGTGCCCACTTGGCTCACGATTTCg       >  2:4632489/4‑151 (MQ=255)
                                                                                                                                         aaaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCgg      <  1:4196428/149‑1 (MQ=255)
                                                                                                                                          aaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGc     >  2:5044641/2‑151 (MQ=255)
                                                                                                                                          aaTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGc     <  2:3967119/150‑1 (MQ=25)
                                                                                                                                           aTCAAAACCATGCGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGcc    <  1:424628/151‑1 (MQ=255)
                                                                                                                                           aTCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGcc    >  1:202987/1‑151 (MQ=255)
                                                                                                                                            tCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCa   <  2:2311334/151‑1 (MQ=255)
                                                                                                                                             cAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAg  <  1:1232330/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAG  >  CP009974/2738369‑2738660

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
CTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGG  >  CP009974/2738364‑2738661
                                                                                                                                                     |                                                                                                                                                    
CTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAA                                                                                                                                                     >  SRR10918230.3284396/1‑151 (MQ=0)
CTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAA                                                                                                                                                     <  SRR10918230.3112913/151‑1 (MQ=0)
CTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAA                                                                                                                                                     >  SRR10918230.1973901/1‑151 (MQ=0)
CTGCGTGGAGCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAA                                                                                                                                                     >  SRR10918230.1971459/1‑151 (MQ=0)
CTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGA                                                                                                                                                     <  SRR10918230.1294778/151‑1 (MQ=0)
 TGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGAC                                                                                                                                                    <  SRR10918230.4598423/151‑1 (MQ=0)
     tggagctAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGACCTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATG                                                                                                                                                <  SRR10918230.1427261/144‑1 (MQ=0)
     tggagctAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATG                                                                                                                                                <  SRR10918230.2689783/144‑1 (MQ=0)
     tggagctAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATG                                                                                                                                                >  SRR10918230.2599307/8‑151 (MQ=0)
     tggagctAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATG                                                                                                                                                <  SRR10918230.637805/144‑1 (MQ=0)
      GGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGA                                                                                                                                               >  SRR10918230.281180/1‑151 (MQ=0)
       gagctAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAG                                                                                                                                              >  SRR10918230.424685/6‑151 (MQ=0)
         GCTAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCG                                                                                                                                            >  SRR10918230.3984559/1‑151 (MQ=0)
          CCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGC                                                                                                                                           >  SRR10918230.986859/1‑151 (MQ=0)
           TAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATagatcgg                                                                                                                                          >  SRR10918230.3875351/1‑144 (MQ=0)
           TAGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGcgcc                                                                                                                                          >  SRR10918230.2037760/1‑147 (MQ=0)
            AGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTG                                                                                                                                         <  SRR10918230.4536937/151‑1 (MQ=0)
            AGCTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGcgccg                                                                                                                                         >  SRR10918230.3310811/1‑146 (MQ=0)
              CTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGC                                                                                                                                       <  SRR10918230.3114769/151‑1 (MQ=0)
              CTCCTACATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGC                                                                                                                                       >  SRR10918230.3669140/1‑151 (MQ=0)
               TCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCG                                                                                                                                      >  SRR10918230.3891025/1‑151 (MQ=0)
               TCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCG                                                                                                                                      <  SRR10918230.2816143/151‑1 (MQ=0)
               TCGTACATCTACGAAGGCAACTAGAAGCTCACCGCCCAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCG                                                                                                                                      <  SRR10918230.1060081/151‑1 (MQ=0)
                gtgctcttccgatctgtAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGcgccggctc                                                                                                                                     <  SRR10918230.537558/134‑10 (MQ=0)
                    CATCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGcgccggctcttgg                                                                                                                                 <  SRR10918230.1567482/151‑14 (MQ=0)
                    CATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGG                                                                                                                                 >  SRR10918230.4252741/1‑151 (MQ=0)
                      TCTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGcgccggctcttgggc                                                                                                                               <  SRR10918230.4889530/151‑16 (MQ=0)
                       CTACGAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGcgccggctcttgggct                                                                                                                              <  SRR10918230.1350092/151‑17 (MQ=0)
                           GAAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGcgccggctcttgggctaagc                                                                                                                          <  SRR10918230.1258987/151‑21 (MQ=0)
                            AAGGTAACTGGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCA                                                                                                                         <  SRR10918230.1375586/151‑1 (MQ=0)
                            AAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCA                                                                                                                         <  SRR10918230.4805898/151‑1 (MQ=0)
                            AAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCA                                                                                                                         >  SRR10918230.502404/1‑151 (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)
                                     GGAAGCTGACCGCCGAGAACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGcgccggctcttgggctaagcagggtggtgg                                                                                                                <  SRR10918230.1274190/151‑31 (MQ=0)
                                        AGCTGACCGCCGAGCACGGCGCCGACGGCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGcggcggctcttgggctaagcagggtggtggttt                                                                                                             <  SRR10918230.4309338/151‑34 (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.4535676/1‑151 (MQ=14)
                                                                   GCTATCATGTTAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCT                                                                                  <  SRR10918230.1749834/151‑1 (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.4880980/144‑1 (MQ=17)
                                                                       tcatgttAGCTCCGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTC                                                                              <  SRR10918230.3625773/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.105024/151‑1 (MQ=33)
                                                                                   CGTGCACTGGAACTACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGG                                                                  >  SRR10918230.4559480/1‑151 (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.968426/151‑1 (MQ=53)
                                                                                                TACGCCGCCACTCAGAACCAGCGCAAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGG                                                     <  SRR10918230.4536127/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.594933/151‑1 (MQ=60)
                                                                                                                        AAGCAGCGCGAGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTG                             >  SRR10918230.2563356/1‑151 (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.144965/151‑1 (MQ=60)
                                                                                                                                  AGGCAGGCGACGAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGC                   <  SRR10918230.144823/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.2648965/1‑151 (MQ=60)
                                                                                                                                             GAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCG        >  SRR10918230.1213137/1‑151 (MQ=60)
                                                                                                                                             GAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCG        <  SRR10918230.4381743/151‑1 (MQ=60)
                                                                                                                                             GAAATCAAAACCATGAGCGCCGGCTCTTGGGCTAAGCAGGGTGGTGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCG        <  SRR10918230.3365186/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)
                                                                                                                                                     |                                                                                                                                                    
CTGCGCGGTGCCAGCTCGTACATCTACGAAGGCAACTGGAAGCTCACCGCCGAAAACGGCGCCGACGGCTACCACGTAAGCTCCGTGCACTGGAACTACGCTGCCACCCAGAACCAGCGCAAGCAGCGCGAGGCGGGTGACGAGATCAAGACCATGAGTGCTGGCGCCTGGGCCAAGCAGGGCGGCGGTTTCTATTCCTTCGACCATGGACACCTGCTGCTCTGGACTCGGTGGGCGAACCCCGAGGATCGCCCAGCTTATGAGCGCCGTGACCAATTGGCTCAAGATTTCGGCCAGG  >  CP009974/2738364‑2738661

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