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 F368 I0 R1 1472 214.4 7116750 95.9% 6824963 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,128,128 T→A 100% K336* (AAA→TAA)  wzc ← colanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,128,1280TA95.2% 725.8 / 19.8 229K336* (AAA→TAA) wzccolanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (3/8);  new base A (95/123);  total (98/131)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 229 total aligned reads displayed.

TGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTT  >  CP009273/2127983‑2128272
                                                                                                                                                 |                                                                                                                                                
tgctgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAgg                                                                                                                                                <  2:2679555/148‑1 (MQ=255)
 ggtgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGgt                                                                                                                                               >  2:1954413/3‑148 (MQ=255)
  ctgACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGgtg                                                                                                                                              <  2:2341388/148‑1 (MQ=255)
    gACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTa                                                                                                                                            <  2:2900047/148‑1 (MQ=255)
    gACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTa                                                                                                                                            <  2:1901272/148‑1 (MQ=255)
    gACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTa                                                                                                                                            <  1:1398136/148‑1 (MQ=255)
     aCCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTAc                                                                                                                                           <  2:394336/148‑1 (MQ=255)
     aCCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTAc                                                                                                                                           <  1:267419/148‑1 (MQ=255)
     aCCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTAc                                                                                                                                           <  1:1728355/148‑1 (MQ=255)
     aCCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTAc                                                                                                                                           >  1:820374/1‑148 (MQ=255)
       cAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAg                                                                                                                                         >  1:1359951/1‑148 (MQ=255)
       cAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAg                                                                                                                                         <  1:1973682/148‑1 (MQ=255)
         gACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCt                                                                                                                                       >  1:1280226/1‑148 (MQ=255)
         gACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCt                                                                                                                                       <  2:3354134/148‑1 (MQ=255)
         gACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCt                                                                                                                                       >  1:793396/1‑148 (MQ=255)
          aCTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCtt                                                                                                                                      >  1:762186/1‑148 (MQ=255)
           cTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTg                                                                                                                                     <  1:1378213/148‑1 (MQ=255)
           cTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTg                                                                                                                                     <  1:880533/148‑1 (MQ=255)
           cTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTg                                                                                                                                     <  1:3147191/148‑1 (MQ=255)
           cTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTg                                                                                                                                     <  1:3147117/148‑1 (MQ=255)
              gACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGgag                                                                                                                                  <  1:870869/148‑1 (MQ=255)
              gACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGgag                                                                                                                                  <  1:3470696/148‑1 (MQ=255)
              gACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGgag                                                                                                                                  <  1:107074/148‑1 (MQ=255)
                cATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGAt                                                                                                                                <  1:3191528/148‑1 (MQ=255)
                  tCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAttt                                                                                                                              <  2:1499547/148‑1 (MQ=255)
                  tCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAttt                                                                                                                              <  2:2224243/148‑1 (MQ=255)
                  tCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAttt                                                                                                                              >  1:412543/1‑148 (MQ=255)
                  tCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAttt                                                                                                                              >  1:412729/1‑148 (MQ=255)
                     cgGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCg                                                                                                                           >  2:941508/1‑148 (MQ=255)
                      gggTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGc                                                                                                                          <  1:2706926/148‑1 (MQ=255)
                       ggTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGcc                                                                                                                         >  2:2595357/1‑148 (MQ=255)
                         tCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCctc                                                                                                                       <  1:3001591/148‑1 (MQ=255)
                             aCGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTa                                                                                                                   >  2:1836874/1‑148 (MQ=255)
                             aCGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTa                                                                                                                   >  2:2678918/1‑148 (MQ=255)
                                gACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAg                                                                                                                <  2:1978808/148‑1 (MQ=255)
                                gACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAg                                                                                                                >  2:3096663/1‑148 (MQ=255)
                                 aCAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAgg                                                                                                               >  2:3228231/1‑148 (MQ=255)
                                  cAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGt                                                                                                              <  2:1446497/148‑1 (MQ=255)
                                  cAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGt                                                                                                              <  2:966199/148‑1 (MQ=255)
                                   aaTCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTc                                                                                                             >  1:465647/1‑148 (MQ=255)
                                      ctcCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgt                                                                                                          >  2:1568765/1‑148 (MQ=255)
                                        ccTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttc                                                                                                        <  1:541741/148‑1 (MQ=255)
                                         ctgctgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGTTTGGAGATTTCCGCCTCTTTAAAGGTCAgttcg                                                                                                       <  2:2173697/148‑1 (MQ=255)
                                            ctgGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAgttcgttc                                                                                                    <  2:784363/148‑1 (MQ=255)
                                              gggTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCaa                                                                                                  <  2:871940/148‑1 (MQ=255)
                                              gggTTTTCGGCATAGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCaa                                                                                                  >  2:600240/1‑148 (MQ=255)
                                                 ttttCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTg                                                                                               >  2:175270/1‑148 (MQ=255)
                                                  tttCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTgc                                                                                              >  2:2827907/1‑148 (MQ=255)
                                                    tCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTgcgc                                                                                            <  2:824479/148‑1 (MQ=255)
                                                        cATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCCGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCg                                                                                        >  1:1029402/1‑148 (MQ=255)
                                                          tCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGAt                                                                                      <  1:1135543/148‑1 (MQ=255)
                                                          tCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCACCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTACGCCTCGTTAAAGGTCAGTTCGTTCAACAGCGCATCGAt                                                                                      >  1:3283858/1‑148 (MQ=255)
                                                             ccGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGtt                                                                                   <  1:1106014/148‑1 (MQ=255)
                                                             ccGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGtt                                                                                   <  2:1137194/148‑1 (MQ=255)
                                                             ccGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGtt                                                                                   <  2:1137508/148‑1 (MQ=255)
                                                              cGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTc                                                                                  >  1:982584/1‑148 (MQ=255)
                                                                 tACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCAcc                                                                               <  1:1682702/148‑1 (MQ=255)
                                                                 tACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCAcc                                                                               <  1:1839340/148‑1 (MQ=255)
                                                                 tACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCAcc                                                                               <  1:38773/148‑1 (MQ=255)
                                                                      gACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcga                                                                          <  1:1464438/148‑1 (MQ=255)
                                                                      gACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcga                                                                          <  2:3268156/148‑1 (MQ=255)
                                                                       aCCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCatcgaa                                                                         >  2:2215269/1‑148 (MQ=255)
                                                                        ccGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCATTAAAGGTCAGTTCGTTCAACTGCGCATCGAAGTTCACCatcgaat                                                                        >  1:2699423/1‑148 (MQ=255)
                                                                            tAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAg                                                                    >  2:1354466/1‑148 (MQ=255)
                                                                            tAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAg                                                                    >  1:2025336/1‑148 (MQ=255)
                                                                                   ggCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCt                                                             <  1:1906665/148‑1 (MQ=255)
                                                                                    gCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCtt                                                            >  2:2475605/1‑148 (MQ=255)
                                                                                      tttttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCtttt                                                          >  2:241414/1‑148 (MQ=255)
                                                                                       ttttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTg                                                         <  2:912111/148‑1 (MQ=255)
                                                                                       ttttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTg                                                         >  2:43916/1‑148 (MQ=255)
                                                                                         ttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCt                                                       >  1:1237985/1‑148 (MQ=255)
                                                                                         ttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCt                                                       <  2:3053175/148‑1 (MQ=255)
                                                                                         ttCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCt                                                       >  2:1666050/1‑148 (MQ=255)
                                                                                          tCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCtt                                                      >  1:3289969/1‑148 (MQ=255)
                                                                                           cGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTc                                                     >  2:1112500/1‑148 (MQ=255)
                                                                                            gTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTcc                                                    <  2:491150/148‑1 (MQ=255)
                                                                                              cTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAg                                                  >  1:898900/1‑148 (MQ=255)
                                                                                              cTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAg                                                  >  1:1740099/1‑148 (MQ=255)
                                                                                                tCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCg                                                >  2:2764448/1‑148 (MQ=255)
                                                                                                    gcgcCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAg                                            >  1:511407/1‑148 (MQ=255)
                                                                                                     cgcCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTACCAGCGGCAGa                                           >  1:1360223/1‑148 (MQ=255)
                                                                                                        cTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCa                                        <  1:584419/148‑1 (MQ=255)
                                                                                                             gTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGa                                   <  2:1029402/148‑1 (MQ=255)
                                                                                                              tttCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGaa                                  <  1:1188275/148‑1 (MQ=255)
                                                                                                               ttCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAAt                                 >  2:2294835/1‑148 (MQ=255)
                                                                                                                 ctcCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCt                               <  2:302884/148‑1 (MQ=255)
                                                                                                                   ccAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCttt                             <  2:469267/148‑1 (MQ=255)
                                                                                                                       cagcGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATcc                         <  2:1830343/148‑1 (MQ=255)
                                                                                                                            tGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACg                    <  2:1016518/148‑1 (MQ=255)
                                                                                                                              cGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGa                  <  1:3063374/148‑1 (MQ=255)
                                                                                                                                 tACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAgg               >  1:3402130/1‑148 (MQ=255)
                                                                                                                                  aCGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGGTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGc              <  1:755746/148‑1 (MQ=255)
                                                                                                                                   cGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa             <  1:1533309/148‑1 (MQ=255)
                                                                                                                                   cGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa             <  1:2104150/148‑1 (MQ=255)
                                                                                                                                   cGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCa             <  1:311289/148‑1 (MQ=255)
                                                                                                                                     ccGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCAtt           >  2:223047/1‑148 (MQ=255)
                                                                                                                                          tGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtt      <  1:2595357/148‑1 (MQ=255)
                                                                                                                                           gAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgttt     >  1:957188/1‑148 (MQ=255)
                                                                                                                                            aaCTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtttg    >  2:621822/1‑148 (MQ=255)
                                                                                                                                             aCTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAgtttgt   >  2:2505823/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTT  >  CP009273/2127983‑2128272

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
CCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTTCGG  >  CP009273/2127981‑2128277
                                                                                                                                                   |                                                                                                                                                     
CCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGG                                                                                                                                                     <  SRR19973067.2411874/150‑1 (MQ=60)
     TGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTAC                                                                                                                                                >  SRR19973067.825987/1‑150 (MQ=60)
      GACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACA                                                                                                                                               <  SRR19973067.1407603/150‑1 (MQ=60)
      GACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACA                                                                                                                                               <  SRR19973067.1407626/150‑1 (MQ=60)
       ACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAG                                                                                                                                              >  SRR19973067.1369240/1‑150 (MQ=60)
       ACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAG                                                                                                                                              >  SRR19973067.1369241/1‑150 (MQ=60)
       ACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAG                                                                                                                                              <  SRR19973067.269621/150‑1 (MQ=60)
       ACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAG                                                                                                                                              >  SRR19973067.462356/1‑150 (MQ=60)
       ACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAG                                                                                                                                              >  SRR19973067.462322/1‑150 (MQ=60)
       ACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAG                                                                                                                                              <  SRR19973067.1740193/150‑1 (MQ=60)
       ACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAG                                                                                                                                              <  SRR19973067.269222/150‑1 (MQ=60)
         CAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCT                                                                                                                                            <  SRR19973067.1986992/150‑1 (MQ=60)
         CAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCT                                                                                                                                            >  SRR19973067.798874/1‑150 (MQ=60)
         CAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCT                                                                                                                                            >  SRR19973067.1288948/1‑150 (MQ=60)
         CAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCT                                                                                                                                            <  SRR19973067.1987055/150‑1 (MQ=60)
          AGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTT                                                                                                                                           <  SRR19973067.487534/150‑1 (MQ=60)
          AGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTT                                                                                                                                           >  SRR19973067.767478/1‑150 (MQ=60)
          AGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTT                                                                                                                                           >  SRR19973067.767636/1‑150 (MQ=60)
           GACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTG                                                                                                                                          >  SRR19973067.3114107/1‑150 (MQ=60)
           GACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTG                                                                                                                                          >  SRR19973067.1333553/1‑150 (MQ=60)
             CTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTTCGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGTA                                                                                                                                        <  SRR19973067.3167085/150‑1 (MQ=60)
             CTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGA                                                                                                                                        <  SRR19973067.1387613/150‑1 (MQ=60)
             CTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGA                                                                                                                                        <  SRR19973067.885456/150‑1 (MQ=60)
             CTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGA                                                                                                                                        <  SRR19973067.886402/150‑1 (MQ=60)
             CTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGA                                                                                                                                        <  SRR19973067.3167087/150‑1 (MQ=60)
             CTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGA                                                                                                                                        <  SRR19973067.3167088/150‑1 (MQ=60)
             CTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGA                                                                                                                                        <  SRR19973067.3167162/150‑1 (MQ=60)
               CGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGA                                                                                                                                      <  SRR19973067.2369639/150‑1 (MQ=60)
                GACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAT                                                                                                                                     <  SRR19973067.876684/150‑1 (MQ=60)
                GACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAT                                                                                                                                     <  SRR19973067.108058/150‑1 (MQ=60)
                GACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAT                                                                                                                                     <  SRR19973067.108128/150‑1 (MQ=60)
                GACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAT                                                                                                                                     <  SRR19973067.3492744/150‑1 (MQ=60)
                 ACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATT                                                                                                                                    <  SRR19973067.1744275/150‑1 (MQ=60)
                  CATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTT                                                                                                                                   <  SRR19973067.3211687/150‑1 (MQ=60)
                  CATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTT                                                                                                                                   <  SRR19973067.3211519/150‑1 (MQ=60)
                  CATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTT                                                                                                                                   >  SRR19973067.415411/1‑150 (MQ=60)
                  CATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTT                                                                                                                                   >  SRR19973067.415497/1‑150 (MQ=60)
                  CATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTT                                                                                                                                   <  SRR19973067.3307705/150‑1 (MQ=60)
                  CATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTT                                                                                                                                   >  SRR19973067.415598/1‑150 (MQ=60)
                    TCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCC                                                                                                                                 <  SRR19973067.210675/150‑1 (MQ=60)
                        GGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCT                                                                                                                             <  SRR19973067.2724349/150‑1 (MQ=60)
                        GGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCT                                                                                                                             <  SRR19973067.2724359/150‑1 (MQ=60)
                        GGGTCAGACGGACAATCTCCGGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCT                                                                                                                             <  SRR19973067.935469/150‑1 (MQ=60)
                         GGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTC                                                                                                                            <  SRR19973067.2664447/150‑1 (MQ=60)
                           TCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTT                                                                                                                          <  SRR19973067.3020644/150‑1 (MQ=60)
                                   ACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTC                                                                                                                  >  SRR19973067.468867/1‑150 (MQ=60)
                                     AATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAG                                                                                                                <  SRR19973067.1431338/150‑1 (MQ=60)
                                          CCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGT                                                                                                           <  SRR19973067.545343/150‑1 (MQ=60)
                                                        GGCATCGCCGTTACGTGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCAGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAATGTCAGTTAGTTCAACAGCGCAACG                                                                                             >  SRR19973067.1237887/1‑150 (MQ=60)
                                                        GGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTATTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTACGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCG                                                                                             >  SRR19973067.1701589/1‑150 (MQ=60)
                                                        GACATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCCGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCG                                                                                             >  SRR19973067.1036487/1‑150 (MQ=60)
                                                        GACATCGCCGTTACGCGACCGGTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTATCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTACGCCTCTTTAAATGTAAGTTAGTTCAACTGCGCAACG                                                                                             >  SRR19973067.1036494/1‑150 (MQ=60)
                                                          CATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCCGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTACGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGAT                                                                                           >  SRR19973067.3305055/1‑150 (MQ=60)
                                                          CATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCACCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTACGCCTCGTTAAAGGTCAGTTCGTTCAACAGCGCATCGAT                                                                                           >  SRR19973067.3304460/1‑150 (MQ=60)
                                                           ATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGTGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATG                                                                                          <  SRR19973067.3163165/150‑1 (MQ=60)
                                                           ATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATG                                                                                          <  SRR19973067.2764121/150‑1 (MQ=60)
                                                            TCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGT                                                                                         <  SRR19973067.1178020/150‑1 (MQ=60)
                                                            TCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGT                                                                                         <  SRR19973067.1143193/150‑1 (MQ=60)
                                                              GCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTC                                                                                       >  SRR19973067.989096/1‑150 (MQ=60)
                                                              GCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTC                                                                                       >  SRR19973067.989331/1‑150 (MQ=60)
                                                               CCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCA                                                                                      <  SRR19973067.1113526/150‑1 (MQ=60)
                                                                   TACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCAT                                                                                  <  SRR19973067.1694307/150‑1 (MQ=60)
                                                                   TACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCAT                                                                                  <  SRR19973067.39328/150‑1 (MQ=60)
                                                                   TACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCAT                                                                                  <  SRR19973067.1851883/150‑1 (MQ=60)
                                                                   TACGCGACCGTTAAGTTTGGCTTTTTCGTCTTACAGCGCCTGACCTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGAATTCCGCCACTTTAAATGTCAGTTCGTTCAACTGCGCATCGATGTTCACCAT                                                                                  >  SRR19973067.2796387/1‑150 (MQ=60)
                                                                       CGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAA                                                                              >  SRR19973067.1098950/1‑150 (MQ=60)
                                                                        GACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAAT                                                                             >  SRR19973067.1454631/1‑150 (MQ=60)
                                                                        GACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAAT                                                                             <  SRR19973067.1474216/150‑1 (MQ=60)
                                                                        GACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAAT                                                                             <  SRR19973067.40017/150‑1 (MQ=60)
                                                                        GACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCATTAAAGGTCAGTTCGTTCAACTGCGCATCGAAGTTCACCATCGAAT                                                                             >  SRR19973067.2716789/1‑150 (MQ=60)
                                                                         ACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATC                                                                            <  SRR19973067.891516/150‑1 (MQ=60)
                                                                           CGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGA                                                                          <  SRR19973067.907785/150‑1 (MQ=60)
                                                                            GTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAG                                                                         >  SRR19973067.2038928/1‑150 (MQ=60)
                                                                                 GTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCG                                                                    >  SRR19973067.2025529/1‑150 (MQ=60)
                                                                                  TTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGC                                                                   >  SRR19973067.1913354/1‑150 (MQ=60)
                                                                                    TGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTT                                                                 >  SRR19973067.153798/1‑150 (MQ=60)
                                                                                    TGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTT                                                                 >  SRR19973067.2419472/1‑150 (MQ=60)
                                                                                     GGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTT                                                                <  SRR19973067.1919664/150‑1 (MQ=60)
                                                                                         TTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCT                                                            <  SRR19973067.6563/150‑1 (MQ=60)
                                                                                         TTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCT                                                            <  SRR19973067.6608/150‑1 (MQ=60)
                                                                                         TTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCT                                                            >  SRR19973067.1246444/1‑150 (MQ=60)
                                                                                         TTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCT                                                            >  SRR19973067.1246295/1‑150 (MQ=60)
                                                                                         TTTTCGTCGTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCT                                                            <  SRR19973067.8241/150‑1 (MQ=60)
                                                                                          TTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTT                                                           >  SRR19973067.3310647/1‑150 (MQ=60)
                                                                                           TTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTC                                                          >  SRR19973067.2795244/1‑150 (MQ=60)
                                                                                           TTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGCCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTC                                                          <  SRR19973067.808155/150‑1 (MQ=60)
                                                                                              GTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAG                                                       >  SRR19973067.904831/1‑150 (MQ=60)
                                                                                              GTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAG                                                       >  SRR19973067.136647/1‑150 (MQ=60)
                                                                                              GTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAG                                                       >  SRR19973067.1751991/1‑150 (MQ=60)
                                                                                              GTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAG                                                       >  SRR19973067.137240/1‑150 (MQ=60)
                                                                                                CTTCCAGCGCCTTACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCG                                                     <  SRR19973067.516913/150‑1 (MQ=60)
                                                                                                CTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCG                                                     <  SRR19973067.1140000/150‑1 (MQ=60)
                                                                                                    CAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAG                                                 >  SRR19973067.514893/1‑150 (MQ=60)
                                                                                                     AGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGA                                                >  SRR19973067.3540361/1‑150 (MQ=60)
                                                                                                     AGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTACCAGCGGCAGA                                                >  SRR19973067.1369514/1‑150 (MQ=60)
                                                                                                        GCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCA                                             <  SRR19973067.3259652/150‑1 (MQ=60)
                                                                                                         CCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAA                                            >  SRR19973067.1195268/1‑150 (MQ=60)
                                                                                                          CTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAAC                                           <  SRR19973067.588160/150‑1 (MQ=60)
                                                                                                            GACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAG                                         >  SRR19973067.839720/1‑150 (MQ=60)
                                                                                                               GTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGGGCGTCGATGTTCACCATCGAATCGAGCACCGCTTTCGCTTCCAGCGGAAGATCAACAGAAT                                      <  SRR19973067.1967709/150‑1 (MQ=60)
                                                                                                                TTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATC                                     <  SRR19973067.1196096/150‑1 (MQ=60)
                                                                                                                TTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATC                                     <  SRR19973067.3069629/150‑1 (MQ=60)
                                                                                                                 TTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCT                                    <  SRR19973067.3424153/150‑1 (MQ=60)
                                                                                                                    TCCAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTA                                 <  SRR19973067.1733527/150‑1 (MQ=60)
                                                                                                                      CAGCAGCGTGCGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATC                               <  SRR19973067.546562/150‑1 (MQ=60)
                                                                                                                                CGGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAG                     <  SRR19973067.3082874/150‑1 (MQ=60)
                                                                                                                                 GGTACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGG                    >  SRR19973067.3423786/1‑150 (MQ=60)
                                                                                                                                    ACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCAT                 <  SRR19973067.761003/150‑1 (MQ=60)
                                                                                                                                    ACGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGGTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCAT                 <  SRR19973067.761002/150‑1 (MQ=60)
                                                                                                                                     CGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATT                <  SRR19973067.3315000/150‑1 (MQ=60)
                                                                                                                                     CGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATT                <  SRR19973067.1543676/150‑1 (MQ=60)
                                                                                                                                     CGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATT                <  SRR19973067.2118114/150‑1 (MQ=60)
                                                                                                                                     CGCCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATT                <  SRR19973067.313337/150‑1 (MQ=60)
                                                                                                                                       CCGGGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCA              >  SRR19973067.1644688/1‑150 (MQ=60)
                                                                                                                                          GGTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTT           >  SRR19973067.3509303/1‑150 (MQ=60)
                                                                                                                                           GTGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTT          >  SRR19973067.963744/1‑150 (MQ=60)
                                                                                                                                            TGAACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTG         <  SRR19973067.2612242/150‑1 (MQ=60)
                                                                                                                                               ACTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTT      >  SRR19973067.116941/1‑150 (MQ=60)
                                                                                                                                                CTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTT     >  SRR19973067.882912/1‑150 (MQ=60)
                                                                                                                                                CTTAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTT     <  SRR19973067.176674/150‑1 (MQ=60)
                                                                                                                                                  TAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTTCG   >  SRR19973067.3184042/1‑150 (MQ=60)
                                                                                                                                                  TAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTTCG   >  SRR19973067.3184038/1‑150 (MQ=60)
                                                                                                                                                  TAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTTCG   <  SRR19973067.981887/150‑1 (MQ=60)
                                                                                                                                                  TAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTTCG   >  SRR19973067.1435827/1‑150 (MQ=60)
                                                                                                                                                  TAGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTTCG   >  SRR19973067.3176347/1‑150 (MQ=60)
                                                                                                                                                   AGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTTCGG  <  SRR19973067.3187334/150‑1 (MQ=60)
                                                                                                                                                   |                                                                                                                                                     
CCTGCTGACCAGACTCGACATCGCGGGTCAGACGGACAATCTCCTGCTGGGTTTTCGGCATCGCCGTTACGCGACCGTTAAGTTTGGCTTTTTCGTCTTCCAGCGCCTGACGTTTCTCCAGCAGCGTGCGGTACGCCGGGTGAACTTTGGTGTACAGCTTGGAGATTTCCGCCTCTTTAAAGGTCAGTTCGTTCAACTGCGCATCGATGTTCACCATCGAATCGAGCACCGCTTTTGCTTCCAGCGGCAGATCAACAGAATCTTTATCCTGACGGAAGGCATTCAGTTTGTTTTCGG  >  CP009273/2127981‑2128277

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