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 3,217,289 A→T 100% H433L (CAC→CTC)  ebgA → evolved beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,217,2890AT99.5% 624.7 / ‑4.6 191H433L (CAC→CTC) ebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (84/106);  total (85/106)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 191 total aligned reads displayed.

AACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGA  >  CP009273/3217145‑3217433
                                                                                                                                                |                                                                                                                                                
aaCTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTcgc                                                                                                                                               >  1:3203829/1‑148 (MQ=255)
       tgATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGaaa                                                                                                                                        <  1:451744/148‑1 (MQ=255)
        gATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGaaaa                                                                                                                                       <  1:222259/148‑1 (MQ=255)
         atatCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGaaaaa                                                                                                                                      >  2:1026211/1‑148 (MQ=255)
         atatCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGaaaaa                                                                                                                                      >  2:2308010/1‑148 (MQ=255)
            tCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCa                                                                                                                                   >  1:1117362/1‑148 (MQ=255)
            tCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCa                                                                                                                                   >  1:909368/1‑148 (MQ=255)
             cTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCAt                                                                                                                                  <  2:2690099/148‑1 (MQ=255)
              tACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATc                                                                                                                                 <  2:1334469/148‑1 (MQ=255)
                   ccTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCg                                                                                                                            <  1:21169/148‑1 (MQ=255)
                       tttGTGATGGAGGAAACAGACGACGAATAGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCATTGGGAAAAATTCTACGTAGAGCGCAATGTACGCCATATCCTCGCGCAGAAAAACCATCAGACGaaca                                                                                                                        >  2:1843641/1‑148 (MQ=255)
                        ttGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcat                                                                                                                       >  1:65887/1‑148 (MQ=255)
                         tgtgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatc                                                                                                                      <  1:2454422/148‑1 (MQ=255)
                           tgATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatcat                                                                                                                    <  2:225708/148‑1 (MQ=255)
                            gATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGatcatcatc                                                                                                                   <  1:2120976/148‑1 (MQ=255)
                               ggCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTgg                                                                                                                >  2:2058296/1‑148 (MQ=255)
                                    aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCt                                                                                                           >  2:1248156/1‑148 (MQ=255)
                                    aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCt                                                                                                           <  2:2661742/148‑1 (MQ=255)
                                    aaaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCt                                                                                                           >  2:1983944/1‑148 (MQ=255)
                                     aaCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTg                                                                                                          <  1:402548/148‑1 (MQ=255)
                                      aCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTgg                                                                                                         >  2:3003972/1‑148 (MQ=255)
                                           cGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAATTCTACTTAGAGCGCATTGTTCGACATATACTAGCGCAGAAAAACCATCCGACGATCATCATCTGGTCGCTGGGCAAt                                                                                                    >  2:3302935/1‑148 (MQ=255)
                                               gAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAAt                                                                                                <  2:1927291/148‑1 (MQ=255)
                                                      aCGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATCTCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTa                                                                                         <  2:535281/148‑1 (MQ=255)
                                                        ggCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATg                                                                                       >  1:3341672/1‑148 (MQ=255)
                                                         gCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATgg                                                                                      <  2:825993/148‑1 (MQ=255)
                                                         gCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATgg                                                                                      >  1:96977/1‑148 (MQ=255)
                                                           tttGCTTATGTCGGCGATTTTAGCCGTATTACCGACGATCCTCAGTGGGAAAAAGTCTACGTCGAGCGCTTTGTTCGCCATATCCTCGCGCAGTAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCt                                                                                    <  1:1470103/148‑1 (MQ=255)
                                                             tGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGt                                                                                  >  1:2268046/1‑148 (MQ=255)
                                                                tAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAAc                                                                               <  1:1659513/148‑1 (MQ=255)
                                                                 aaTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACa                                                                              >  2:3217025/1‑148 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             <  2:3372133/148‑1 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             <  1:2925182/148‑1 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             <  1:2899011/148‑1 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             >  2:207309/1‑148 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             >  2:3331567/1‑148 (MQ=255)
                                                                  aTGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAt                                                                             <  1:1147286/148‑1 (MQ=255)
                                                                   tGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCCGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATc                                                                            >  1:2400915/1‑148 (MQ=255)
                                                                   tGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATc                                                                            >  2:3333882/1‑148 (MQ=255)
                                                                   tGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATc                                                                            >  1:461555/1‑148 (MQ=255)
                                                                       ggCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCcgcg                                                                        <  1:2722654/148‑1 (MQ=255)
                                                                       ggCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCcgcg                                                                        >  2:3257714/1‑148 (MQ=255)
                                                                          gATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGa                                                                     <  1:1002445/148‑1 (MQ=255)
                                                                           atatTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGAt                                                                    >  1:3499278/1‑148 (MQ=255)
                                                                            tatTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATg                                                                   <  2:1141078/148‑1 (MQ=255)
                                                                            tatTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATg                                                                   >  1:1614805/1‑148 (MQ=255)
                                                                             atTAGCCGTATTCCCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGt                                                                  <  2:2016578/148‑1 (MQ=255)
                                                                             atTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGt                                                                  >  1:1838623/1‑148 (MQ=255)
                                                                                aGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTAcc                                                               <  1:3065082/148‑1 (MQ=255)
                                                                                  ccGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCAt                                                             >  2:1750476/1‑148 (MQ=255)
                                                                                    gTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATgc                                                           >  2:1606789/1‑148 (MQ=255)
                                                                                          ccGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGaa                                                     <  2:3137470/148‑1 (MQ=255)
                                                                                          ccGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGaa                                                     <  2:1503295/148‑1 (MQ=255)
                                                                                           cgacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGaaa                                                    <  1:2857201/148‑1 (MQ=255)
                                                                                           cgacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGaaa                                                    <  1:2532958/148‑1 (MQ=255)
                                                                                           cgacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGaaa                                                    >  2:2523646/1‑148 (MQ=255)
                                                                                            gacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   <  2:2911099/148‑1 (MQ=255)
                                                                                            gacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   <  2:1161441/148‑1 (MQ=255)
                                                                                            gacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   <  2:1787491/148‑1 (MQ=255)
                                                                                            gacgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   <  1:936313/148‑1 (MQ=255)
                                                                                              cgaTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAgcg                                                 <  1:1361664/148‑1 (MQ=255)
                                                                                                aTCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCt                                               >  2:1837806/1‑148 (MQ=255)
                                                                                                  ccGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTgg                                             >  1:1499488/1‑148 (MQ=255)
                                                                                                  ccGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTgg                                             <  2:2464675/148‑1 (MQ=255)
                                                                                                   cGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGa                                            <  2:1747945/148‑1 (MQ=255)
                                                                                                   cGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGa                                            <  2:753592/148‑1 (MQ=255)
                                                                                                   cGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGa                                            <  2:1041028/148‑1 (MQ=255)
                                                                                                   cGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGa                                            <  1:1983944/148‑1 (MQ=255)
                                                                                                      aGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGa                                         <  2:2010169/148‑1 (MQ=255)
                                                                                                       gTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   <  1:626536/137‑1 (MQ=255)
                                                                                                       gTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAg                                                   >  2:626736/1‑137 (MQ=255)
                                                                                                        tGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGaca                                       >  1:2241452/1‑148 (MQ=255)
                                                                                                        tGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGaca                                       >  2:2777340/1‑148 (MQ=255)
                                                                                                        tGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGaca                                       <  2:2653080/148‑1 (MQ=255)
                                                                                                         gggAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGacac                                      >  1:1863368/1‑148 (MQ=255)
                                                                                                          ggAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACAcg                                     <  2:1840434/148‑1 (MQ=255)
                                                                                                           gAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACAcgc                                    >  1:2152149/1‑148 (MQ=255)
                                                                                                               aaGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACt                                >  1:2632253/1‑148 (MQ=255)
                                                                                                                aGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTg                               <  1:1671287/148‑1 (MQ=255)
                                                                                                                 gTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTgg                              <  2:3499278/148‑1 (MQ=255)
                                                                                                                      cGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCAt                         <  1:359868/148‑1 (MQ=255)
                                                                                                                        tCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGGGGCGAAAGCGCTGGTTGACACGCGACTGGTGCATTa                       <  1:1811304/148‑1 (MQ=255)
                                                                                                                        tCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTCGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTa                       >  2:828867/1‑148 (MQ=255)
                                                                                                                            gcgcATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACgaa                   >  2:463447/1‑148 (MQ=255)
                                                                                                                            gcgcATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACgaa                   <  1:3305855/148‑1 (MQ=255)
                                                                                                                             cgcATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACgaag                  <  2:912458/148‑1 (MQ=255)
                                                                                                                             cgcATTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACgaag                  <  1:2201260/148‑1 (MQ=255)
                                                                                                                                aTTGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGaagaag               <  2:362625/148‑1 (MQ=255)
                                                                                                                                 ttGTTCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAagaaga              <  2:350053/148‑1 (MQ=255)
                                                                                                                                     tCGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcgc          >  2:2470667/1‑148 (MQ=255)
                                                                                                                                      cGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcgcg         <  2:928856/148‑1 (MQ=255)
                                                                                                                                      cGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcgcg         <  2:3198349/148‑1 (MQ=255)
                                                                                                                                      cGCCATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATcgcg         <  1:3390641/148‑1 (MQ=255)
                                                                                                                                        ccATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGAt       >  1:2617742/1‑148 (MQ=255)
                                                                                                                                         cATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATg      >  2:1656499/1‑148 (MQ=255)
                                                                                                                                         cATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATg      >  1:3237347/1‑148 (MQ=255)
                                                                                                                                         cATATCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATg      <  1:1026211/148‑1 (MQ=255)
                                                                                                                                          atatCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGACAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGc     >  2:2867240/1‑148 (MQ=255)
                                                                                                                                          atatCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGc     >  2:824129/1‑148 (MQ=255)
                                                                                                                                             tCCTCGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGa  <  2:2752279/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAGAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGA  >  CP009273/3217145‑3217433

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
TTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGAAGTG  >  CP009273/3217140‑3217437
                                                                                                                                                     |                                                                                                                                                     
TTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCT                                                                                                                                                       >  SRR19973065.39589/1‑150 (MQ=60)
   CGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGC                                                                                                                                                    >  SRR19973065.3223114/1‑150 (MQ=60)
      ACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGTCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCA                                                                                                                                                 >  SRR19973065.1148455/1‑150 (MQ=60)
            TGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCTATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAA                                                                                                                                          <  SRR19973065.424611/150‑1 (MQ=60)
            TGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAA                                                                                                                                          <  SRR19973065.2233411/150‑1 (MQ=60)
            TGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAA                                                                                                                                          <  SRR19973065.454874/150‑1 (MQ=60)
             GATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAAC                                                                                                                                         <  SRR19973065.223913/150‑1 (MQ=60)
             GATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAAC                                                                                                                                         <  SRR19973065.223908/150‑1 (MQ=60)
              ATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCTCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCTAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACC                                                                                                                                        >  SRR19973065.227366/1‑150 (MQ=60)
               TATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCA                                                                                                                                       >  SRR19973065.1124655/1‑150 (MQ=60)
               TATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCA                                                                                                                                       >  SRR19973065.915267/1‑150 (MQ=60)
                ATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCAT                                                                                                                                      >  SRR19973065.1498059/1‑150 (MQ=60)
                     CGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTC                                                                                                                                 <  SRR19973065.39867/150‑1 (MQ=60)
                     CGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTC                                                                                                                                 <  SRR19973065.1829053/150‑1 (MQ=60)
                     CGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTC                                                                                                                                 >  SRR19973065.352503/1‑150 (MQ=60)
                        CCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGAT                                                                                                                              <  SRR19973065.21456/150‑1 (MQ=60)
                        CCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGAT                                                                                                                              <  SRR19973065.21416/150‑1 (MQ=60)
                        CCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGAT                                                                                                                              <  SRR19973065.593228/150‑1 (MQ=60)
                          TGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATACTCGCGCAG‑AAAAACCATCCGTCGATCA                                                                                                                            <  SRR19973065.8865/150‑1 (MQ=60)
                           GTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCAT                                                                                                                           >  SRR19973065.66591/1‑150 (MQ=60)
                            TTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATC                                                                                                                          >  SRR19973065.1415452/1‑150 (MQ=60)
                              TGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCAT                                                                                                                        <  SRR19973065.2469891/150‑1 (MQ=60)
                                 GATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTG                                                                                                                     <  SRR19973065.2134620/150‑1 (MQ=60)
                                 GATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTG                                                                                                                     <  SRR19973065.2134629/150‑1 (MQ=60)
                                 GATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTG                                                                                                                     <  SRR19973065.2135332/150‑1 (MQ=60)
                                   TGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGT                                                                                                                   <  SRR19973065.2940655/150‑1 (MQ=60)
                                         AAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGG                                                                                                             >  SRR19973065.3234172/1‑150 (MQ=60)
                                          AACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGG                                                                                                            <  SRR19973065.405369/150‑1 (MQ=60)
                                           ACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGC                                                                                                           >  SRR19973065.2448935/1‑150 (MQ=60)
                                                     AATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCG                                                                                                 >  SRR19973065.3419018/1‑150 (MQ=60)
                                                        CGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCT                                                                                              >  SRR19973065.1223099/1‑150 (MQ=60)
                                                        CGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCT                                                                                              >  SRR19973065.918403/1‑150 (MQ=60)
                                                         GCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTA                                                                                             <  SRR19973065.3288948/150‑1 (MQ=60)
                                                           ACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATG                                                                                           >  SRR19973065.3362090/1‑150 (MQ=60)
                                                            CGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGG                                                                                          >  SRR19973065.97877/1‑150 (MQ=60)
                                                                TTTGCTTATGTCGGCGATTTTAGCCGTATTACCGACGATCCTCAGTGGGAAAAAGTCTACGTCGAGCGCTTTGTTCGCCATATCCTCGCGCAG‑TAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGT                                                                                      <  SRR19973065.1479526/150‑1 (MQ=60)
                                                                TTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGT                                                                                      >  SRR19973065.2282372/1‑150 (MQ=60)
                                                                TTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGT                                                                                      >  SRR19973065.1051474/1‑150 (MQ=60)
                                                                     TAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAT                                                                                 <  SRR19973065.1670614/150‑1 (MQ=60)
                                                                     TAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACAT                                                                                 <  SRR19973065.1670631/150‑1 (MQ=60)
                                                                      AATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCCG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATC                                                                                >  SRR19973065.2416012/1‑150 (MQ=60)
                                                                      AATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATC                                                                                >  SRR19973065.464731/1‑150 (MQ=60)
                                                                      AATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATC                                                                                >  SRR19973065.464740/1‑150 (MQ=60)
                                                                       ATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCC                                                                               <  SRR19973065.2942953/150‑1 (MQ=60)
                                                                       ATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCC                                                                               <  SRR19973065.2916700/150‑1 (MQ=60)
                                                                       ATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCC                                                                               <  SRR19973065.1154685/150‑1 (MQ=60)
                                                                          TCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCG                                                                            <  SRR19973065.3036985/150‑1 (MQ=60)
                                                                            GGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCG                                                                          <  SRR19973065.2739590/150‑1 (MQ=60)
                                                                              CGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGAT                                                                        >  SRR19973065.3520535/1‑150 (MQ=60)
                                                                              CGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGAT                                                                        >  SRR19973065.1810255/1‑150 (MQ=60)
                                                                               GATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATG                                                                       <  SRR19973065.1009102/150‑1 (MQ=60)
                                                                               GATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATG                                                                       >  SRR19973065.1625641/1‑150 (MQ=60)
                                                                                ATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGT                                                                      >  SRR19973065.1850746/1‑150 (MQ=60)
                                                                                     AGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCAT                                                                 <  SRR19973065.3083822/150‑1 (MQ=60)
                                                                                                CGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGC                                                      <  SRR19973065.2548856/150‑1 (MQ=60)
                                                                                                CGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGC                                                      <  SRR19973065.2874754/150‑1 (MQ=60)
                                                                                                 ttccgatctCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCG                                                     <  SRR19973065.630607/141‑1 (MQ=60)
                                                                                                 ttccgatctCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCG                                                     <  SRR19973065.630808/141‑1 (MQ=60)
                                                                                                 GACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCG                                                     <  SRR19973065.942499/150‑1 (MQ=60)
                                                                                                 GACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCG                                                     >  SRR19973065.2715748/1‑150 (MQ=60)
                                                                                                   CGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCT                                                   <  SRR19973065.3126518/150‑1 (MQ=60)
                                                                                                   CGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCT                                                   <  SRR19973065.1370679/150‑1 (MQ=60)
                                                                                                     ATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAGAAAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCAATGGGTGTAAAAACCCCGCGAGGTACCATGCGGCGAAAGCGCTG                                                  >  SRR19973065.1510155/1‑150 (MQ=60)
                                                                                                     ATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGG                                                 >  SRR19973065.1509157/1‑150 (MQ=60)
                                                                                                        CGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATG                                              <  SRR19973065.1996827/150‑1 (MQ=60)
                                                                                                        CGCAGTGGGAAAAAGTCTACGTCCAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTCTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATG                                              >  SRR19973065.2172077/1‑150 (MQ=60)
                                                                                                           AGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACA                                           >  SRR19973065.2255694/1‑150 (MQ=60)
                                                                                                           AGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGAGAAAGCGCTGGATGACA                                           <  SRR19973065.2312544/150‑1 (MQ=60)
                                                                                                            GTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACAC                                          >  SRR19973065.1875631/1‑150 (MQ=60)
                                                                                                              GGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGC                                        >  SRR19973065.2165987/1‑150 (MQ=60)
                                                                                                                  AAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACT                                    >  SRR19973065.2648521/1‑150 (MQ=60)
                                                                                                                  AAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACT                                    >  SRR19973065.2648474/1‑150 (MQ=60)
                                                                                                                  AAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACT                                    >  SRR19973065.3401454/1‑150 (MQ=60)
                                                                                                                     AGTCTACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGT                                 <  SRR19973065.1682466/150‑1 (MQ=60)
                                                                                                                          ACGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATT                            <  SRR19973065.2815795/150‑1 (MQ=60)
                                                                                                                           CGTCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTA                           <  SRR19973065.362397/150‑1 (MQ=60)
                                                                                                                             TCGAGCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGGGGCGAAAGCGCTGGTTGACACGCGACTGGTGCATTACG                         <  SRR19973065.1823223/150‑1 (MQ=60)
                                                                                                                                 GCGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGA                     <  SRR19973065.3325998/150‑1 (MQ=60)
                                                                                                                                  CGCATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAA                    <  SRR19973065.2215356/150‑1 (MQ=60)
                                                                                                                                     ATTGTTCGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGagatcggaagagcacacg                 >  SRR19973065.1678776/1‑132 (MQ=60)
                                                                                                                                           CGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGAT           <  SRR19973065.3411381/150‑1 (MQ=60)
                                                                                                                                           CGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGAT           >  SRR19973065.2633955/1‑150 (MQ=60)
                                                                                                                                           CGCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGAT           <  SRR19973065.2573759/150‑1 (MQ=60)
                                                                                                                                            GCCATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATG          >  SRR19973065.3256840/1‑150 (MQ=60)
                                                                                                                                              CATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCT        <  SRR19973065.1032999/150‑1 (MQ=60)
                                                                                                                                               ATATCCTCGCGCAT‑AAAAACCATCCGTAGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTG       <  SRR19973065.1993787/150‑1 (MQ=60)
                                                                                                                                               ATATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTG       >  SRR19973065.2896958/1‑150 (MQ=60)
                                                                                                                                                 ATCCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGAA     >  SRR19973065.3200207/1‑150 (MQ=60)
                                                                                                                                                   CCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGAAGT   >  SRR19973065.1618737/1‑150 (MQ=60)
                                                                                                                                                   CCTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGAAGT   <  SRR19973065.2071434/150‑1 (MQ=60)
                                                                                                                                                    CTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGAAGTG  <  SRR19973065.1691880/150‑1 (MQ=60)
                                                                                                                                                    CTCGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGAAGTG  <  SRR19973065.200147/150‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                                     
TTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCAGTGGGAAAAAGTCTACGTCGAGCGCATTGTTCGCCATATCCACGCGCAG‑AAAAACCATCCGTCGATCATCATCTGGTCGCTGGGCAATGAATCCGGCTATGGCTGTAACATCCGCGCGATGTACCATGCGGCGAAAGCGCTGGATGACACGCGACTGGTGCATTACGAAGAAGATCGCGATGCTGAAGTG  >  CP009273/3217140‑3217437

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