Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F434 I0 R1 1130 215.9 7031662 97.1% 6827743 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,639,116 C→A 100% A523A (GCG→GCT)  pbpC ← penicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,639,1160CA99.3% 559.2 / ‑3.9 186A523A (GCG→GCT) pbpCpenicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (102/81);  total (103/83)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 186 total aligned reads displayed.

CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGC  >  CP009273/2638972‑2639260
                                                                                                                                                |                                                                                                                                                
  gATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCa                                                                                                                                             >  1:133911/1‑148 (MQ=255)
  gATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCa                                                                                                                                             >  1:134006/1‑148 (MQ=255)
  gATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCa                                                                                                                                             <  2:3069665/148‑1 (MQ=255)
   aTTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAg                                                                                                                                            >  2:27283/1‑148 (MQ=255)
         aCTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCa                                                                                                                                      <  1:2443979/148‑1 (MQ=255)
          cTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAc                                                                                                                                     >  2:1233922/1‑148 (MQ=255)
               cgACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAg                                                                                                                                >  2:3448660/1‑148 (MQ=255)
               cgACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAg                                                                                                                                <  1:1771281/148‑1 (MQ=255)
                gACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGc                                                                                                                               <  2:1567612/148‑1 (MQ=255)
                  cAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCa                                                                                                                             >  1:3154309/1‑148 (MQ=255)
                   aGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCAt                                                                                                                            <  2:3434214/148‑1 (MQ=255)
                    gTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATa                                                                                                                           >  2:241720/1‑148 (MQ=255)
                    gTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATa                                                                                                                           >  1:719005/1‑148 (MQ=255)
                     tAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAg                                                                                                                          >  1:2266563/1‑148 (MQ=255)
                     tAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAg                                                                                                                          >  1:594757/1‑148 (MQ=255)
                     tAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAg                                                                                                                          <  2:1634032/148‑1 (MQ=255)
                      aaGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGc                                                                                                                         >  1:2877241/1‑148 (MQ=255)
                      aaGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGc                                                                                                                         >  1:435979/1‑148 (MQ=255)
                       aGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCt                                                                                                                        <  2:2574279/148‑1 (MQ=255)
                              attGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGccc                                                                                                                 >  1:3210624/1‑148 (MQ=255)
                                    cTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTc                                                                                                           >  1:1268911/1‑148 (MQ=255)
                                     tGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTcc                                                                                                          <  1:1990148/148‑1 (MQ=255)
                                       aTTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCAt                                                                                                        >  1:1733534/1‑148 (MQ=255)
                                        ttCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATg                                                                                                       >  1:3366604/1‑148 (MQ=255)
                                         tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                      >  1:3500805/1‑148 (MQ=255)
                                         tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                      >  1:944231/1‑148 (MQ=255)
                                         tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                      >  2:2647305/1‑148 (MQ=255)
                                         tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                      >  2:2747501/1‑148 (MQ=255)
                                           aacaaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCa                                                                                                    <  1:243889/148‑1 (MQ=255)
                                            acaaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAg                                                                                                   >  2:1492599/1‑148 (MQ=255)
                                               aTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTgg                                                                                                >  1:455609/1‑148 (MQ=255)
                                               aTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTgg                                                                                                <  2:349823/148‑1 (MQ=255)
                                               aTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTgg                                                                                                <  2:350668/148‑1 (MQ=255)
                                                tGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTggg                                                                                               >  1:2877103/1‑148 (MQ=255)
                                                tGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTggg                                                                                               >  1:827907/1‑148 (MQ=255)
                                                  gTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGc                                                                                             >  2:887954/1‑148 (MQ=255)
                                                     cGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAc                                                                                          >  2:1534187/1‑148 (MQ=255)
                                                        cACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcgcg                                                                                       >  1:107597/1‑148 (MQ=255)
                                                        cACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcgcg                                                                                       <  2:3015576/148‑1 (MQ=255)
                                                        cACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcgcg                                                                                       >  1:1883035/1‑148 (MQ=255)
                                                          cTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcgcgcg                                                                                     >  2:118136/1‑148 (MQ=255)
                                                             gCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCa                                                                                  <  1:525842/148‑1 (MQ=255)
                                                              cAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCaa                                                                                 <  2:199920/148‑1 (MQ=255)
                                                              cAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCaa                                                                                 >  2:1267901/1‑148 (MQ=255)
                                                              cAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCaa                                                                                 <  2:594757/148‑1 (MQ=255)
                                                                   ccAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAc                                                                            >  1:1057940/1‑148 (MQ=255)
                                                                    cAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACt                                                                           <  1:871074/148‑1 (MQ=255)
                                                                      aaCTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTAt                                                                         <  2:3238683/148‑1 (MQ=255)
                                                                      aaCTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTAt                                                                         <  2:197898/148‑1 (MQ=255)
                                                                       aCTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATc                                                                        <  2:1604109/148‑1 (MQ=255)
                                                                            ccAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCa                                                                   <  1:1582911/148‑1 (MQ=255)
                                                                              aacaacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAAGAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAg                                                                 <  2:2083389/148‑1 (MQ=255)
                                                                               acaacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAgg                                                                >  1:3180192/1‑148 (MQ=255)
                                                                               acaacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAgg                                                                >  1:2241794/1‑148 (MQ=255)
                                                                                caacGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAggg                                                               <  1:2610463/148‑1 (MQ=255)
                                                                                   cGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTg                                                            <  2:2226015/148‑1 (MQ=255)
                                                                                   cGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTg                                                            <  2:2226084/148‑1 (MQ=255)
                                                                                   cGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTg                                                            >  1:1502282/1‑148 (MQ=255)
                                                                                    gggCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgc                                                           <  2:2444719/148‑1 (MQ=255)
                                                                                     ggCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcg                                                          >  2:3140561/1‑148 (MQ=255)
                                                                                     ggCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcg                                                          <  2:2951790/148‑1 (MQ=255)
                                                                                      gCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcgc                                                         >  1:194520/1‑148 (MQ=255)
                                                                                          gCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtca                                                     >  1:1758033/1‑148 (MQ=255)
                                                                                              tCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAg                                                 >  1:1969699/1‑148 (MQ=255)
                                                                                               ccGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGc                                                >  2:942818/1‑148 (MQ=255)
                                                                                               ccGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGc                                                >  1:803622/1‑148 (MQ=255)
                                                                                                cGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGcc                                               <  2:1268911/148‑1 (MQ=255)
                                                                                                    cTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAAGAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATaa                                           <  1:117951/148‑1 (MQ=255)
                                                                                                      gCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGGCTCTGGGGCGTGGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATc                                         <  1:3448218/148‑1 (MQ=255)
                                                                                                      gCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATc                                         >  2:756664/1‑148 (MQ=255)
                                                                                                      gCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATc                                         <  1:3490171/148‑1 (MQ=255)
                                                                                                        cAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCACTAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCg                                       <  2:2528599/148‑1 (MQ=255)
                                                                                                         aGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCgg                                      <  1:2472897/148‑1 (MQ=255)
                                                                                                          gTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGc                                     >  1:1669540/1‑148 (MQ=255)
                                                                                                          gTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGc                                     >  1:209162/1‑148 (MQ=255)
                                                                                                            ccAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGa                                   <  2:1909918/148‑1 (MQ=255)
                                                                                                            ccAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGa                                   >  1:2742956/1‑148 (MQ=255)
                                                                                                            ccAGATCCCAATGACAGAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGa                                   <  1:457837/148‑1 (MQ=255)
                                                                                                                aTCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGa                               >  1:573180/1‑148 (MQ=255)
                                                                                                                 tCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGAt                              <  1:3030995/148‑1 (MQ=255)
                                                                                                                   ccAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATcc                            <  2:3118059/148‑1 (MQ=255)
                                                                                                                    cAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCa                           >  1:233501/1‑148 (MQ=255)
                                                                                                                      aTGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcg                         <  1:27283/148‑1 (MQ=255)
                                                                                                                      aTGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcg                         <  1:442074/148‑1 (MQ=255)
                                                                                                                        gACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcgcc                       >  1:49602/1‑148 (MQ=255)
                                                                                                                         aCATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcgccg                      >  1:2908180/1‑148 (MQ=255)
                                                                                                                             agcgagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGcccc                  >  2:2516036/1‑148 (MQ=255)
                                                                                                                                gagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcga               <  1:3238553/148‑1 (MQ=255)
                                                                                                                                gagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcga               <  1:2142740/148‑1 (MQ=255)
                                                                                                                                gagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcga               <  1:3238557/148‑1 (MQ=255)
                                                                                                                                  gcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcgacg             <  1:999005/148‑1 (MQ=255)
                                                                                                                                    gTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCgacgac           <  2:134006/148‑1 (MQ=255)
                                                                                                                                    gTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCgacgac           <  2:435979/148‑1 (MQ=255)
                                                                                                                                    gTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCgacgac           <  2:133911/148‑1 (MQ=255)
                                                                                                                                      tAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAt         >  2:3177138/1‑148 (MQ=255)
                                                                                                                                      tAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAt         >  2:3483741/1‑148 (MQ=255)
                                                                                                                                      tAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAt         >  2:3037442/1‑148 (MQ=255)
                                                                                                                                      tAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAt         >  2:1747324/1‑148 (MQ=255)
                                                                                                                                            cAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAgg   <  2:719005/148‑1 (MQ=255)
                                                                                                                                            cAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAgg   >  1:2837098/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGC  >  CP009273/2638972‑2639260

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
TCCGGCAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGTT  >  CP009273/2638967‑2639264
                                                                                                                                                     |                                                                                                                                                    
TCCGGCAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATA                                                                                                                                                      >  SRR19973065.3375269/1‑150 (MQ=60)
 CCGGCAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAG                                                                                                                                                     <  SRR19973065.901015/150‑1 (MQ=60)
     CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCA                                                                                                                                                 >  SRR19973065.135060/1‑150 (MQ=60)
     CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCA                                                                                                                                                 >  SRR19973065.135156/1‑150 (MQ=60)
     CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCA                                                                                                                                                 >  SRR19973065.135088/1‑150 (MQ=60)
     CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCA                                                                                                                                                 >  SRR19973065.2266249/1‑150 (MQ=60)
       GATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGG                                                                                                                                               <  SRR19973065.1288404/150‑1 (MQ=60)
        ATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGC                                                                                                                                              >  SRR19973065.1716148/1‑150 (MQ=60)
              ACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACG                                                                                                                                        <  SRR19973065.2459380/150‑1 (MQ=60)
                    CGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCC                                                                                                                                  <  SRR19973065.1782878/150‑1 (MQ=60)
                    CGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCC                                                                                                                                  <  SRR19973065.1782855/150‑1 (MQ=60)
                     GACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCA                                                                                                                                 >  SRR19973065.3173443/1‑150 (MQ=60)
                       CAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATA                                                                                                                               >  SRR19973065.1208263/1‑150 (MQ=60)
                       CAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATA                                                                                                                               >  SRR19973065.2969773/1‑150 (MQ=60)
                       CAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATA                                                                                                                               >  SRR19973065.723819/1‑150 (MQ=60)
                        AGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAG                                                                                                                              >  SRR19973065.598439/1‑150 (MQ=60)
                        AGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAG                                                                                                                              >  SRR19973065.2280886/1‑150 (MQ=60)
                         GTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGC                                                                                                                             >  SRR19973065.2894861/1‑150 (MQ=60)
                         GTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGC                                                                                                                             >  SRR19973065.439629/1‑150 (MQ=60)
                         GTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGC                                                                                                                             >  SRR19973065.439007/1‑150 (MQ=60)
                             GATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGT                                                                                                                         >  SRR19973065.864993/1‑150 (MQ=60)
                                 TTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCC                                                                                                                     >  SRR19973065.3229928/1‑150 (MQ=60)
                                 TTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCC                                                                                                                     >  SRR19973065.3229937/1‑150 (MQ=60)
                                 TTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCC                                                                                                                     >  SRR19973065.1944419/1‑150 (MQ=60)
                                     TGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTT                                                                                                                 >  SRR19973065.1589926/1‑150 (MQ=60)
                                      GACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTT                                                                                                                <  SRR19973065.2320384/150‑1 (MQ=60)
                                      GACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTAAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTT                                                                                                                <  SRR19973065.794588/150‑1 (MQ=60)
                                       ACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTC                                                                                                               >  SRR19973065.1277285/1‑150 (MQ=60)
                                          TGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCAT                                                                                                            >  SRR19973065.1745042/1‑150 (MQ=60)
                                          TGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCAT                                                                                                            >  SRR19973065.3400092/1‑150 (MQ=60)
                                          TGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCAT                                                                                                            <  SRR19973065.2003057/150‑1 (MQ=60)
                                          TGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCAT                                                                                                            >  SRR19973065.1745008/1‑150 (MQ=60)
                                           GATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATG                                                                                                           >  SRR19973065.3387190/1‑150 (MQ=60)
                                            ATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGC                                                                                                          >  SRR19973065.3522066/1‑150 (MQ=60)
                                            ATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGC                                                                                                          >  SRR19973065.3522060/1‑150 (MQ=60)
                                            ATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGC                                                                                                          >  SRR19973065.950468/1‑150 (MQ=60)
                                                AACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGT                                                                                                      <  SRR19973065.245620/150‑1 (MQ=60)
                                                 ACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTG                                                                                                     >  SRR19973065.2181948/1‑150 (MQ=60)
                                                  CAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGG                                                                                                    >  SRR19973065.458759/1‑150 (MQ=60)
                                                   AATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGG                                                                                                   >  SRR19973065.833326/1‑150 (MQ=60)
                                                   AATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGG                                                                                                   >  SRR19973065.2894723/1‑150 (MQ=60)
                                                    ATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGG                                                                                                  <  SRR19973065.2091292/150‑1 (MQ=60)
                                                        TACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAC                                                                                              >  SRR19973065.1165784/1‑150 (MQ=60)
                                                        TACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAC                                                                                              >  SRR19973065.2738030/1‑150 (MQ=60)
                                                        TACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCAATACCAGATCCCAATGACATAGCGAGCGTTAACACCAATAGACCAGGCGTCTCGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAC                                                                                              >  SRR19973065.3256547/1‑150 (MQ=60)
                                                           GGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCG                                                                                           >  SRR19973065.1895431/1‑150 (MQ=60)
                                                           GGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCG                                                                                           >  SRR19973065.108577/1‑150 (MQ=60)
                                                             CACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCG                                                                                         <  SRR19973065.2909464/150‑1 (MQ=60)
                                                                  GCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAG                                                                                    <  SRR19973065.529294/150‑1 (MQ=60)
                                                                   CAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGG                                                                                   >  SRR19973065.2884607/1‑150 (MQ=60)
                                                                      AGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAC                                                                                >  SRR19973065.1064952/1‑150 (MQ=60)
                                                                        CCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTA                                                                              >  SRR19973065.989506/1‑150 (MQ=60)
                                                                         CAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGTCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTAT                                                                             >  SRR19973065.3261220/1‑150 (MQ=60)
                                                                         CAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTAT                                                                             <  SRR19973065.876650/150‑1 (MQ=60)
                                                                           AACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCC                                                                           <  SRR19973065.612274/150‑1 (MQ=60)
                                                                                 CCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAG                                                                     <  SRR19973065.1593469/150‑1 (MQ=60)
                                                                                  CAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGG                                                                    >  SRR19973065.3199409/1‑150 (MQ=60)
                                                                                  CAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGG                                                                    >  SRR19973065.3199412/1‑150 (MQ=60)
                                                                                  CAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGG                                                                    >  SRR19973065.1963162/1‑150 (MQ=60)
                                                                                  CAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGG                                                                    >  SRR19973065.2256037/1‑150 (MQ=60)
                                                                                     CAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTT                                                                 <  SRR19973065.2626648/150‑1 (MQ=60)
                                                                                      AACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTG                                                                >  SRR19973065.1511995/1‑150 (MQ=60)
                                                                                         GGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGC                                                             >  SRR19973065.196033/1‑150 (MQ=60)
                                                                                             GTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCA                                                         >  SRR19973065.897151/1‑150 (MQ=60)
                                                                                             GTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCA                                                         >  SRR19973065.1769585/1‑150 (MQ=60)
                                                                                                 CGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAG                                                     >  SRR19973065.1982543/1‑150 (MQ=60)
                                                                                                  GTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGC                                                    >  SRR19973065.808934/1‑150 (MQ=60)
                                                                                                    CCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCA                                                  <  SRR19973065.38768/150‑1 (MQ=60)
                                                                                                      GGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATA                                                >  SRR19973065.1783837/1‑150 (MQ=60)
                                                                                                         CTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAAGAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATC                                             <  SRR19973065.118986/150‑1 (MQ=60)
                                                                                                           GCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGGCTCTGGGGCGTGGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCG                                           <  SRR19973065.3469242/150‑1 (MQ=60)
                                                                                                           GCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCG                                           <  SRR19973065.3469689/150‑1 (MQ=60)
                                                                                                           GCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCG                                           <  SRR19973065.3511395/150‑1 (MQ=60)
                                                                                                             CAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGC                                         >  SRR19973065.1680707/1‑150 (MQ=60)
                                                                                                             CAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGC                                         >  SRR19973065.210745/1‑150 (MQ=60)
                                                                                                              AGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCG                                        <  SRR19973065.2488464/150‑1 (MQ=60)
                                                                                                              AGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCG                                        <  SRR19973065.2488756/150‑1 (MQ=60)
                                                                                                               GTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGA                                       >  SRR19973065.2760022/1‑150 (MQ=60)
                                                                                                                 CCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAAT                                     <  SRR19973065.460531/150‑1 (MQ=60)
                                                                                                                 CCAGATCCCAATGACAGAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAAT                                     <  SRR19973065.460997/150‑1 (MQ=60)
                                                                                                                   AGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGA                                   >  SRR19973065.577316/1‑150 (MQ=60)
                                                                                                                   AGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGA                                   >  SRR19973065.576792/1‑150 (MQ=60)
                                                                                                                   AGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGA                                   >  SRR19973065.576596/1‑150 (MQ=60)
                                                                                                                     ATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATC                                 >  SRR19973065.1425433/1‑150 (MQ=60)
                                                                                                                      TCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCC                                <  SRR19973065.3049562/150‑1 (MQ=60)
                                                                                                                       CCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCA                               >  SRR19973065.235193/1‑150 (MQ=60)
                                                                                                                       CCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCA                               >  SRR19973065.2140384/1‑150 (MQ=60)
                                                                                                                         CAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACG                             <  SRR19973065.3159651/150‑1 (MQ=60)
                                                                                                                           ATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCC                           >  SRR19973065.50165/1‑150 (MQ=60)
                                                                                                                           ATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCC                           >  SRR19973065.50274/1‑150 (MQ=60)
                                                                                                                           ATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCC                           <  SRR19973065.27648/150‑1 (MQ=60)
                                                                                                                           ATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCC                           <  SRR19973065.445137/150‑1 (MQ=60)
                                                                                                                           ATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCC                           <  SRR19973065.27643/150‑1 (MQ=60)
                                                                                                                            TGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCG                          >  SRR19973065.2925891/1‑150 (MQ=60)
                                                                                                                                  AGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTACCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATACGGCGAATGATCCACTCCGCCCCCG                    >  SRR19973065.489925/1‑150 (MQ=60)
                                                                                                                                     GAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACG                 <  SRR19973065.3258060/150‑1 (MQ=60)
                                                                                                                                     GAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACG                 <  SRR19973065.2156523/150‑1 (MQ=60)
                                                                                                                                     GAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACG                 <  SRR19973065.3258056/150‑1 (MQ=60)
                                                                                                                                      AGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGA                >  SRR19973065.2836056/1‑150 (MQ=60)
                                                                                                                                       GCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGAC               <  SRR19973065.1005643/150‑1 (MQ=60)
                                                                                                                                               CCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGG       >  SRR19973065.2854563/1‑150 (MQ=60)
                                                                                                                                                CCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGC      <  SRR19973065.241261/150‑1 (MQ=60)
                                                                                                                                                  AATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCG    <  SRR19973065.2469754/150‑1 (MQ=60)
                                                                                                                                                    TAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGTT  <  SRR19973065.893592/150‑1 (MQ=60)
                                                                                                                                                    TAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGTT  <  SRR19973065.893593/150‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                                    
TCCGGCAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGTT  >  CP009273/2638967‑2639264

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