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 F394 I0 R1 1342 215.0 7134094 96.2% 6862998 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,461,041 A→T 95.0% pseudogene (1393/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,461,0410AT95.0% 546.2 / 13.7 178pseudogene (1393/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base A (5/4);  new base T (104/65);  total (109/69)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 179 total aligned reads displayed.

CGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCA  >  CP009273/1460896‑1461188
                                                                                                                                                 |                                                                                                                                                   
cGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTcc                                                                                                                                                   >  2:530159/1‑148 (MQ=255)
      tgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGAc                                                                                                                                             >  1:3310691/1‑148 (MQ=255)
       gaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACa                                                                                                                                            <  2:2615747/148‑1 (MQ=255)
          gCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAgcg                                                                                                                                         <  1:2522986/148‑1 (MQ=255)
          gCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAgcg                                                                                                                                         <  2:1885679/148‑1 (MQ=255)
           ccACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAgcgc                                                                                                                                        >  2:1498920/1‑148 (MQ=255)
            cACGGCAAACAACATCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGATATTGCTGGAATTAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCa                                                                                                                                       <  1:2667217/148‑1 (MQ=14)
            cACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCa                                                                                                                                       >  1:1866634/1‑148 (MQ=255)
             aCGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCaa                                                                                                                                      <  1:2625861/148‑1 (MQ=255)
                gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACgg                                                                                                                                   <  2:508026/148‑1 (MQ=255)
                gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACgg                                                                                                                                   <  2:1866634/148‑1 (MQ=255)
                 cAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGt                                                                                                                                  <  2:481954/148‑1 (MQ=255)
                    acaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGa                                                                                                                               <  2:3481932/148‑1 (MQ=255)
                    acaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGa                                                                                                                               >  1:3481932/1‑148 (MQ=255)
                    acaacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGa                                                                                                                               >  1:1306298/1‑148 (MQ=255)
                      aacaacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATa                                                                                                                             <  2:1562329/148‑1 (MQ=255)
                        caacGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAAc                                                                                                                           >  2:2549847/1‑148 (MQ=255)
                         acccgCAACACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACa                                                                                                                          <  2:2719196/144‑1 (MQ=255)
                          acGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACaa                                                                                                                         >  2:3004141/1‑148 (MQ=255)
                           cGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAg                                                                                                                        >  2:316342/1‑148 (MQ=255)
                           cGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAg                                                                                                                        >  2:2864991/1‑148 (MQ=255)
                              cAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGCCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATGGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGc                                                                                                                     <  1:3435720/148‑1 (MQ=255)
                                aaaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCAc                                                                                                                   <  2:1255718/148‑1 (MQ=255)
                                 aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCAcc                                                                                                                  >  2:660157/1‑148 (MQ=255)
                                 aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCAcc                                                                                                                  <  1:2870678/148‑1 (MQ=255)
                                 aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCAcc                                                                                                                  <  1:277878/148‑1 (MQ=255)
                                 aaCCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGAGAACAAGGGCAcc                                                                                                                  >  1:1622604/1‑148 (MQ=255)
                                  accaccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCa                                                                                                                 >  2:1715480/1‑148 (MQ=255)
                                   ccaccGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAt                                                                                                                >  2:2971804/1‑148 (MQ=255)
                                     accGTTGACGGCAAGGATTCCACCGGTACGGAACTTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGa                                                                                                              <  2:3230987/148‑1 (MQ=255)
                                         ttGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGt                                                                                                          >  2:3211754/1‑148 (MQ=255)
                                          tGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTc                                                                                                         >  2:2784573/1‑148 (MQ=255)
                                          tGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTc                                                                                                         >  2:2062102/1‑148 (MQ=255)
                                           gACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCa                                                                                                        <  1:292752/148‑1 (MQ=255)
                                           gACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCa                                                                                                        <  1:3155036/148‑1 (MQ=255)
                                           gACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCa                                                                                                        <  1:3155046/148‑1 (MQ=255)
                                            aCGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAc                                                                                                       >  2:3493131/1‑148 (MQ=255)
                                                cAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                                                                                   >  1:2701961/1‑148 (MQ=255)
                                                cAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                                                                                   >  1:1723220/1‑148 (MQ=255)
                                                  aGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATcc                                                                                                 >  1:241057/1‑148 (MQ=255)
                                                   ggATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                >  1:1877680/1‑148 (MQ=255)
                                                    gATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                               >  1:42236/1‑148 (MQ=255)
                                                    gATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                               >  1:1247149/1‑148 (MQ=255)
                                                    gATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                               >  1:1247142/1‑148 (MQ=255)
                                                    gATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                               >  1:1845584/1‑148 (MQ=255)
                                                    gATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                               <  1:2084198/148‑1 (MQ=255)
                                                     aTTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                                                                              <  1:2615325/148‑1 (MQ=255)
                                                       tCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                                                                            >  1:79413/1‑148 (MQ=255)
                                                          aCCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCa                                                                                         <  2:228991/148‑1 (MQ=255)
                                                               tACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGt                                                                                    >  2:3310018/1‑148 (MQ=255)
                                                               tACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGt                                                                                    >  2:3309962/1‑148 (MQ=255)
                                                                aCGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTa                                                                                   >  1:1680189/1‑148 (MQ=255)
                                                                 cGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAt                                                                                  <  2:421868/148‑1 (MQ=255)
                                                                 cGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAt                                                                                  >  1:3158359/1‑148 (MQ=255)
                                                                  ggAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATc                                                                                 >  2:1584517/1‑148 (MQ=255)
                                                                        tGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGAtt                                                                           >  2:2915847/1‑148 (MQ=255)
                                                                         gCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTg                                                                          <  2:1362999/148‑1 (MQ=255)
                                                                           tGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGAc                                                                        <  1:2784573/148‑1 (MQ=255)
                                                                             gCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACgg                                                                      >  2:483720/1‑148 (MQ=255)
                                                                                      ggAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGc                                                             >  1:2878913/1‑148 (MQ=38)
                                                                                      ggAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGc                                                             >  1:886090/1‑148 (MQ=38)
                                                                                         aGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCAt                                                          <  1:3004141/148‑1 (MQ=38)
                                                                                         aGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCAt                                                          <  1:1490989/148‑1 (MQ=38)
                                                                                          ggTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATc                                                         >  1:1912886/1‑148 (MQ=38)
                                                                                          ggTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATc                                                         >  1:1912747/1‑148 (MQ=38)
                                                                                           gTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCg                                                        <  2:2277902/148‑1 (MQ=38)
                                                                                           gTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAt                                                                                                                >  1:3391882/1‑92 (MQ=33)
                                                                                           gTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAt                                                                                                                <  2:3391882/92‑1 (MQ=33)
                                                                                             gATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTc                                                      >  1:1153160/1‑148 (MQ=38)
                                                                                             gATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTc                                                      >  1:1152048/1‑148 (MQ=38)
                                                                                             gATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTc                                                      >  1:10533/1‑148 (MQ=38)
                                                                                                tCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAAt                                                   >  2:13534/1‑148 (MQ=255)
                                                                                                  aGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATaa                                                 >  1:1501373/1‑148 (MQ=255)
                                                                                                  aGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATaa                                                 >  1:3272713/1‑148 (MQ=255)
                                                                                                    gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACg                                               >  2:205884/1‑148 (MQ=255)
                                                                                                    gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACg                                               >  2:2006979/1‑148 (MQ=255)
                                                                                                    gACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACg                                               >  1:1979533/1‑148 (MQ=255)
                                                                                                       ggCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAg                                            <  1:2062102/148‑1 (MQ=255)
                                                                                                         cGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGc                                          >  2:220425/1‑148 (MQ=255)
                                                                                                           aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACCGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCga                                        >  2:1311923/1‑148 (MQ=255)
                                                                                                           aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGaa                                             <  2:2517594/143‑1 (MQ=255)
                                                                                                           aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGaa                                             >  1:2517594/1‑143 (MQ=255)
                                                                                                                gATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCACTAACGAAGGCGAGAGCa                                   >  1:677100/1‑148 (MQ=255)
                                                                                                                gATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCa                                   >  2:2445240/1‑148 (MQ=255)
                                                                                                                gATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCa                                   >  2:2444354/1‑148 (MQ=255)
                                                                                                                 aTGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAgg                                                              <  2:1744631/120‑1 (MQ=17)
                                                                                                                 aTGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAc                                  >  1:2494238/1‑148 (MQ=255)
                                                                                                                 aTGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAc                                  >  1:668780/1‑148 (MQ=255)
                                                                                                                       gcggcggcgGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAc                            >  2:459833/1‑148 (MQ=255)
                                                                                                                           cggcggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAAt                        >  1:1025986/1‑148 (MQ=255)
                                                                                                                            ggcggTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATg                       <  2:103275/148‑1 (MQ=255)
                                                                                                                                 tCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcggc                  >  2:3271285/1‑148 (MQ=255)
                                                                                                                                 tCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcggc                  >  2:1499628/1‑148 (MQ=255)
                                                                                                                                 tCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATggcggc                  >  1:1614699/1‑148 (MQ=255)
                                                                                                                                       tATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggc            <  1:2553516/148‑1 (MQ=255)
                                                                                                                                        aTTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcaccggca           <  2:668780/148‑1 (MQ=255)
                                                                                                                                          tGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACt         >  2:3117226/1‑148 (MQ=255)
                                                                                                                                           gATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTc        <  2:1247142/148‑1 (MQ=255)
                                                                                                                                              atCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGa     >  1:2818510/1‑148 (MQ=255)
                                                                                                                                                 aCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGCTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCa  <  2:149271/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCA  >  CP009273/1460896‑1461188

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
TAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATC  >  CP009273/1460893‑1461187
                                                                                                                                                    |                                                                                                                                                  
TAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCTCCGTTGACGGCCAGTATTCCACCGGTACGGTAATTTCTGGCAATAACGGTAATGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTC                                                                                                                                                   <  SRR19973066.2698727/150‑1 (MQ=51)
       GATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGAC                                                                                                                                            >  SRR19973066.100318/1‑150 (MQ=60)
       GATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGAC                                                                                                                                            >  SRR19973066.3332608/1‑150 (MQ=60)
       GATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGAC                                                                                                                                            >  SRR19973066.3332653/1‑150 (MQ=60)
             GCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCA                                                                                                                                      >  SRR19973066.1879549/1‑150 (MQ=60)
             GCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCA                                                                                                                                      >  SRR19973066.1879967/1‑150 (MQ=60)
             GCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCA                                                                                                                                      <  SRR19973066.2540146/150‑1 (MQ=60)
             GCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCA                                                                                                                                      <  SRR19973066.2541276/150‑1 (MQ=60)
               CACGGCAAACAACATCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGATATTGCTGGAATTAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAAC                                                                                                                                    <  SRR19973066.2685061/150‑1 (MQ=60)
                ACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACG                                                                                                                                   <  SRR19973066.2643564/150‑1 (MQ=60)
                   GCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTG                                                                                                                                <  SRR19973066.3283552/150‑1 (MQ=60)
                    CTAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGG                                                                                                                               <  SRR19973066.343907/150‑1 (MQ=60)
                    CAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATGGCGGGCAATAACGGTAATGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGG                                                                                                                               <  SRR19973066.343465/150‑1 (MQ=38)
                     AAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGA                                                                                                                              >  SRR19973066.3505162/1‑150 (MQ=60)
                     AAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGA                                                                                                                              >  SRR19973066.1315771/1‑150 (MQ=60)
                          CCAACGGCAACACCACCGTTGACGGCAAGGAGTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTGTGTAGGGCGTTCTGTATGTCAGCGGCGTCGTTCAAGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACA                                                                                                                         <  SRR19973066.3331249/150‑1 (MQ=60)
                               GGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGC                                                                                                                    <  SRR19973066.1707127/150‑1 (MQ=60)
                               GGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGC                                                                                                                    <  SRR19973066.3081627/150‑1 (MQ=60)
                                GCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCA                                                                                                                   <  SRR19973066.1663622/150‑1 (MQ=60)
                                 CAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGCCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATGGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCAC                                                                                                                  <  SRR19973066.3458709/150‑1 (MQ=60)
                                 CAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCAC                                                                                                                  <  SRR19973066.3066300/150‑1 (MQ=60)
                                  AAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACC                                                                                                                 >  SRR19973066.283898/1‑150 (MQ=60)
                                  AAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGAGAACAAGGGCACC                                                                                                                 >  SRR19973066.1634364/1‑150 (MQ=60)
                                   AAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCA                                                                                                                >  SRR19973066.1914615/1‑150 (MQ=60)
                                    AACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT                                                                                                               <  SRR19973066.279905/150‑1 (MQ=60)
                                    AACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT                                                                                                               <  SRR19973066.1097029/150‑1 (MQ=60)
                                    AACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT                                                                                                               <  SRR19973066.280010/150‑1 (MQ=60)
                                    AACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCAT                                                                                                               <  SRR19973066.2889930/150‑1 (MQ=60)
                                        ccccgagacggcaaggcatccaccggtaagcccagatcgcatttgtagtgtctggagttcagccgtgtggtcttcccttcCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACC                                                                                                           <  SRR19973066.533911/70‑1 (MQ=0)
                                           GTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTC                                                                                                        >  SRR19973066.2474237/1‑150 (MQ=60)
                                            TTGACGGCAAGGAGTCCACCGGTACGGAAATTGCTGGCAATACCGGGAAGGTGATTCCGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCA                                                                                                       <  SRR19973066.1596914/150‑1 (MQ=60)
                                              GACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACC                                                                                                     <  SRR19973066.294958/150‑1 (MQ=60)
                                              GACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACC                                                                                                     <  SRR19973066.294845/150‑1 (MQ=60)
                                              GACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACC                                                                                                     <  SRR19973066.3176069/150‑1 (MQ=60)
                                              GACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACC                                                                                                     <  SRR19973066.3176078/150‑1 (MQ=60)
                                              GACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACC                                                                                                     <  SRR19973066.3176079/150‑1 (MQ=60)
                                                 GGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAT                                                                                                  >  SRR19973066.1735692/1‑150 (MQ=60)
                                                 GGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGAT                                                                                                  >  SRR19973066.2720053/1‑150 (MQ=60)
                                                  gcaaccccaccaccgcgacccaaatcgctggcaatcacgggaaggtgattcaggagcgctttgtgcttgactggcgttgaggcgtgtgctctgactaaACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATC                                                                                                 <  SRR19973066.2134837/52‑1 (MQ=0)
                                                   CAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCC                                                                                                >  SRR19973066.242944/1‑150 (MQ=60)
                                                   CAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCC                                                                                                >  SRR19973066.242774/1‑150 (MQ=60)
                                                    AAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCG                                                                                               >  SRR19973066.1891088/1‑150 (MQ=60)
                                                     AGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGG                                                                                              >  SRR19973066.1858785/1‑150 (MQ=60)
                                                     AGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGG                                                                                              >  SRR19973066.1256089/1‑150 (MQ=60)
                                                     AGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGG                                                                                              >  SRR19973066.1256096/1‑150 (MQ=60)
                                                     AGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGG                                                                                              >  SRR19973066.42798/1‑150 (MQ=60)
                                                     AGGATTCCACCGGTACGGAAATGGCTGGCCATTACGGGTAGGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGG                                                                                              <  SRR19973066.1226153/150‑1 (MQ=29)
                                                       GATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAG                                                                                            <  SRR19973066.2098575/150‑1 (MQ=60)
                                                        ATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGT                                                                                           <  SRR19973066.2632997/150‑1 (MQ=60)
                                                        ATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGT                                                                                           >  SRR19973066.80329/1‑150 (MQ=60)
                                                        ATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGT                                                                                           <  SRR19973066.2632972/150‑1 (MQ=60)
                                                            gagcgcgcattgtttgtgtctggagatctgtcgtgtgatggtccgcgcgGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAT                                                                                       <  SRR19973066.1320933/101‑1 (MQ=0)
                                                            gagatcgcattggttgtgactgtagttcagacgtgtgctcttccgatgtGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAT                                                                                       <  SRR19973066.1320891/101‑1 (MQ=0)
                                                                 GTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTA                                                                                  >  SRR19973066.1692417/1‑150 (MQ=60)
                                                                  TACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAT                                                                                 >  SRR19973066.3179406/1‑150 (MQ=60)
                                                                   ACTGAAATTGCTCGCAATATCGGGAAGGTGATTCAGGGCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACTGTGGATAAAAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATC                                                                                <  SRR19973066.639095/150‑1 (MQ=31)
                                                                     GGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCA                                                                              <  SRR19973066.2523352/150‑1 (MQ=59)
                                                                      gacaggcCTGGCAATACCGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGCTATCGCCGGCGACAGCGCACCGGTGGATTACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAG                                                                             <  SRR19973066.3022086/143‑1 (MQ=29)
                                                                      GAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAG                                                                             >  SRR19973066.274345/1‑150 (MQ=57)
                                                                       AAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGA                                                                            <  SRR19973066.3334985/150‑1 (MQ=57)
                                                                        AATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGAT                                                                           >  SRR19973066.1562340/1‑150 (MQ=57)
                                                                         ATTGCTGGCCATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTTTTGATATCTCCGTCTCCAGCGCAAAGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATT                                                                          <  SRR19973066.2342532/150‑1 (MQ=39)
                                                                              TGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGG                                                                     <  SRR19973066.2803321/150‑1 (MQ=50)
                                                                                GCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCG                                                                   <  SRR19973066.3483532/150‑1 (MQ=46)
                                                                                  AATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGAC                                                                 >  SRR19973066.2918580/1‑150 (MQ=42)
                                                                                     AACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAG                                                              >  SRR19973066.679156/1‑150 (MQ=36)
                                                                                       CGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGC                                                            >  SRR19973066.2898202/1‑150 (MQ=32)
                                                                                       CGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGC                                                            >  SRR19973066.892303/1‑150 (MQ=32)
                                                                                       CGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGC                                                            >  SRR19973066.892385/1‑150 (MQ=32)
                                                                                          GAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCAT                                                         <  SRR19973066.198163/150‑1 (MQ=32)
                                                                                           AAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATC                                                        >  SRR19973066.1926368/1‑150 (MQ=32)
                                                                                           AAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATC                                                        >  SRR19973066.1186728/1‑150 (MQ=32)
                                                                                           AAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATC                                                        >  SRR19973066.1926367/1‑150 (MQ=32)
                                                                                           AAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATC                                                        >  SRR19973066.1926509/1‑150 (MQ=32)
                                                                                            AGGTGTTTCAGGAGGGCGATGTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCG                                                       <  SRR19973066.3024119/150‑1 (MQ=21)
                                                                                            AGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCG                                                       >  SRR19973066.1691255/1‑150 (MQ=32)
                                                                                            AGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCG                                                       <  SRR19973066.1501421/150‑1 (MQ=32)
                                                                                              GTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTC                                                     >  SRR19973066.10705/1‑150 (MQ=32)
                                                                                              GTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTC                                                     >  SRR19973066.1161353/1‑150 (MQ=32)
                                                                                              GTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTC                                                     >  SRR19973066.1161372/1‑150 (MQ=32)
                                                                                              GTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTC                                                     >  SRR19973066.1160256/1‑150 (MQ=32)
                                                                                                 ATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTagatcggaagagcacacgtctcaactccagtcactacacatgcgatctcgtgttggccggggggcaggcggtggggggc                                                  >  SRR19973066.2247683/1‑71 (MQ=0)
                                                                                                   TCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAA                                                >  SRR19973066.1511832/1‑150 (MQ=36)
                                                                                                   TCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAA                                                >  SRR19973066.3294189/1‑150 (MQ=36)
                                                                                                     AGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACG                                              >  SRR19973066.1993512/1‑150 (MQ=40)
                                                                                                     AGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACG                                              >  SRR19973066.2662137/1‑150 (MQ=40)
                                                                                                          GGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGC                                         <  SRR19973066.2076394/150‑1 (MQ=50)
                                                                                                            CGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAagatcggaagagcacacgtctgaactccagtcactaccaatgc                                       >  SRR19973066.1068276/1‑107 (MQ=0)
                                                                                                            CGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGAGA                                       >  SRR19973066.2534716/1‑150 (MQ=49)
                                                                                                                 TGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCACTAACGAAGGCGAGAGCA                                  >  SRR19973066.682057/1‑150 (MQ=53)
                                                                                                                  GGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAC                                 >  SRR19973066.2511241/1‑150 (MQ=60)
                                                                                                                  GGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAC                                 >  SRR19973066.672547/1‑150 (MQ=60)
                                                                                                                  GGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCAC                                 >  SRR19973066.673664/1‑150 (MQ=60)
                                                                                                                  GGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCggagcggacgagcacacctctgaaca                                 >  SRR19973066.1757192/1‑124 (MQ=11)
                                                                                                                    ATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTA                               >  SRR19973066.1278082/1‑150 (MQ=60)
                                                                                                                       TCAGCGGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCAGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCA                            <  SRR19973066.37978/150‑1 (MQ=60)
                                                                                                                            GGCGGCGGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAAT                       >  SRR19973066.1033539/1‑150 (MQ=60)
                                                                                                                                  GGTCACGGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGC                 >  SRR19973066.1626413/1‑150 (MQ=60)
                                                                                                                                   GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCagatcggaagagcacacgtctgaactccagtcactcccaatgcgagctcgttggcgagagcgggctctccaacggggggg                >  SRR19973066.1672441/1‑70 (MQ=0)
                                                                                                                                        GGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGC           >  SRR19973066.13804/1‑150 (MQ=60)
                                                                                                                                        GGTATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGC           >  SRR19973066.2016263/1‑150 (MQ=60)
                                                                                                                                        GGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGTTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGC           >  SRR19973066.150674/1‑150 (MQ=60)
                                                                                                                                          TATTGATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCAC         <  SRR19973066.2570857/150‑1 (MQ=60)
                                                                                                                                           ATTGGTATCTCCGGCGACAGCGCAACGGTGGATAAGAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACT        <  SRR19973066.2634583/150‑1 (MQ=60)
                                                                                                                                               ATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGA    >  SRR19973066.2840537/1‑150 (MQ=60)
                                                                                                                                               ATATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGA    >  SRR19973066.2837453/1‑150 (MQ=60)
                                                                                                                                                 ATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATC  >  SRR19973066.1175107/1‑150 (MQ=60)
                                                                                                                                                 ATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATC  >  SRR19973066.1174245/1‑150 (MQ=60)
                                                                                                                                                 ATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATC  >  SRR19973066.1174241/1‑150 (MQ=60)
                                                                                                                                                 ATCTCCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATC  >  SRR19973066.1636729/1‑150 (MQ=60)
                                                                                                                                                    |                                                                                                                                                  
TAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATC  >  CP009273/1460893‑1461187

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