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 F434 I0 R1 1130 215.9 7031662 97.1% 6827743 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,559,507 T→A 100% N162Y (AAT→TAT)  dosP ← oxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,559,5070TA98.5% 661.1 / ‑0.4 204N162Y (AAT→TAT) dosPoxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base A (98/103);  total (98/106)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 215 total aligned reads displayed.

TCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCC  >  CP009273/1559362‑1559650
                                                                                                                                                 |                                                                                                                                               
tCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATagc                                                                                                                                               >  1:965483/1‑148 (MQ=255)
 ccTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATagca                                                                                                                                              >  1:207532/1‑148 (MQ=255)
 ccTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATagca                                                                                                                                              <  2:1252186/148‑1 (MQ=255)
  cTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATagcac                                                                                                                                             <  1:2326414/148‑1 (MQ=255)
    gATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATagcactg                                                                                                                                           >  1:2840590/1‑148 (MQ=255)
        gcgggcggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGcaca                                                                                                                                       <  2:2840381/148‑1 (MQ=255)
         cgggcggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACaa                                                                                                                                      <  2:130516/148‑1 (MQ=255)
          gggcggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAat                                                                                                                                     <  1:264346/148‑1 (MQ=255)
          gggcggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAat                                                                                                                                     <  2:3198623/148‑1 (MQ=255)
             cggTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATTTTCAGGAGTGTATCGGTCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACAATTCGGTAAATTCGCGATAGCACAGCACAATATg                                                                                                                                  >  1:2101108/1‑148 (MQ=255)
                ttttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCg                                                                                                                               >  2:2510165/1‑148 (MQ=255)
                ttttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCg                                                                                                                               >  1:1911501/1‑148 (MQ=255)
                ttttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCg                                                                                                                               <  1:2055901/148‑1 (MQ=255)
                ttttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCg                                                                                                                               >  2:1461425/1‑148 (MQ=255)
                 tttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGt                                                                                                                              <  2:3491990/148‑1 (MQ=255)
                  ttCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGtt                                                                                                                             <  1:2146079/148‑1 (MQ=255)
                   tCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGTGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTc                                                                                                                            <  1:658016/148‑1 (MQ=255)
                    ccATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTcc                                                                                                                           >  2:3272428/1‑148 (MQ=255)
                      aTAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCgg                                                                                                                         <  1:1556757/148‑1 (MQ=255)
                       gaGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGa                                                                                                                        <  2:2405593/147‑1 (MQ=255)
                           aaCTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGa                                                                                                                    <  2:595830/148‑1 (MQ=255)
                                ttgtAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACa                                                                                                               >  1:2021511/1‑148 (MQ=255)
                                 tgtAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACaa                                                                                                              <  2:2688377/148‑1 (MQ=255)
                                  gtAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAAt                                                                                                             <  1:731120/148‑1 (MQ=255)
                                      aCGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTGTCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCAcc                                                                                                         >  1:30369/1‑148 (MQ=255)
                                       cGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcg                                                                                                        <  1:1536609/148‑1 (MQ=255)
                                         aaTGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggt                                                                                                      >  1:51246/1‑148 (MQ=255)
                                          aTGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtc                                                                                                     <  1:2944691/148‑1 (MQ=255)
                                            gCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtcgg                                                                                                   <  1:814762/148‑1 (MQ=255)
                                            gCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtcgg                                                                                                   >  1:1448189/1‑148 (MQ=255)
                                             cGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtcggt                                                                                                  >  1:2700692/1‑148 (MQ=255)
                                             cGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACcggtcggt                                                                                                  <  1:2728442/148‑1 (MQ=255)
                                                ttATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCga                                                                                               <  1:762048/148‑1 (MQ=255)
                                                 tATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCgag                                                                                              <  2:2563036/148‑1 (MQ=255)
                                                 tATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCgag                                                                                              <  1:1471002/148‑1 (MQ=255)
                                                 tATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCgag                                                                                              <  1:1627376/148‑1 (MQ=255)
                                                 tATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCgag                                                                                              <  1:3367653/148‑1 (MQ=255)
                                                  aTCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCgaga                                                                                             <  2:1913898/148‑1 (MQ=255)
                                                  aTCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCgaga                                                                                             <  2:1805800/148‑1 (MQ=255)
                                                  aTCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCgaga                                                                                             >  1:1241491/1‑148 (MQ=255)
                                                  aTCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCgaga                                                                                             <  2:1694370/148‑1 (MQ=255)
                                                   tCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGAt                                                                                            >  2:887245/1‑148 (MQ=255)
                                                      gCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGt                                                                                         >  1:3493140/1‑148 (MQ=255)
                                                         gggAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCaa                                                                                      >  1:1816661/1‑148 (MQ=255)
                                                          ggAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAAc                                                                                     >  2:234731/1‑148 (MQ=255)
                                                          ggAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAAc                                                                                     >  2:1098327/1‑148 (MQ=255)
                                                          ggAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAAc                                                                                     >  2:281341/1‑148 (MQ=255)
                                                                  ggAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCaataa                                                                             >  2:1313001/1‑148 (MQ=255)
                                                                   gAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCaataat                                                                            >  1:1752801/1‑148 (MQ=255)
                                                                     aTGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCa                                                                          <  1:330210/148‑1 (MQ=255)
                                                                      tGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCaa                                                                         <  1:72372/148‑1 (MQ=255)
                                                                      tGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCaa                                                                         <  1:2983742/148‑1 (MQ=255)
                                                                       gTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAAt                                                                        <  2:1072923/148‑1 (MQ=255)
                                                                        ttCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAAtt                                                                       <  2:3384497/148‑1 (MQ=255)
                                                                        ttCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAAtt                                                                       <  2:1865743/148‑1 (MQ=255)
                                                                               gtgtATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGt                                                                <  2:243505/148‑1 (MQ=255)
                                                                                  tATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTg                                                             >  2:2960794/1‑148 (MQ=255)
                                                                                   aTCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCACACATTTCGTTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGAAACTCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGt                                                            <  2:116802/148‑1 (MQ=255)
                                                                                    tCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGtt                                                           <  1:1314558/148‑1 (MQ=255)
                                                                                     cGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTc                                                          >  2:2624765/1‑148 (MQ=255)
                                                                                        gCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCttt                                                       <  2:1065536/148‑1 (MQ=255)
                                                                                        gCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCttt                                                       <  2:645159/148‑1 (MQ=255)
                                                                                        gCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCttt                                                       >  1:550804/1‑148 (MQ=255)
                                                                                           gCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTg                                                    >  1:3093425/1‑148 (MQ=255)
                                                                                              tACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTgcgc                                                 >  2:3088853/1‑148 (MQ=255)
                                                                                                ccGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCa                                               <  2:2315570/148‑1 (MQ=255)
                                                                                                   cTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCAttt                                            <  2:705121/148‑1 (MQ=255)
                                                                                                   cTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCAttt                                            <  2:2463476/148‑1 (MQ=255)
                                                                                                    tGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTc                                           <  2:2232196/148‑1 (MQ=255)
                                                                                                      gCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTa                                         >  2:413281/1‑148 (MQ=255)
                                                                                                         tCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGc                                      <  1:1031179/148‑1 (MQ=255)
                                                                                                            cTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTgg                                   >  1:2514276/1‑148 (MQ=255)
                                                                                                            cTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTgg                                   >  1:1828556/1‑148 (MQ=255)
                                                                                                             tAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGc                                  >  2:1895497/1‑148 (MQ=255)
                                                                                                              aaTGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCa                                 >  2:2976482/1‑148 (MQ=255)
                                                                                                               aTGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCAt                                >  2:3415609/1‑148 (MQ=255)
                                                                                                               aTGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCAt                                >  2:3415541/1‑148 (MQ=255)
                                                                                                                tGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATc                               >  1:497699/1‑148 (MQ=255)
                                                                                                                 gCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATcc                              >  2:1622276/1‑148 (MQ=255)
                                                                                                                    gTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGt                           >  2:957433/1‑148 (MQ=255)
                                                                                                                    gTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGt                           <  1:2485497/148‑1 (MQ=255)
                                                                                                                     tAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTa                          >  2:10767/1‑148 (MQ=255)
                                                                                                                     tAACCAAACATTCCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATTCCGTa                          >  2:1357978/1‑148 (MQ=255)
                                                                                                                       aCCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTAcc                        >  2:2682061/1‑148 (MQ=255)
                                                                                                                        ccAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCa                       >  2:303164/1‑148 (MQ=255)
                                                                                                                        ccAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCa                       >  1:2296477/1‑148 (MQ=255)
                                                                                                                         cAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAg                      >  2:211816/1‑148 (MQ=255)
                                                                                                                          aaaCATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAgc                     >  2:1964327/1‑148 (MQ=255)
                                                                                                                           aaCATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAgcg                    >  1:1963578/1‑148 (MQ=255)
                                                                                                                            aCATTTCGGTAAATGCGCGATAGCACTGCACGATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGCTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAgcgc                   <  1:165853/148‑1 (MQ=255)
                                                                                                                            aCATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAgcgc                   <  2:1658001/148‑1 (MQ=255)
                                                                                                                             cATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGcc                  <  1:481389/148‑1 (MQ=255)
                                                                                                                             cATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGcc                  <  1:2198259/148‑1 (MQ=255)
                                                                                                                                ttCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCACTAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAgg               >  2:276602/1‑148 (MQ=255)
                                                                                                                                   ggTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGtaa            <  2:2891557/148‑1 (MQ=255)
                                                                                                                                     tAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGtaata          <  1:591545/148‑1 (MQ=255)
                                                                                                                                       aaTGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATaaa        >  2:342884/1‑148 (MQ=255)
                                                                                                                                          gcgcGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACtt     <  1:2118732/148‑1 (MQ=255)
                                                                                                                                          gcgcGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACtt     <  1:1734139/148‑1 (MQ=255)
                                                                                                                                             cgATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTcc  >  2:2964324/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCC  >  CP009273/1559362‑1559650

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCG  >  CP009273/1559358‑1559655
                                                                                                                                                     |                                                                                                                                                    
TTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATA                                                                                                                                                      >  SRR19973065.709817/1‑150 (MQ=60)
TTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATA                                                                                                                                                      >  SRR19973065.2479003/1‑150 (MQ=60)
  CGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGC                                                                                                                                                    >  SRR19973065.971891/1‑150 (MQ=60)
  CGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGC                                                                                                                                                    <  SRR19973065.2446566/150‑1 (MQ=60)
  CGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGC                                                                                                                                                    <  SRR19973065.1446919/150‑1 (MQ=60)
  CGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGC                                                                                                                                                    <  SRR19973065.1446903/150‑1 (MQ=60)
   GTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCA                                                                                                                                                   >  SRR19973065.209107/1‑150 (MQ=60)
      CTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTG                                                                                                                                                >  SRR19973065.2858078/1‑150 (MQ=60)
      CTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTG                                                                                                                                                >  SRR19973065.2726685/1‑150 (MQ=60)
      CTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTG                                                                                                                                                <  SRR19973065.2340973/150‑1 (MQ=60)
       TGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGC                                                                                                                                               >  SRR19973065.3043701/1‑150 (MQ=60)
        GATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCA                                                                                                                                              <  SRR19973065.2524395/150‑1 (MQ=60)
              GGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATAT                                                                                                                                        <  SRR19973065.266157/150‑1 (MQ=60)
               GGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATTTTCAGGAGTGTATCGGTCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACAATTCGGTAAATTCGCGATAGCACAGCACAATATG                                                                                                                                       >  SRR19973065.2114621/1‑150 (MQ=60)
               GGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATG                                                                                                                                       >  SRR19973065.3127053/1‑150 (MQ=60)
                GCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGG                                                                                                                                      <  SRR19973065.2980857/150‑1 (MQ=60)
                GCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGG                                                                                                                                      <  SRR19973065.2288352/150‑1 (MQ=60)
                  GGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCG                                                                                                                                    >  SRR19973065.1924052/1‑150 (MQ=60)
                  GGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCG                                                                                                                                    <  SRR19973065.270790/150‑1 (MQ=60)
                    TTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTT                                                                                                                                  <  SRR19973065.2069028/150‑1 (MQ=60)
                    TTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTT                                                                                                                                  <  SRR19973065.2069043/150‑1 (MQ=60)
                      TTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCC                                                                                                                                <  SRR19973065.2159880/150‑1 (MQ=60)
                       TCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGTGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCG                                                                                                                               <  SRR19973065.662361/150‑1 (MQ=60)
                       TCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCG                                                                                                                               <  SRR19973065.3267032/150‑1 (MQ=60)
                          ATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGAT                                                                                                                            <  SRR19973065.1567087/150‑1 (MQ=60)
                                  TGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACA                                                                                                                    >  SRR19973065.2034516/1‑150 (MQ=60)
                                      GTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCA                                                                                                                <  SRR19973065.736016/150‑1 (MQ=60)
                                       TAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCAC                                                                                                               <  SRR19973065.1416887/150‑1 (MQ=60)
                                        AAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTGTCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACC                                                                                                              >  SRR19973065.30753/1‑150 (MQ=60)
                                        AAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACC                                                                                                              >  SRR19973065.533794/1‑150 (MQ=60)
                                        AAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACC                                                                                                              <  SRR19973065.2037516/150‑1 (MQ=60)
                                           CGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGT                                                                                                           >  SRR19973065.51857/1‑150 (MQ=60)
                                           CGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGT                                                                                                           >  SRR19973065.51814/1‑150 (MQ=60)
                                           CGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGT                                                                                                           <  SRR19973065.1546863/150‑1 (MQ=60)
                                           CGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGT                                                                                                           <  SRR19973065.1546744/150‑1 (MQ=60)
                                             AATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCG                                                                                                         <  SRR19973065.586271/150‑1 (MQ=60)
                                              ATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGG                                                                                                        <  SRR19973065.2962585/150‑1 (MQ=60)
                                              ATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGG                                                                                                        >  SRR19973065.1457544/1‑150 (MQ=60)
                                               TGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGT                                                                                                       >  SRR19973065.2717470/1‑150 (MQ=60)
                                                GCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTC                                                                                                      <  SRR19973065.666525/150‑1 (MQ=60)
                                                GCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTC                                                                                                      <  SRR19973065.820183/150‑1 (MQ=60)
                                                GCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTC                                                                                                      <  SRR19973065.820131/150‑1 (MQ=60)
                                                 CGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCG                                                                                                     <  SRR19973065.2745412/150‑1 (MQ=60)
                                                  GGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGA                                                                                                    >  SRR19973065.166629/1‑150 (MQ=60)
                                                    TTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGA                                                                                                  >  SRR19973065.1249655/1‑150 (MQ=60)
                                                    TTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGA                                                                                                  <  SRR19973065.767124/150‑1 (MQ=60)
                                                     TATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGAT                                                                                                 <  SRR19973065.3388245/150‑1 (MQ=60)
                                                     TATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGAT                                                                                                 <  SRR19973065.641277/150‑1 (MQ=60)
                                                     TATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGAT                                                                                                 <  SRR19973065.1480429/150‑1 (MQ=60)
                                                     TATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGAT                                                                                                 <  SRR19973065.1638283/150‑1 (MQ=60)
                                                        CGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGT                                                                                              >  SRR19973065.3514376/1‑150 (MQ=60)
                                                         GGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTC                                                                                             >  SRR19973065.2632248/1‑150 (MQ=60)
                                                           CAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAA                                                                                           >  SRR19973065.1828625/1‑150 (MQ=60)
                                                           CAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAA                                                                                           >  SRR19973065.1830929/1‑150 (MQ=60)
                                                                AATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCA                                                                                      >  SRR19973065.771086/1‑150 (MQ=60)
                                                                AATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCA                                                                                      >  SRR19973065.168499/1‑150 (MQ=60)
                                                                     AGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAAT                                                                                 >  SRR19973065.1764338/1‑150 (MQ=60)
                                                                       GAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCA                                                                               <  SRR19973065.1823530/150‑1 (MQ=60)
                                                                        AATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCACTAATCAA                                                                              >  SRR19973065.502646/1‑150 (MQ=60)
                                                                         ATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAAT                                                                             <  SRR19973065.332521/150‑1 (MQ=60)
                                                                          TGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATT                                                                            <  SRR19973065.73128/150‑1 (MQ=60)
                                                                          TGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATT                                                                            <  SRR19973065.3002007/150‑1 (MQ=60)
                                                                              CAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCG                                                                        >  SRR19973065.1746092/1‑150 (MQ=60)
                                                                               AGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGG                                                                       <  SRR19973065.2120731/150‑1 (MQ=60)
                                                                                 GAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGT                                                                     >  SRR19973065.3483149/1‑150 (MQ=60)
                                                                                  AGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTC                                                                    >  SRR19973065.3357800/1‑150 (MQ=60)
                                                                                  AGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTC                                                                    <  SRR19973065.1470817/150‑1 (MQ=60)
                                                                                      TATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTT                                                                >  SRR19973065.2152771/1‑150 (MQ=60)
                                                                                        TCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCT                                                              <  SRR19973065.2092694/150‑1 (MQ=60)
                                                                                        TCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCT                                                              <  SRR19973065.1323250/150‑1 (MQ=60)
                                                                                          GGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTT                                                            >  SRR19973065.554346/1‑150 (MQ=60)
                                                                                            GCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTT                                                          <  SRR19973065.2141582/150‑1 (MQ=60)
                                                                                             CTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTG                                                         >  SRR19973065.3112312/1‑150 (MQ=60)
                                                                                               GCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCG                                                       <  SRR19973065.2663391/150‑1 (MQ=60)
                                                                                                CATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGC                                                      >  SRR19973065.1846374/1‑150 (MQ=60)
                                                                                                     CGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTT                                                 >  SRR19973065.1932931/1‑150 (MQ=60)
                                                                                                     CGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTT                                                 <  SRR19973065.117411/150‑1 (MQ=60)
                                                                                                     CGATGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACCGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTT                                                 >  SRR19973065.1933004/1‑150 (MQ=60)
                                                                                                         GGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTAC                                             <  SRR19973065.3079456/150‑1 (MQ=60)
                                                                                                             TCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTG                                         <  SRR19973065.1038010/150‑1 (MQ=60)
                                                                                                              CGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGG                                        >  SRR19973065.758011/1‑150 (MQ=60)
                                                                                                              CGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGG                                        >  SRR19973065.2530074/1‑150 (MQ=60)
                                                                                                              CGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGG                                        >  SRR19973065.1840601/1‑150 (MQ=60)
                                                                                                               GCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGC                                       <  SRR19973065.2864866/150‑1 (MQ=60)
                                                                                                                 TAATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCAT                                     >  SRR19973065.1513457/1‑150 (MQ=60)
                                                                                                                  AATGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATC                                    >  SRR19973065.501040/1‑150 (MQ=60)
                                                                                                                    TGCAGTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCC                                  <  SRR19973065.2176873/150‑1 (MQ=60)
                                                                                                                        GTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTAC                              <  SRR19973065.730540/150‑1 (MQ=60)
                                                                                                                        GTAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTAC                              <  SRR19973065.2501134/150‑1 (MQ=60)
                                                                                                                         TAACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACC                             >  SRR19973065.44867/1‑150 (MQ=60)
                                                                                                                          AACCAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCA                            >  SRR19973065.2310938/1‑150 (MQ=60)
                                                                                                                             CAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCG                         >  SRR19973065.1976654/1‑150 (MQ=60)
                                                                                                                             CAAACATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCG                         >  SRR19973065.1976396/1‑150 (MQ=60)
                                                                                                                                ACATTTCGGTAAATGCGCGATAGCACTGCACGATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGCTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCA                      <  SRR19973065.167207/150‑1 (MQ=60)
                                                                                                                                 CATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAG                     <  SRR19973065.2212337/150‑1 (MQ=60)
                                                                                                                                 CATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAG                     <  SRR19973065.484751/150‑1 (MQ=60)
                                                                                                                                 CATTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAG                     <  SRR19973065.484666/150‑1 (MQ=60)
                                                                                                                                   TTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGT                   >  SRR19973065.1575211/1‑150 (MQ=60)
                                                                                                                                   GTTCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCAACGGTCGGTCGATATGGTCAACGGCAATAATCAATTGTCGGGTCTGTACTTTTTGCTCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGT                   >  SRR19973065.1575215/1‑150 (MQ=60)
                                                                                                                                     TCGGTAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAAGCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAA                 >  SRR19973065.2828197/1‑150 (MQ=60)
                                                                                                                                         TAAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAA             <  SRR19973065.595206/150‑1 (MQ=60)
                                                                                                                                          AAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTCATAAAC            <  SRR19973065.231918/150‑1 (MQ=60)
                                                                                                                                          AAATGCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAAC            <  SRR19973065.1297447/150‑1 (MQ=60)
                                                                                                                                              GCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTC        <  SRR19973065.1745614/150‑1 (MQ=60)
                                                                                                                                              GCGCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTC        <  SRR19973065.2132365/150‑1 (MQ=60)
                                                                                                                                                GCGATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCC      <  SRR19973065.971828/150‑1 (MQ=60)
                                                                                                                                                 CGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCC     <  SRR19973065.2423627/150‑1 (MQ=60)
                                                                                                                                                   ATAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTC   >  SRR19973065.663052/1‑150 (MQ=60)
                                                                                                                                                    TAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCG  <  SRR19973065.1172473/150‑1 (MQ=60)
                                                                                                                                                    TAGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCG  <  SRR19973065.2690129/150‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                                    
TTCGTCCTGATCGCGGGCGGTTTTCCATAGCAACTGTTGTAAACGAATGCGGTTATCGGCAGGGAATTCAGGAATGTTCAGGAGTGTATCGGGCTGCATACCGCTGGCTTCGCTAATGCAGTAACCAAACATTTCGGTAAATGCGCGATTGCACTGCACAATATGGCGTTCCGGATCGAGGACAATCACCGGTCGGTCGAGATGGTCAACGGCAATAATCAATTGTCGGGTCTGTTCTTTTTGCGCCATTTCTACGCTGGCATCCCGTACCAGCGCCAGGTAATAAACTTTCCCCTCG  >  CP009273/1559358‑1559655

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