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 F361 I0 R1 1353 250.7 8230166 96.0% 7900959 148.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,664,138 A→C 58.7% 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,1380AC58.7% 212.2 / 225.1 206E130D (GAA→GAChcaF3‑phenylpropionate dioxygenase, small (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (39/46);  new base C (48/73);  total (87/119)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 211 total aligned reads displayed.

GAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGT  >  CP009273/2663992‑2664285
                                                                                                                                                  |                                                                                                                                                   
gAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGacac                                                                                                                                                   >  1:3783284/1‑149 (MQ=255)
  aaCGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACAt                                                                                                                                                 <  2:3913327/149‑1 (MQ=255)
      ggCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTa                                                                                                                                             >  2:2961877/1‑149 (MQ=255)
         aTGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAGATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGt                                                                                                                                          <  2:3184426/149‑1 (MQ=255)
         aTGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGt                                                                                                                                          >  2:1111139/1‑149 (MQ=255)
         aTGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGt                                                                                                                                          >  2:3163758/1‑149 (MQ=255)
          tGGTCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGtt                                                                                                                                         <  1:492624/149‑1 (MQ=255)
          tGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCACACGGATGTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGtt                                                                                                                                         <  2:3215631/149‑1 (MQ=255)
          tGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGtt                                                                                                                                         <  2:2107944/149‑1 (MQ=255)
          tGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGtt                                                                                                                                         <  1:945556/149‑1 (MQ=255)
          tGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGtt                                                                                                                                         >  2:186559/1‑149 (MQ=255)
              cTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGaa                                                                                                                                     >  1:3062326/1‑149 (MQ=255)
              cTGGGCAGAAGAGCCGCCGTCACGCAACAGTCACTTAATCAGCAACTGCCAGATAAGCGACACCGACATAACAAACGTATTTGATGTGCGGGAAAATAAACTGCTTTATCGGGCCCAAAAAGAGCGCGATGACAGCGTCTCTGTTGGaa                                                                                                                                     >  2:261659/1‑149 (MQ=255)
               tGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAAc                                                                                                                                    >  2:2409099/1‑149 (MQ=255)
                gggCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAAcg                                                                                                                                   >  1:1447055/1‑149 (MQ=255)
                gggCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAAcg                                                                                                                                   >  2:156555/1‑149 (MQ=255)
                    agaagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGtt                                                                                                                               >  1:1335287/1‑149 (MQ=255)
                      aagaGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTc                                                                                                                             <  1:3571787/149‑1 (MQ=255)
                          gccgccGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACa                                                                                                                         >  1:1712027/1‑149 (MQ=255)
                            cgccgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACaaa                                                                                                                       <  2:553977/149‑1 (MQ=255)
                             gccgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAg                                                                                                                      >  1:1477648/1‑149 (MQ=255)
                             gccgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAg                                                                                                                      >  1:1044997/1‑149 (MQ=255)
                               cgTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGtt                                                                                                                    <  1:786951/149‑1 (MQ=255)
                                gTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTc                                                                                                                   <  2:2168652/149‑1 (MQ=255)
                                 tCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCg                                                                                                                  >  1:560099/1‑149 (MQ=255)
                                  cACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGc                                                                                                                 <  1:2940456/149‑1 (MQ=255)
                                  cACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGc                                                                                                                 <  2:1244594/149‑1 (MQ=255)
                                       aCCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCt                                                                                                            >  1:2564980/1‑149 (MQ=255)
                                             cACTTAATCAGCAACTGCCAGATAATCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGa                                                                                                      >  1:3219987/1‑149 (MQ=255)
                                             cACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGa                                                                                                      >  1:1164420/1‑149 (MQ=255)
                                                       gCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTgg                                                                                            >  1:1764334/1‑149 (MQ=255)
                                                        cAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGgc                                                                                           <  2:45980/149‑1 (MQ=255)
                                                           cTGCCAGATAAGCGAAACCGACATCCCACACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCt                                                                                        <  2:2822310/149‑1 (MQ=255)
                                                           cTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCt                                                                                        <  1:2266399/149‑1 (MQ=255)
                                                           cTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCt                                                                                        <  1:1978925/149‑1 (MQ=255)
                                                            tGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCtt                                                                                       <  2:2868413/149‑1 (MQ=255)
                                                             gCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTg                                                                                      >  2:1367987/1‑149 (MQ=255)
                                                              ccAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGc                                                                                     <  2:1638003/149‑1 (MQ=255)
                                                                  aTAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGa                                                                                 >  2:3147293/1‑149 (MQ=255)
                                                                      gCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACgg                                                                             <  2:807991/149‑1 (MQ=255)
                                                                      gCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACgg                                                                             >  1:2046350/1‑149 (MQ=255)
                                                                      gCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACgg                                                                             <  1:606395/149‑1 (MQ=255)
                                                                        gAAACCGACATCCCACACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGa                                                                           <  1:3144449/149‑1 (MQ=255)
                                                                        gAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGa                                                                           <  1:1857764/149‑1 (MQ=255)
                                                                          aaCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGata                                                                         <  1:53130/149‑1 (MQ=255)
                                                                           aCCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGatat                                                                        <  2:3770991/149‑1 (MQ=255)
                                                                              gACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGt                                                                     <  2:2639440/149‑1 (MQ=255)
                                                                               aCATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                    >  1:1652124/1‑149 (MQ=255)
                                                                               aCATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                    <  2:2371845/149‑1 (MQ=255)
                                                                               aCATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                    <  1:605208/149‑1 (MQ=255)
                                                                               aCATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTc                                                                    >  2:3066653/1‑149 (MQ=255)
                                                                                cATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTcc                                                                   <  1:4010237/149‑1 (MQ=255)
                                                                                cATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTcc                                                                   <  1:3906595/149‑1 (MQ=255)
                                                                                  tCCCAAACTTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTg                                                                 <  2:3296361/149‑1 (MQ=255)
                                                                                     cAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGAt                                                              <  1:2389458/149‑1 (MQ=255)
                                                                                      aaaCGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATc                                                             <  2:360793/149‑1 (MQ=255)
                                                                                          gTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGc                                                         <  1:683745/149‑1 (MQ=255)
                                                                                          gTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGc                                                         >  2:2990616/1‑149 (MQ=255)
                                                                                          gTATTTGCTGTGCGGGTAAATTAACTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGAACAAGc                                                         <  2:3696903/149‑1 (MQ=255)
                                                                                            aTTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCgg                                                       <  2:4108382/149‑1 (MQ=255)
                                                                                              ttGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAGAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTa                                                     >  2:3359344/1‑149 (MQ=255)
                                                                                               tGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTaa                                                    <  1:1928177/149‑1 (MQ=255)
                                                                                               tGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTaa                                                    <  2:3410862/149‑1 (MQ=255)
                                                                                                 cTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATc                                                  >  1:1421662/1‑149 (MQ=255)
                                                                                                  tgtgCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCa                                                 <  2:3486168/149‑1 (MQ=255)
                                                                                                   gtgCGGTTAATTTATCTTCTTTATCGGTCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGAGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCAc                                                <  2:2609333/149‑1 (MQ=255)
                                                                                                   gtgCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCAc                                                <  2:2609379/149‑1 (MQ=255)
                                                                                                         gTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCa                                          >  1:2151731/1‑149 (MQ=255)
                                                                                                          tAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCAt                                         <  2:1193014/149‑1 (MQ=255)
                                                                                                            aaTTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATaa                                       <  1:2680523/149‑1 (MQ=255)
                                                                                                            aaTTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATaa                                       <  2:154212/149‑1 (MQ=255)
                                                                                                                aTCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTg                                   >  2:2015765/1‑149 (MQ=255)
                                                                                                                  cTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAg                                 <  2:560099/149‑1 (MQ=255)
                                                                                                                    gCTTTATCGGGCACAAAACGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAgtg                               <  2:3973127/149‑1 (MQ=255)
                                                                                                                    gCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAgtg                               <  2:2112934/149‑1 (MQ=255)
                                                                                                                     cTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAgtgt                              <  1:2628981/149‑1 (MQ=255)
                                                                                                                      tttATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTa                             <  2:1907064/149‑1 (MQ=255)
                                                                                                                        tataGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACt                           >  2:3439797/1‑149 (MQ=255)
                                                                                                                         atagggCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTg                          <  1:2711982/146‑1 (MQ=255)
                                                                                                                         aTCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTg                          >  2:2378667/1‑149 (MQ=255)
                                                                                                                          tCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGt                         >  1:705618/1‑149 (MQ=255)
                                                                                                                             ggCACAAACAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCt                      <  2:1973138/149‑1 (MQ=255)
                                                                                                                                 caAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatg                  >  2:2623496/1‑149 (MQ=255)
                                                                                                                                 caAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatg                  >  1:4052756/1‑149 (MQ=255)
                                                                                                                                  aaaaaGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCtgatga                 <  1:2961877/149‑1 (MQ=255)
                                                                                                                                    aaaGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAAt               <  2:1426191/149‑1 (MQ=255)
                                                                                                                                     aaGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATc              >  1:3046421/1‑149 (MQ=255)
                                                                                                                                      agagCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCg             >  1:138356/1‑149 (MQ=255)
                                                                                                                                      agagCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCg             >  1:2212956/1‑149 (MQ=255)
                                                                                                                                       gagCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGa            >  1:319074/1‑149 (MQ=255)
                                                                                                                                         gcgcGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAt          <  2:2390761/149‑1 (MQ=255)
                                                                                                                                          cgcgATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAtt         <  2:1501452/149‑1 (MQ=255)
                                                                                                                                          cgcgATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAtt         <  2:4065421/149‑1 (MQ=255)
                                                                                                                                           gcgATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAAttt        >  2:2230095/1‑149 (MQ=255)
                                                                                                                                             gATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAGTTTAt      <  1:3359344/149‑1 (MQ=255)
                                                                                                                                             gATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTAt      <  1:1367987/149‑1 (MQ=255)
                                                                                                                                              aTGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATg     >  1:581024/1‑149 (MQ=255)
                                                                                                                                               tGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGc    <  2:3635965/149‑1 (MQ=255)
                                                                                                                                               tGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGc    >  2:1692778/1‑149 (MQ=255)
                                                                                                                                                 aaaCATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCgt  <  2:4051624/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGT  >  CP009273/2663992‑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
TGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGT  >  CP009273/2663990‑2664285
                                                                                                                                                    |                                                                                                                                                   
TGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACA                                                                                                                                                    >  SRR19973068.3849455/1‑150 (MQ=60)
 GGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACAC                                                                                                                                                   >  SRR19973068.3805865/1‑150 (MQ=60)
  TAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATAGGGCACAAAAAGAGCGCGATGAAACA                                                                                                                                                  <  SRR19973068.1546071/150‑1 (MQ=60)
  GAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACA                                                                                                                                                  <  SRR19973068.88194/150‑1 (MQ=60)
           ATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTT                                                                                                                                         <  SRR19973068.2242336/150‑1 (MQ=60)
            TGGTCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTG                                                                                                                                        <  SRR19973068.496072/150‑1 (MQ=60)
            TGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTG                                                                                                                                        <  SRR19973068.951890/150‑1 (MQ=60)
               CCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAA                                                                                                                                     >  SRR19973068.1347534/1‑150 (MQ=60)
               CCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAA                                                                                                                                     >  SRR19973068.3081058/1‑150 (MQ=60)
                 TGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACG                                                                                                                                   >  SRR19973068.1456546/1‑150 (MQ=60)
                  GGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCACACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGC                                                                                                                                  <  SRR19973068.2902035/150‑1 (MQ=60)
                     CAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTT                                                                                                                               >  SRR19973068.1343925/1‑150 (MQ=60)
                      AGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTT                                                                                                                              >  SRR19973068.957580/1‑150 (MQ=60)
                        AAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCG                                                                                                                            <  SRR19973068.4024189/150‑1 (MQ=60)
                        AAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCG                                                                                                                            <  SRR19973068.3593946/150‑1 (MQ=60)
                        AAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCG                                                                                                                            <  SRR19973068.3593480/150‑1 (MQ=60)
                          GAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGAC                                                                                                                          >  SRR19973068.1388684/1‑150 (MQ=60)
                           AGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACA                                                                                                                         >  SRR19973068.1722787/1‑150 (MQ=60)
                           AGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACA                                                                                                                         >  SRR19973068.1722943/1‑150 (MQ=60)
                              CGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAG                                                                                                                      >  SRR19973068.1051722/1‑150 (MQ=60)
                              CGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAG                                                                                                                      >  SRR19973068.1487362/1‑150 (MQ=60)
                                 CGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTC                                                                                                                   <  SRR19973068.792308/150‑1 (MQ=60)
                                  GTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCG                                                                                                                  >  SRR19973068.562720/1‑150 (MQ=60)
                                  GTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCG                                                                                                                  >  SRR19973068.563956/1‑150 (MQ=60)
                                  GTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCG                                                                                                                  >  SRR19973068.2641431/1‑150 (MQ=60)
                                    CACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCC                                                                                                                <  SRR19973068.2958711/150‑1 (MQ=60)
                                        CACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCT                                                                                                            >  SRR19973068.2581229/1‑150 (MQ=60)
                                             GTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAG                                                                                                       >  SRR19973068.2529252/1‑150 (MQ=60)
                                              TCACTTAATCAGCAACTGCCAGATAATCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGA                                                                                                      >  SRR19973068.3239783/1‑150 (MQ=60)
                                              TCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGA                                                                                                      >  SRR19973068.1172201/1‑150 (MQ=60)
                                              TCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGA                                                                                                      >  SRR19973068.1442181/1‑150 (MQ=60)
                                              TCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGA                                                                                                      >  SRR19973068.3235930/1‑150 (MQ=60)
                                               CACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGAT                                                                                                     >  SRR19973068.252965/1‑150 (MQ=60)
                                                  TTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGAC                                                                                                  <  SRR19973068.2972746/150‑1 (MQ=60)
                                                        AGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGG                                                                                            >  SRR19973068.1775601/1‑150 (MQ=60)
                                                        AGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGG                                                                                            >  SRR19973068.1775726/1‑150 (MQ=60)
                                                        AGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGG                                                                                            >  SRR19973068.2699049/1‑150 (MQ=60)
                                                             CTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTT                                                                                       <  SRR19973068.1992017/150‑1 (MQ=60)
                                                             CTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTT                                                                                       <  SRR19973068.2281240/150‑1 (MQ=60)
                                                             CTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTT                                                                                       <  SRR19973068.2281110/150‑1 (MQ=60)
                                                             CTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTT                                                                                       <  SRR19973068.1992038/150‑1 (MQ=60)
                                                                CCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCT                                                                                    >  SRR19973068.520559/1‑150 (MQ=60)
                                                                CCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCT                                                                                    >  SRR19973068.520458/1‑150 (MQ=60)
                                                                      AAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACG                                                                              >  SRR19973068.701205/1‑150 (MQ=60)
                                                                       AGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGG                                                                             >  SRR19973068.2059698/1‑150 (MQ=60)
                                                                        GCGGAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGG                                                                            <  SRR19973068.610543/150‑1 (MQ=60)
                                                                        GCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGG                                                                            <  SRR19973068.610457/150‑1 (MQ=60)
                                                                          GAAACCGACATCCCACACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGAT                                                                          <  SRR19973068.3163754/150‑1 (MQ=60)
                                                                          GAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGAT                                                                          <  SRR19973068.1870045/150‑1 (MQ=60)
                                                                          GAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGAT                                                                          <  SRR19973068.3914618/150‑1 (MQ=60)
                                                                            AACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATAT                                                                        <  SRR19973068.53766/150‑1 (MQ=60)
                                                                               CGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGT                                                                     >  SRR19973068.2981036/1‑150 (MQ=60)
                                                                                GACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTC                                                                    >  SRR19973068.1662838/1‑150 (MQ=60)
                                                                                 ACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCC                                                                   <  SRR19973068.609218/150‑1 (MQ=60)
                                                                                 ACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCC                                                                   <  SRR19973068.609279/150‑1 (MQ=60)
                                                                                 ACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCC                                                                   <  SRR19973068.609267/150‑1 (MQ=60)
                                                                                  CATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCT                                                                  <  SRR19973068.4034458/150‑1 (MQ=60)
                                                                                  CATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCT                                                                  <  SRR19973068.3930229/150‑1 (MQ=60)
                                                                                       CAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATC                                                             <  SRR19973068.2404688/150‑1 (MQ=60)
                                                                                        AAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACAGTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCA                                                            >  SRR19973068.3192379/1‑150 (MQ=60)
                                                                                         AACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAA                                                           >  SRR19973068.1740420/1‑150 (MQ=60)
                                                                                            GTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCG                                                        <  SRR19973068.688076/150‑1 (MQ=60)
                                                                                               TTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTA                                                     >  SRR19973068.2227203/1‑150 (MQ=60)
                                                                                                 TGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAAT                                                   <  SRR19973068.1240820/150‑1 (MQ=60)
                                                                                                 TGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAAT                                                   <  SRR19973068.1940994/150‑1 (MQ=60)
                                                                                                  GCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATC                                                  >  SRR19973068.1430916/1‑150 (MQ=60)
                                                                                                      TGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTT                                              <  SRR19973068.578083/150‑1 (MQ=60)
                                                                                                       GCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTC                                             <  SRR19973068.1837982/150‑1 (MQ=60)
                                                                                                          GGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCA                                          >  SRR19973068.2165829/1‑150 (MQ=60)
                                                                                                          GGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCA                                          <  SRR19973068.3256689/150‑1 (MQ=60)
                                                                                                           GTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCAT                                         >  SRR19973068.4038480/1‑150 (MQ=60)
                                                                                                           GTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCAT                                         <  SRR19973068.2663941/150‑1 (MQ=60)
                                                                                                           GTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCAT                                         <  SRR19973068.1421281/150‑1 (MQ=60)
                                                                                                             AAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAA                                       <  SRR19973068.1971353/150‑1 (MQ=60)
                                                                                                              AATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAAC                                      <  SRR19973068.2697226/150‑1 (MQ=60)
                                                                                                                    CTGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGT                                >  SRR19973068.3080423/1‑150 (MQ=60)
                                                                                                                     TGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTG                               <  SRR19973068.3468164/150‑1 (MQ=60)
                                                                                                                     TGCTTTATCGGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTG                               >  SRR19973068.315609/1‑150 (MQ=60)
                                                                                                                      GCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGT                              >  SRR19973068.749560/1‑150 (MQ=60)
                                                                                                                       CTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTA                             <  SRR19973068.2645508/150‑1 (MQ=60)
                                                                                                                           ATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGT                         >  SRR19973068.710142/1‑150 (MQ=60)
                                                                                                                           ATAGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGT                         <  SRR19973068.2728785/150‑1 (MQ=60)
                                                                                                                              GGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCT                      <  SRR19973068.1482891/150‑1 (MQ=60)
                                                                                                                              GGGCACAAAAAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCT                      <  SRR19973068.2163436/150‑1 (MQ=60)
                                                                                                                                  ACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATG                  >  SRR19973068.4077142/1‑150 (MQ=60)
                                                                                                                                  ACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATG                  >  SRR19973068.4077342/1‑150 (MQ=60)
                                                                                                                                    AAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAA                <  SRR19973068.2980245/150‑1 (MQ=60)
                                                                                                                                      AAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATC              >  SRR19973068.3065092/1‑150 (MQ=60)
                                                                                                                                       AAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCG             >  SRR19973068.139592/1‑150 (MQ=60)
                                                                                                                                       AAGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCG             >  SRR19973068.2227419/1‑150 (MQ=60)
                                                                                                                                        AGAGCGCGATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGA            >  SRR19973068.321144/1‑150 (MQ=60)
                                                                                                                                         GAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAA           <  SRR19973068.188032/150‑1 (MQ=60)
                                                                                                                                               GATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATG     >  SRR19973068.584993/1‑150 (MQ=60)
                                                                                                                                               GATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATG     <  SRR19973068.3183213/150‑1 (MQ=60)
                                                                                                                                               GATGACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATG     <  SRR19973068.1118556/150‑1 (MQ=60)
                                                                                                                                               GATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAGTTTATG     <  SRR19973068.3379776/150‑1 (MQ=60)
                                                                                                                                               GATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATG     <  SRR19973068.1376735/150‑1 (MQ=60)
                                                                                                                                                  GACACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGT  <  SRR19973068.2724927/150‑1 (MQ=60)
                                                                                                                                                    |                                                                                                                                                   
TGGAAACGGGCATGGCCTGGGCAGAAGAGCCGCCGTCACGCACCCGTCACTTAATCAGCAACTGCCAGATAAGCGAAACCGACATCCCAAACGTATTTGCTGTGCGGGTAAATTATCTGCTTTATCGGGCACAAAAAGAGCGCGATGAAACATTCTATGTTGGAACGCGTTTCGACAAAGTTCGCCGTCTGGAAGATGACAACTGGCGCTTGCTGGAACGGGATATCGTCCTGGATCAAGCGGTAATCACTTCCCATAACCTGAGTGTACTGTTCTGATGAATCGAATTTATGCGT  >  CP009273/2663990‑2664285

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: