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 F342 I0 R1 1331 225.9 7641038 95.8% 7320114 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,664,138 A→C 49.4% E130D (GAA→GAC)  hcaF → 3‑phenylpropionate dioxygenase, small (beta) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,664,1380AC49.4% 93.9 / inf 243E130D (GAA→GAC) hcaF3‑phenylpropionate dioxygenase, small (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (67/56);  new base C (59/61);  total (126/117)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 246 total aligned reads displayed.

AAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGT  >  CP009273/2663993‑2664285
                                                                                                                                                 |                                                                                                                                                   
    gggCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTc                                                                                                                                               <  1:2570928/148‑1 (MQ=255)
    gggCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTc                                                                                                                                               <  2:2327814/148‑1 (MQ=255)
     ggCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCt                                                                                                                                              <  2:3811661/148‑1 (MQ=255)
     ggCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCt                                                                                                                                              <  2:2469327/148‑1 (MQ=255)
      gCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTa                                                                                                                                             <  1:3598973/148‑1 (MQ=255)
             cTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGa                                                                                                                                      <  1:762614/148‑1 (MQ=255)
              tGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGaa                                                                                                                                     <  2:3025232/148‑1 (MQ=255)
                ggCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAAcg                                                                                                                                   >  1:3399053/1‑148 (MQ=255)
                  cAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAAcgcg                                                                                                                                 <  1:3406948/148‑1 (MQ=255)
                     aagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGttt                                                                                                                              <  2:944780/148‑1 (MQ=255)
                     aagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGttt                                                                                                                              <  1:360489/148‑1 (MQ=255)
                     aagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGttt                                                                                                                              <  1:3738194/148‑1 (MQ=255)
                        agCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGa                                                                                                                           >  2:2014186/1‑148 (MQ=255)
                          ccgccgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACa                                                                                                                         >  1:3478054/1‑148 (MQ=255)
                           cgccgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACaa                                                                                                                        <  2:178086/148‑1 (MQ=255)
                                tCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTc                                                                                                                   >  2:336921/1‑148 (MQ=255)
                                     cACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGt                                                                                                              >  2:3779734/1‑148 (MQ=255)
                                      aCCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTc                                                                                                             <  2:238600/148‑1 (MQ=255)
                                       cccGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCt                                                                                                            >  1:2705639/1‑148 (MQ=255)
                                          gTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGa                                                                                                         <  1:957780/148‑1 (MQ=255)
                                            cACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAg                                                                                                       >  2:2061986/1‑148 (MQ=255)
                                            cACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAg                                                                                                       >  2:3431617/1‑148 (MQ=255)
                                            cACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAg                                                                                                       >  2:1841210/1‑148 (MQ=255)
                                             aCTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGa                                                                                                      <  1:720586/148‑1 (MQ=255)
                                               ttAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATg                                                                                                    <  1:3675021/148‑1 (MQ=255)
                                                tAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGa                                                                                                   >  2:2584392/1‑148 (MQ=255)
                                                tAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGa                                                                                                   >  1:2029686/1‑148 (MQ=255)
                                                   tCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACaa                                                                                                <  1:1896769/148‑1 (MQ=255)
                                                   tCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACaa                                                                                                >  1:922928/1‑148 (MQ=255)
                                                       cAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTgg                                                                                            >  1:1296899/1‑148 (MQ=255)
                                                        aaCTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGgc                                                                                           >  1:3464161/1‑148 (MQ=255)
                                                         aCTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGgcg                                                                                          >  1:791753/1‑148 (MQ=255)
                                                         aCTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGgcg                                                                                          >  2:3694525/1‑148 (MQ=255)
                                                          cTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGgcgc                                                                                         <  2:2308853/148‑1 (MQ=255)
                                                          cTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAATGTTCGCCGTCTGGAAGATGACAACTGgcgc                                                                                         >  1:3012246/1‑148 (MQ=255)
                                                          cTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGgcgc                                                                                         >  1:2878094/1‑148 (MQ=255)
                                                           tGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCt                                                                                        >  2:1123634/1‑148 (MQ=255)
                                                             ccAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTg                                                                                      <  2:114884/148‑1 (MQ=255)
                                                             ccAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTg                                                                                      <  2:48116/148‑1 (MQ=255)
                                                                 aTAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTgg                                                                                  >  1:3227874/1‑148 (MQ=255)
                                                                  tAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGa                                                                                 >  1:3185438/1‑148 (MQ=255)
                                                                   aaGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGaa                                                                                >  2:1656512/1‑148 (MQ=255)
                                                                      cGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACgg                                                                             >  1:3454846/1‑148 (MQ=255)
                                                                        aaaCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGa                                                                           <  1:1666273/148‑1 (MQ=255)
                                                                        aaaCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGa                                                                           <  1:2791040/148‑1 (MQ=255)
                                                                        aaaCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCCGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGa                                                                           >  2:935759/1‑148 (MQ=255)
                                                                          aCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGata                                                                         >  2:2243420/1‑148 (MQ=255)
                                                                           ccGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGatat                                                                        <  1:2315383/148‑1 (MQ=255)
                                                                               cATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                    <  2:2262224/148‑1 (MQ=255)
                                                                               cATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                    >  1:2978202/1‑148 (MQ=255)
                                                                               cATCCCAAACGTATGTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                    <  2:2261115/148‑1 (MQ=255)
                                                                                aTCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTcc                                                                   >  1:498990/1‑148 (MQ=255)
                                                                                 tCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCt                                                                  >  2:3418903/1‑148 (MQ=255)
                                                                                  cccAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTg                                                                 <  2:2542788/148‑1 (MQ=255)
                                                                                      aaCGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATc                                                             >  2:319155/1‑148 (MQ=255)
                                                                                      aaCGGATTTGCTGTGCGTGTAAATTATCTGCTTTATCGGGCACAAACAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATc                                                             <  2:2144632/148‑1 (MQ=255)
                                                                                       aCGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCa                                                            >  1:292134/1‑148 (MQ=255)
                                                                                       aCGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCa                                                            >  2:2780110/1‑148 (MQ=255)
                                                                                       aCGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCa                                                            >  2:864669/1‑148 (MQ=255)
                                                                                        cGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCaa                                                           >  2:2279685/1‑148 (MQ=255)
                                                                                         gTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAg                                                          <  2:2042360/148‑1 (MQ=255)
                                                                                         gTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAg                                                          <  2:1366194/148‑1 (MQ=255)
                                                                                         gTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAg                                                          <  2:1475134/148‑1 (MQ=255)
                                                                                           aTTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCg                                                        >  2:1573055/1‑148 (MQ=255)
                                                                                            tttGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCgg                                                       >  1:2272481/1‑148 (MQ=255)
                                                                                             ttGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGt                                                      >  2:213268/1‑148 (MQ=255)
                                                                                                cTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAAt                                                   >  1:1872021/1‑148 (MQ=255)
                                                                                                 tgtgCGGGTAAATTATCTGCTTTATCGGGCACAACAAGAGCGCTCTGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATc                                                  <  1:2210753/148‑1 (MQ=255)
                                                                                                       ggTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTcc                                            >  2:3116845/1‑148 (MQ=255)
                                                                                                            aTTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATaa                                       >  1:818325/1‑148 (MQ=255)
                                                                                                             ttATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAc                                      >  1:3373420/1‑148 (MQ=255)
                                                                                                             ttATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAc                                      <  1:414350/148‑1 (MQ=255)
                                                                                                               aTCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCt                                    >  1:260111/1‑148 (MQ=255)
                                                                                                                tCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTg                                   <  1:2536996/148‑1 (MQ=255)
                                                                                                                tCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTg                                   >  2:799365/1‑148 (MQ=255)
                                                                                                                   tcTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAgt                                >  1:1164247/2‑148 (MQ=255)
                                                                                                                    cTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAgtg                               >  2:1237150/1‑148 (MQ=255)
                                                                                                                    cTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAgtg                               <  1:78614/148‑1 (MQ=255)
                                                                                                                      ttATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTa                             >  2:1448290/1‑148 (MQ=255)
                                                                                                                        aTCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACt                           <  2:1088942/148‑1 (MQ=255)
                                                                                                                          cGGGCACAAAAAGAGCGCGATGACACATTCGATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCAGAACCTGAGTGTACTGt                         >  2:3412462/1‑148 (MQ=255)
                                                                                                                            ggCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTc                       >  2:2244866/1‑148 (MQ=255)
                                                                                                                              cacaAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtg                     <  1:3440801/148‑1 (MQ=255)
                                                                                                                                 aaaaaGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatg                  <  1:2061986/148‑1 (MQ=255)
                                                                                                                                  aaaaGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatga                 >  2:2855051/1‑148 (MQ=255)
                                                                                                                                  aaaaGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatga                 >  2:1344844/1‑148 (MQ=255)
                                                                                                                                  aaaaGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatga                 <  1:2444961/148‑1 (MQ=255)
                                                                                                                                  aaaaGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatga                 <  1:525774/148‑1 (MQ=255)
                                                                                                                                  aaaaGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatga                 >  2:758706/1‑148 (MQ=255)
                                                                                                                                    aaGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAAt               <  2:3168721/148‑1 (MQ=255)
                                                                                                                                      gagCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCg             >  2:2678811/1‑148 (MQ=255)
                                                                                                                                      gagCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCg             >  1:1096975/1‑148 (MQ=255)
                                                                                                                                       agCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGa            >  1:1533033/1‑148 (MQ=255)
                                                                                                                                       agCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGa            >  2:2272487/1‑148 (MQ=255)
                                                                                                                                        gcgcGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGaa           <  1:2571118/148‑1 (MQ=255)
                                                                                                                                           cgATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAttt        >  2:3261627/1‑148 (MQ=255)
                                                                                                                                            gATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTa       <  2:3634596/148‑1 (MQ=255)
                                                                                                                                             aTGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACACCTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTAt      <  2:2520419/148‑1 (MQ=255)
                                                                                                                                              tGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATg     >  2:2578112/1‑148 (MQ=255)
                                                                                                                                               gAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGc    >  1:2314251/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGT  >  CP009273/2663993‑2664285

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
CTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGT  >  CP009273/2663989‑2664287
                                                                                                                                                     |                                                                                                                                                     
CTGGAAACGGGCCTGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAC                                                                                                                                                       >  SRR19973069.3828900/1‑150 (MQ=60)
CTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAC                                                                                                                                                       <  SRR19973069.3056230/150‑1 (MQ=60)
CTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAC                                                                                                                                                       <  SRR19973069.3056483/150‑1 (MQ=60)
   GAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACCTCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACA                                                                                                                                                    <  SRR19973069.2981401/150‑1 (MQ=60)
    AAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACAT                                                                                                                                                   <  SRR19973069.1870809/150‑1 (MQ=60)
        GGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTA                                                                                                                                               <  SRR19973069.2587960/150‑1 (MQ=60)
        GGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTA                                                                                                                                               <  SRR19973069.2588015/150‑1 (MQ=60)
          GCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATG                                                                                                                                             <  SRR19973069.512277/150‑1 (MQ=60)
          GCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATG                                                                                                                                             <  SRR19973069.3622657/150‑1 (MQ=60)
          GCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATG                                                                                                                                             <  SRR19973069.3622465/150‑1 (MQ=60)
            ATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTT                                                                                                                                           <  SRR19973069.2705350/150‑1 (MQ=60)
               GCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGA                                                                                                                                        <  SRR19973069.2844185/150‑1 (MQ=60)
                 CTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAAC                                                                                                                                      <  SRR19973069.768000/150‑1 (MQ=60)
                  TGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACG                                                                                                                                     >  SRR19973069.3421233/1‑150 (MQ=60)
                      CAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTT                                                                                                                                 <  SRR19973069.3429160/150‑1 (MQ=60)
                         TAGAGCCGGCGTCACGCACCCGTCAATTAATCCGCAACTGACAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGAGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCG                                                                                                                              <  SRR19973069.3761872/150‑1 (MQ=60)
                         AAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCG                                                                                                                              <  SRR19973069.363084/150‑1 (MQ=60)
                         AAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCG                                                                                                                              <  SRR19973069.363131/150‑1 (MQ=60)
                         AAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCG                                                                                                                              <  SRR19973069.3762558/150‑1 (MQ=60)
                         AAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCG                                                                                                                              >  SRR19973069.2559612/1‑150 (MQ=60)
                          AGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGA                                                                                                                             >  SRR19973069.3412398/1‑150 (MQ=60)
                            AGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACA                                                                                                                           >  SRR19973069.3500528/1‑150 (MQ=60)
                            AGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACA                                                                                                                           >  SRR19973069.1149680/1‑150 (MQ=60)
                             GCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACAGTCTATGTTGGAACGCGTTTCGACAA                                                                                                                          >  SRR19973069.99766/1‑150 (MQ=60)
                                 CCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTT                                                                                                                      <  SRR19973069.1606954/150‑1 (MQ=60)
                                   GTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCG                                                                                                                    >  SRR19973069.729639/1‑150 (MQ=60)
                                     CACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCC                                                                                                                  >  SRR19973069.2986629/1‑150 (MQ=60)
                                     CACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCC                                                                                                                  >  SRR19973069.2084263/1‑150 (MQ=60)
                                     CACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCC                                                                                                                  <  SRR19973069.2598820/150‑1 (MQ=60)
                                         CACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCT                                                                                                              >  SRR19973069.2723778/1‑150 (MQ=60)
                                           CCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGG                                                                                                            >  SRR19973069.2537179/1‑150 (MQ=60)
                                           CCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGG                                                                                                            >  SRR19973069.2537195/1‑150 (MQ=60)
                                              GTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAG                                                                                                         <  SRR19973069.964192/150‑1 (MQ=60)
                                              GTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAG                                                                                                         <  SRR19973069.490753/150‑1 (MQ=60)
                                                 ACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATG                                                                                                      <  SRR19973069.723922/150‑1 (MQ=60)
                                                 ACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATG                                                                                                      <  SRR19973069.725666/150‑1 (MQ=60)
                                                 ACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATG                                                                                                      <  SRR19973069.2710115/150‑1 (MQ=60)
                                                  CTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGA                                                                                                     <  SRR19973069.3441465/150‑1 (MQ=60)
                                                  CTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGA                                                                                                     >  SRR19973069.2043785/1‑150 (MQ=60)
                                                   TTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGAC                                                                                                    <  SRR19973069.3699062/150‑1 (MQ=60)
                                                     AATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAA                                                                                                  >  SRR19973069.929214/1‑150 (MQ=60)
                                                      ATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAAC                                                                                                 >  SRR19973069.1284570/1‑150 (MQ=60)
                                                       TCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACT                                                                                                <  SRR19973069.1909900/150‑1 (MQ=60)
                                                         AGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGG                                                                                              >  SRR19973069.1305635/1‑150 (MQ=60)
                                                          GCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGC                                                                                             >  SRR19973069.3486557/1‑150 (MQ=60)
                                                          GCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGC                                                                                             >  SRR19973069.3486584/1‑150 (MQ=60)
                                                           CAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCG                                                                                            >  SRR19973069.797374/1‑150 (MQ=60)
                                                            AACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAATGTTCGCCGTCTGGAAGATGACAACTGGCGC                                                                                           >  SRR19973069.3032237/1‑150 (MQ=60)
                                                            AACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGC                                                                                           >  SRR19973069.2896974/1‑150 (MQ=60)
                                                             ACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCT                                                                                          <  SRR19973069.662674/150‑1 (MQ=60)
                                                              CTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTT                                                                                         >  SRR19973069.2060657/1‑150 (MQ=60)
                                                                   AGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGG                                                                                    >  SRR19973069.3248981/1‑150 (MQ=60)
                                                                    GATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGA                                                                                   >  SRR19973069.348518/1‑150 (MQ=60)
                                                                    GATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGA                                                                                   >  SRR19973069.3206138/1‑150 (MQ=60)
                                                                     ATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAA                                                                                  >  SRR19973069.2359692/1‑150 (MQ=60)
                                                                        AGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGG                                                                               >  SRR19973069.3476081/1‑150 (MQ=60)
                                                                        AGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGG                                                                               >  SRR19973069.3477234/1‑150 (MQ=60)
                                                                        AGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGG                                                                               >  SRR19973069.2992306/1‑150 (MQ=60)
                                                                          CGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGA                                                                             <  SRR19973069.3239103/150‑1 (MQ=60)
                                                                           GAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGAT                                                                            >  SRR19973069.803652/1‑150 (MQ=60)
                                                                           GAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGAT                                                                            >  SRR19973069.1953364/1‑150 (MQ=60)
                                                                            AAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATA                                                                           <  SRR19973069.1677825/150‑1 (MQ=60)
                                                                            AAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATA                                                                           <  SRR19973069.364036/150‑1 (MQ=60)
                                                                            AAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATA                                                                           <  SRR19973069.361972/150‑1 (MQ=60)
                                                                            AAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATA                                                                           <  SRR19973069.2809601/150‑1 (MQ=60)
                                                                               CCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCG                                                                        <  SRR19973069.2331108/150‑1 (MQ=60)
                                                                               CCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCG                                                                        <  SRR19973069.2331144/150‑1 (MQ=60)
                                                                                 GACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTC                                                                      <  SRR19973069.2321300/150‑1 (MQ=60)
                                                                                 GACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTC                                                                      >  SRR19973069.2997976/1‑150 (MQ=60)
                                                                                  ACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCC                                                                     >  SRR19973069.502571/1‑150 (MQ=60)
                                                                                   CATCCCAAACGTATTTGCTGTGCGTGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCT                                                                    <  SRR19973069.770075/150‑1 (MQ=60)
                                                                                   CATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCT                                                                    <  SRR19973069.770029/150‑1 (MQ=60)
                                                                                         AAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCA                                                              >  SRR19973069.294143/1‑150 (MQ=60)
                                                                                         AAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCA                                                              <  SRR19973069.956691/150‑1 (MQ=60)
                                                                                         AAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCA                                                              >  SRR19973069.3189367/1‑150 (MQ=60)
                                                                                           ACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAG                                                            <  SRR19973069.801520/150‑1 (MQ=60)
                                                                                             GTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCG                                                          >  SRR19973069.3195930/1‑150 (MQ=60)
                                                                                              TATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGG                                                         >  SRR19973069.2770017/1‑150 (MQ=60)
                                                                                              TATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGG                                                         >  SRR19973069.2288124/1‑150 (MQ=60)
                                                                                              TATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGG                                                         >  SRR19973069.2287898/1‑150 (MQ=60)
                                                                                                  TGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAAT                                                     >  SRR19973069.1885050/1‑150 (MQ=60)
                                                                                                  TGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAAT                                                     >  SRR19973069.3787581/1‑150 (MQ=60)
                                                                                                     TGTGCGGGTAAATTATCTGCTTTATCGGGCACAACAAGAGCGCTCTGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCAC                                                  <  SRR19973069.2225732/150‑1 (MQ=60)
                                                                                                     TGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCAC                                                  <  SRR19973069.2225799/150‑1 (MQ=60)
                                                                                                          GGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCC                                             <  SRR19973069.3764248/150‑1 (MQ=60)
                                                                                                            GTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCAT                                           <  SRR19973069.931016/150‑1 (MQ=60)
                                                                                                             TAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATA                                          >  SRR19973069.342381/1‑150 (MQ=60)
                                                                                                             TAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATA                                          >  SRR19973069.341930/1‑150 (MQ=60)
                                                                                                              AAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAA                                         >  SRR19973069.824106/1‑150 (MQ=60)
                                                                                                               AATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAC                                        >  SRR19973069.3394251/1‑150 (MQ=60)
                                                                                                               AATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAC                                        >  SRR19973069.3395477/1‑150 (MQ=60)
                                                                                                                 TTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCT                                      >  SRR19973069.261986/1‑150 (MQ=60)
                                                                                                                 TTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCT                                      <  SRR19973069.417391/150‑1 (MQ=60)
                                                                                                                    TCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAG                                   <  SRR19973069.2553801/150‑1 (MQ=60)
                                                                                                                     CTTCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGT                                  >  SRR19973069.1172336/1‑150 (MQ=60)
                                                                                                                      TGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTG                                 >  SRR19973069.894410/1‑150 (MQ=60)
                                                                                                                        CTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTA                               <  SRR19973069.79470/150‑1 (MQ=60)
                                                                                                                        CTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTA                               <  SRR19973069.79491/150‑1 (MQ=60)
                                                                                                                                 GCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGA                      <  SRR19973069.2028191/150‑1 (MQ=60)
                                                                                                                                  CACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGAT                     <  SRR19973069.3463136/150‑1 (MQ=60)
                                                                                                                                   ACACAAAGAGAGCTATGAAACATTCTATGTTGGATCGCTTTTAGAAAAAGTTCGCCGTCTGGAAGATGACTACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATAAATTCCCATAACCTGAGTGTACTGTTCTGATG                    <  SRR19973069.1889538/150‑1 (MQ=60)
                                                                                                                                    CAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGA                   >  SRR19973069.3021155/1‑150 (MQ=60)
                                                                                                                                     AAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAA                  <  SRR19973069.3453923/150‑1 (MQ=60)
                                                                                                                                     AAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAA                  <  SRR19973069.132638/150‑1 (MQ=60)
                                                                                                                                     AAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAA                  <  SRR19973069.2076282/150‑1 (MQ=60)
                                                                                                                                     AAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAA                  >  SRR19973069.2217669/1‑150 (MQ=60)
                                                                                                                                      AAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAAT                 <  SRR19973069.529495/150‑1 (MQ=60)
                                                                                                                                      AAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAAT                 <  SRR19973069.2461420/150‑1 (MQ=60)
                                                                                                                                      AAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAAT                 <  SRR19973069.1773697/150‑1 (MQ=60)
                                                                                                                                      AAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAAT                 <  SRR19973069.292837/150‑1 (MQ=60)
                                                                                                                                        AAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCG               >  SRR19973069.3557440/1‑150 (MQ=60)
                                                                                                                                        AAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCG               >  SRR19973069.1104641/1‑150 (MQ=60)
                                                                                                                                         AGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGA              >  SRR19973069.1543441/1‑150 (MQ=60)
                                                                                                                                         AGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGA              >  SRR19973069.1543565/1‑150 (MQ=60)
                                                                                                                                           AGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAT            <  SRR19973069.2829199/150‑1 (MQ=60)
                                                                                                                                            GCGCGCTGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATT           <  SRR19973069.2588142/150‑1 (MQ=60)
                                                                                                                                            GCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATT           <  SRR19973069.2588152/150‑1 (MQ=60)
                                                                                                                                               CGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTAT        >  SRR19973069.3491152/1‑150 (MQ=60)
                                                                                                                                                 ATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGC      >  SRR19973069.2330003/1‑150 (MQ=60)
                                                                                                                                                    ACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTG   <  SRR19973069.1111376/150‑1 (MQ=60)
                                                                                                                                                     AACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGT  >  SRR19973069.2532961/1‑150 (MQ=60)
                                                                                                                                                     AACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGT  >  SRR19973069.2532831/1‑150 (MQ=60)
                                                                                                                                                     |                                                                                                                                                     
CTGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGTGT  >  CP009273/2663989‑2664287

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: ‑