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 F247 I0 R1 1496 299.8 9582018 97.8% 9371213 149.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,004,533 G→T 89.3% T151T (ACG→ACT)  yedL → putative acyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,004,5330GT89.3% 447.3 / 34.7 202T151T (ACG→ACT) yedLputative acyl transferase
Reads supporting (aligned to +/- strand):  ref base G (15/9);  new base T (61/104);  total (89/113)
Fisher's exact test for biased strand distribution p-value = 2.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 204 total aligned reads displayed.

AGATGCCTTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCAGTGA  >  CP009273/2004387‑2004678
                                                                                                                                                  |                                                                                                                                                  
aGATGCCTTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCAATGActct                                                                                                                                                  >  2:2923066/1‑149 (MQ=255)
 gATGCCTTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTa                                                                                                                                                 >  2:3844534/1‑149 (MQ=255)
   tGCCTTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTTTCATTGAAGCAATg                                                                                                                                               >  1:846564/1‑149 (MQ=255)
      cTTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAATTTACTCTATGCGa                                                                                                                                            >  1:3550522/1‑149 (MQ=255)
       ttCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGaa                                                                                                                                           <  1:1096405/149‑1 (MQ=255)
         cATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGaaaa                                                                                                                                         <  1:1816013/149‑1 (MQ=255)
          aTCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAaat                                                                                                                                        <  1:4753209/149‑1 (MQ=255)
             gCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAaataat                                                                                                                                     <  2:4564333/149‑1 (MQ=255)
             gCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAaataat                                                                                                                                     <  2:4008972/149‑1 (MQ=255)
             gCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCAATGTCTCTATTCGAAaataat                                                                                                                                     >  1:2387283/1‑149 (MQ=255)
              ccTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATg                                                                                                                                    <  2:4727492/149‑1 (MQ=255)
               cTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGa                                                                                                                                   <  1:151517/149‑1 (MQ=255)
                tGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGaa                                                                                                                                  <  1:698772/149‑1 (MQ=255)
                tGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCGAAAATAATGaa                                                                                                                                  <  2:2857901/149‑1 (MQ=255)
                 gAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAAt                                                                                                                                 <  1:1465662/149‑1 (MQ=255)
                    cTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGc                                                                                                                              <  1:3402343/149‑1 (MQ=255)
                    cTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGc                                                                                                                              <  1:2707144/149‑1 (MQ=255)
                       aGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CAt                                                                                                                          <  1:4273744/149‑1 (MQ=255)
                         aaaCGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGa                                                                                                                        <  1:3245853/149‑1 (MQ=255)
                         aaaCGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAATTTTATCTAATCGAAAAAAATTAATTGA‑CATGa                                                                                                                        >  1:4668275/1‑149 (MQ=255)
                           aCGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAATTTACAATATTCGAAAAAAATTAATTTC‑CATGaac                                                                                                                      >  2:2072156/1‑147 (MQ=255)
                            cGGGGACACAGGCTGGATTTACGGCCTGCCATCAAATTTATTTGAAGCTTGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGcc                                                                                                                     <  1:3615026/149‑1 (MQ=255)
                             gggTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCa                                                                                                                    <  2:2904729/149‑1 (MQ=255)
                             gggTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCa                                                                                                                    <  2:3676441/149‑1 (MQ=255)
                             gggTACACAGGCTGGAGTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCa                                                                                                                    <  2:2904196/149‑1 (MQ=255)
                                tACACAGGCTGGATTTACGGCCTGCCCTCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTGTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCCTTGCCACACTGCGAAAATAATGAATTGC‑CAAGAGCCCGAc                                                                                                                 >  2:3262584/1‑149 (MQ=255)
                                tACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGAc                                                                                                                 <  1:1267633/149‑1 (MQ=255)
                                tACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGAc                                                                                                                 <  1:1266120/149‑1 (MQ=255)
                                 acacAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCAATGTCACAAAGCGAAAAAAAATAATTGA‑CATGATCCacaaa                                                                                                                >  2:4409231/1‑145 (MQ=255)
                                 acacAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGAGTTTTTTAATTGACGCGATTCGAAAAAAATTAATTTC‑CCATAGCCAGAcc                                                                                                                >  1:565794/1‑148 (MQ=255)
                                   acaGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACgca                                                                                                              <  2:4256375/149‑1 (MQ=255)
                                      ggCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTTTTATTTTTACTATTCGAAAAAAAAGGATTTC‑CATGTGCCCGACcccgca                                                                                                           >  1:1651733/1‑149 (MQ=255)
                                          ggATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGAGTTTTATAATTGACACTATGCGAAAATAATGAATTTCACAT‑AGCCAGACGCAGCACAtt                                                                                                       >  2:3419597/1‑149 (MQ=255)
                                            aTTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAAAGAATTGC‑CATGAGCCAGACGCAGCACATTCt                                                                                                     >  2:197357/1‑149 (MQ=255)
                                            aTTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAAATGACTCTATTCGAAAAAAAAGAATTGC‑CATGAGCCAGACGCAGcacaatat                                                                                                     >  1:305814/1‑145 (MQ=255)
                                             tttACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCtt                                                                                                    <  2:2415626/149‑1 (MQ=255)
                                              ttACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTg                                                                                                   <  2:164375/149‑1 (MQ=255)
                                              ttACGGCCTGCCAGCAACTTGATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTg                                                                                                   <  2:4514589/149‑1 (MQ=255)
                                                 cGGCCTGCCATCAACTTTATTTGAAGCCTGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCAt                                                                                                <  2:2913320/149‑1 (MQ=255)
                                                      tGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGAc                                                                                           <  2:3656040/149‑1 (MQ=255)
                                                      tGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGAc                                                                                           <  1:4785012/149‑1 (MQ=255)
                                                      tGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAAATGAAAATATGCGAAAAAAATTAATTTA‑CAAGAGCCAGACGCAGCACAAAATTGCATTCGAc                                                                                           >  1:1434804/1‑149 (MQ=255)
                                                       gCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTTTAAATTACACTATTCGAAAATAATTAAATTC‑CAGGAGCCCGACGCAGCACAATCTTTCAATAGACg                                                                                          >  1:2367608/1‑149 (MQ=255)
                                                       gCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTTAATCGATTCGAAAATAATTAATTTC‑CATGAGCCCGACGCAGCACAATAATTCATTAGACg                                                                                          >  1:3476579/1‑149 (MQ=255)
                                                       gCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCAAACAGTCGATTTTTATCATTGACGCTATTCGAAAATAATTAATTTA‑CATGAGCCAGACCCAGCACAATCTTGCATTCGaag                                                                                          >  2:4545063/1‑147 (MQ=255)
                                                         cATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAAATTACTCAATTCGAAAAAAATGAAATTG‑CAAGAGCCAGACCCAGCCCAATCTTTCAATAGAccgg                                                                                        >  1:2633307/1‑146 (MQ=255)
                                                           tCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACACTATGCGAAAAAAATGAAATTC‑CATGAGCCAGACGCAGCACATTAATGCATTCGACGtgca                                                                                      >  2:3540668/1‑148 (MQ=255)
                                                                   aTTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCttt                                                                              <  1:2216022/149‑1 (MQ=255)
                                                                   aTTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCttt                                                                              <  1:2216342/149‑1 (MQ=255)
                                                                   aTTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCttt                                                                              <  1:4025278/149‑1 (MQ=255)
                                                                   aTTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAAATGAATCTAATCGAAAAAAATGAATTTC‑CATGAGCCAGTAGCAGCACAATCATTCCTTAGACGTTCTGCGCCgtt                                                                              >  1:4354796/1‑149 (MQ=255)
                                                                    tttGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTTAATAGATGCGAAAATAATGAATTTC‑CATGAGCCAGACGCAGCACAATAATGCATTAGACCTGCTGCGCCtttt                                                                             >  2:2888612/1‑148 (MQ=255)
                                                                      tGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCtttatt                                                                           <  1:881422/149‑1 (MQ=255)
                                                                         aGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTttatttat                                                                        <  2:1852802/149‑1 (MQ=255)
                                                                         aGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTttatttat                                                                        <  2:2246215/149‑1 (MQ=255)
                                                                          gCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAATTGACGCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATc                                                                       >  2:2762801/1‑149 (MQ=255)
                                                                           cATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCa                                                                      <  1:3844780/149‑1 (MQ=255)
                                                                            aTGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATTAATTTC‑CATGATCCAGACGCAGCACATTCTTGGATTCGGCGTGCTGCGCCGTTATTTATccc                                                                     >  2:1387372/1‑147 (MQ=255)
                                                                             tGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGAATCTATGCGAAAATAATTTATTGG‑CATGAGCCAGACGCCGCCCATTCTTTCCTTCGACGTGGTGCGCCTTTATTTATAAcc                                                                    >  2:603692/1‑149 (MQ=255)
                                                                              ggTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACca                                                                   <  2:1176026/149‑1 (MQ=255)
                                                                               gTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATGATTGACTAGAATCGAAAATAATGAATTTA‑CATGTTCCATAAGCAGCACATTCTTGCCTTCGAAGTGCTGGGCCATTATTTATAAAcag                                                                  >  2:200776/1‑149 (MQ=255)
                                                                               gTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTAATCGAAAATAATGAAATGC‑CATGATCCAGAAGCAGCACAATATTGCAGTCGACGTGCTGCGCCGttatttataaacac                                                                  >  2:3682907/1‑143 (MQ=255)
                                                                                   cGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAATTGACTCTATGCGAAAATAATTAATTGC‑CATGAGCCAGACGCAGCACAATATTTCATTCGACCTGCTGCGCCATTATTTATAAACAACAgg                                                              >  2:3339065/1‑149 (MQ=255)
                                                                                      tGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTTACTCTATTCGAAAAAAATGAATTGC‑CATGAGCCAGACGCAGCACAATATTGCATTCGAAGTGCTGCGCCGTTATTTATCACCAGCAGGaaa                                                           >  2:644354/1‑149 (MQ=255)
                                                                                       gATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAAc                                                          <  2:1089271/149‑1 (MQ=255)
                                                                                          tGCGATCCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGcc                                                       <  2:3527452/149‑1 (MQ=255)
                                                                                           gCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCt                                                      <  2:86846/149‑1 (MQ=255)
                                                                                             gAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTg                                                    <  2:1661601/149‑1 (MQ=255)
                                                                                             gAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTg                                                    <  1:4597101/149‑1 (MQ=255)
                                                                                              aaCCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGt                                                   <  1:4644661/149‑1 (MQ=255)
                                                                                              aaCCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAAATGACTCAAAGCGAAAAAAAATAAATGC‑CCAGAGCCACACCCAGCACAATCTTGCATTAGACCTGCTGCGCCATTATTTTTAACCAACACGAAACGCCCTGt                                                   >  2:1215376/1‑149 (MQ=255)
                                                                                                ccGTTTGCTGATTATCGCCTTGATCCACCCAGTCGATTTTTATCAGTGACTCTATGCGAAAATAAAGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTcc                                                 <  2:3957187/149‑1 (MQ=255)
                                                                                                ccGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTcc                                                 <  1:2307820/149‑1 (MQ=255)
                                                                                                 cGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCa                                                <  1:4502225/149‑1 (MQ=255)
                                                                                                   tttGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAATTTAATATAATCGAAAATAATTTAATGA‑CAAGATCCAGACGCAGCACAATGTTGAATTAGTCGTGCTGCGCCTTTATTTATAAACAGCAGGAAACGCCTTGTCCATa                                              >  1:1174833/1‑149 (MQ=255)
                                                                                                    ttGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATTCGAAAATAATGAATTTC‑CATGAGCCAGACGCAGCACATTAATGCATTCGACGTGCTGCGCCTTTATTTATAACCAGCAGGAAACGCCTTGTCCAAAg                                             >  2:1558746/1‑149 (MQ=255)
                                                                                                    ttGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAg                                             <  2:3750224/149‑1 (MQ=255)
                                                                                                     tGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGa                                            <  2:2686794/149‑1 (MQ=255)
                                                                                                     tGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTAATCGAAAATAATTAAATGC‑CAAGAGCCAGACGCAGCAAATTCTTGCATTCGACGTTCTGCGCCTTTTTTTATCAACAGCCGGAAACGCCGTGTCAATAGa                                            >  2:4684707/1‑149 (MQ=255)
                                                                                                           ttATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGcccc                                      <  1:3310374/149‑1 (MQ=255)
                                                                                                             aTCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAAAAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATAACCAGCAGGAAACGCCTTTTCCATAGACGccccgt                                    >  2:776531/1‑147 (MQ=255)
                                                                                                              tCGACTTGATCCACACAGTCGATTTTTATAATTGACTCGAAGCGAAAATAATTAAATGC‑CATGTGCCAGACGCAGCACATTCTTTCAGTCGACGTGCTTCGCCGTTATTTATCACCACCAGGAAACGCCGTGTCCATAGACGCCCCATc                                   >  1:4342585/1‑149 (MQ=255)
                                                                                                                   ttGATCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACAt                              <  2:850174/149‑1 (MQ=255)
                                                                                                                     gATCCACACAGTCGATTTTTATAATTTACGCTATGCGAAAATAATGAATTTC‑CATGAGCCAGACGCAGCACAATCTTTCATTCGACGTGCAGCGCCTTTATTTATCAACAGCAGGAAACGCCTTGTACATATACTCCCCATCCACATgc                            >  2:3159275/1‑149 (MQ=255)
                                                                                                                       tCCACACAGTCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATgcgc                          <  2:1023829/149‑1 (MQ=255)
                                                                                                                          acacaGTCGATTTTTATCAATGACGCTAAGCGAAAATAATTAATTGC‑CATGAGCCAGAAGCAGCACATTATTGAATTCGACGTGCTGCGCCTTTATTTATCAACAGCAGGAAAAGCCTTGTACATATAAGCCCCTTACACAAGCGCcac                       >  2:4203452/1‑149 (MQ=255)
                                                                                                                                tCGATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGaaa                 >  1:1000583/1‑149 (MQ=255)
                                                                                                                                  gATTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAAcc               >  2:2171844/1‑149 (MQ=255)
                                                                                                                                   aTTTTTATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACct              >  1:618610/1‑149 (MQ=255)
                                                                                                                                    tttttATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCGCATGCGCCACAAGAAACctc             <  2:3978505/149‑1 (MQ=255)
                                                                                                                                    tttttATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACctc             >  1:2743550/1‑149 (MQ=255)
                                                                                                                                     ttttATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACctct            <  1:2911393/149‑1 (MQ=255)
                                                                                                                                      tttATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTa           <  1:4201876/149‑1 (MQ=255)
                                                                                                                                        tATCATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTAtt         >  2:4374938/1‑149 (MQ=255)
                                                                                                                                           cATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCa      >  1:2308299/1‑149 (MQ=255)
                                                                                                                                           cATTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCa      >  1:4698490/1‑149 (MQ=255)
                                                                                                                                            aTTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCAg     <  2:2884025/149‑1 (MQ=255)
                                                                                                                                            aTTGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCAg     <  2:2883974/149‑1 (MQ=255)
                                                                                                                                              tGACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTAGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCAGTg   >  1:1346097/1‑149 (MQ=255)
                                                                                                                                               gACTCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCAGTGa  >  2:4499082/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGATGCCTTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCGAAAATAATGAATTGC‑CATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCAGTGA  >  CP009273/2004387‑2004678

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ 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
AGATGCCTTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCAGTGACAC  >  CP009273/2004387‑2004681
                                                                                                                                                  |                                                                                                                                                     
AGATGCCTTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTA                                                                                                                                                    >  SRR19973071.3125636/1‑150 (MQ=60)
  ATGCCTTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTTTCATTGAAGCAATG                                                                                                                                                  >  SRR19973071.849565/1‑150 (MQ=60)
     CCTTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAATTTACTCTATGCG‑A                                                                                                                                              >  SRR19973071.3562911/1‑150 (MQ=60)
       TTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAA                                                                                                                                            <  SRR19973071.1100275/150‑1 (MQ=60)
       TTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAA                                                                                                                                            <  SRR19973071.2888978/150‑1 (MQ=60)
       TTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCAACACAGTCGATTTTTATCATTGACTCTATGCG‑AAA                                                                                                                                            <  SRR19973071.489260/150‑1 (MQ=60)
         CATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAAT                                                                                                                                          <  SRR19973071.1822640/150‑1 (MQ=60)
         CATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAAT                                                                                                                                          <  SRR19973071.1822771/150‑1 (MQ=60)
         CATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAAT                                                                                                                                          <  SRR19973071.1822792/150‑1 (MQ=60)
          ATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATA                                                                                                                                         <  SRR19973071.4770017/150‑1 (MQ=60)
            CGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAAT                                                                                                                                       <  SRR19973071.4706236/150‑1 (MQ=60)
            CGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCAATGTCTCTATTCG‑AAAATAAT                                                                                                                                       >  SRR19973071.2395879/1‑150 (MQ=60)
               CTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAA                                                                                                                                    <  SRR19973071.4399026/150‑1 (MQ=60)
               CTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAA                                                                                                                                    <  SRR19973071.152206/150‑1 (MQ=60)
                TGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAAT                                                                                                                                   <  SRR19973071.701237/150‑1 (MQ=60)
                 GAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATT                                                                                                                                  <  SRR19973071.1471006/150‑1 (MQ=60)
                 GAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAATTTAAACTATGAG‑AAAATAATGAATT                                                                                                                                  >  SRR19973071.3650377/1‑150 (MQ=60)
                  AGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGACTTTTATCATTGACTCTATGCG‑AAAATAATGAATTG                                                                                                                                 <  SRR19973071.1207668/150‑1 (MQ=60)
                    CTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCC                                                                                                                               <  SRR19973071.2716794/150‑1 (MQ=60)
                    CTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCC                                                                                                                               <  SRR19973071.3414395/150‑1 (MQ=60)
                    CTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCC                                                                                                                               <  SRR19973071.3421992/150‑1 (MQ=60)
                    CTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCC                                                                                                                               <  SRR19973071.3415118/150‑1 (MQ=60)
                       AGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATG                                                                                                                            <  SRR19973071.4288887/150‑1 (MQ=60)
                         AAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAG                                                                                                                          <  SRR19973071.3257301/150‑1 (MQ=60)
                           ACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCC                                                                                                                        <  SRR19973071.1494727/150‑1 (MQ=60)
                            CGGGGACACAGGCTGGATTTACGGCCTGCCATCAAATTTATTTGAAGCTTGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCA                                                                                                                       <  SRR19973071.3627677/150‑1 (MQ=60)
                                TACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTTTATCTATTCG‑AAAATAATTAATTTCCATGAGCCAGACG                                                                                                                   >  SRR19973071.1461972/1‑150 (MQ=60)
                                TACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACG                                                                                                                   <  SRR19973071.1272209/150‑1 (MQ=60)
                                TACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACG                                                                                                                   <  SRR19973071.1270685/150‑1 (MQ=60)
                                TACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACG                                                                                                                   <  SRR19973071.3913611/150‑1 (MQ=60)
                                TACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGAGTTTTTTAATTGACGCgattcgaaaaaaattaatttcccatagccagacc                                                                                                                    >  SRR19973071.567916/1‑116 (MQ=60)
                                           GATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAAATGACTCTATTCG‑AAAAAAAAGAATTGCCATGAGCCAGACGCAGCACAatat                                                                                                        >  SRR19973071.306891/1‑146 (MQ=60)
                                            ATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTT                                                                                                       <  SRR19973071.2737865/150‑1 (MQ=60)
                                                ACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCAT                                                                                                   <  SRR19973071.872255/150‑1 (MQ=60)
                                                 CGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATT                                                                                                  <  SRR19973071.2667050/150‑1 (MQ=60)
                                                  GGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTC                                                                                                 <  SRR19973071.1169177/150‑1 (MQ=60)
                                                  GGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAATTTAAGCAATGCG‑AAAATAATGAATTtccaagatccagaacccgcacaatattggatta                                                                                                 >  SRR19973071.1271452/1‑117 (MQ=60)
                                                     CTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGAC                                                                                              <  SRR19973071.2249278/150‑1 (MQ=60)
                                                     CTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAAATGAAAATATGCG‑AAAAAAATTAATTTACAAGAGCCAGACGCAGCACAAAATTGCATTCGAC                                                                                              >  SRR19973071.1440073/1‑150 (MQ=60)
                                                      TGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTTAATCGATTCG‑AAAATAATTAATTTCCATGAGCCCGACGCAGCACAataattcattagacg                                                                                             >  SRR19973071.3488823/1‑135 (MQ=60)
                                                      TGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACG                                                                                             <  SRR19973071.4801878/150‑1 (MQ=60)
                                                       GCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGT                                                                                            <  SRR19973071.134493/150‑1 (MQ=60)
                                                                  TATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTTACTCTATGCG‑AAAAAAATTTATTGCCATGTGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTT                                                                                 >  SRR19973071.4370077/1‑150 (MQ=60)
                                                                   ATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTA                                                                                <  SRR19973071.4039448/150‑1 (MQ=60)
                                                                   ATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTA                                                                                <  SRR19973071.2224531/150‑1 (MQ=60)
                                                                   ATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTA                                                                                <  SRR19973071.2224362/150‑1 (MQ=60)
                                                                   ATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTA                                                                                <  SRR19973071.2224211/150‑1 (MQ=60)
                                                                      TGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATGATTTAATCTATGCGAAAAAAAATAATTTGCCA‑GAGCCAGACGCACCACCAACTTTCATTCGACGTGGTGCGCCGTTATTT                                                                             >  SRR19973071.104081/1‑150 (MQ=60)
                                                                      TGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTT                                                                             <  SRR19973071.884584/150‑1 (MQ=60)
                                                                           CATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCAC                                                                        <  SRR19973071.3858405/150‑1 (MQ=60)
                                                                            ATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAAGTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACC                                                                       <  SRR19973071.1418911/150‑1 (MQ=60)
                                                                              GGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAG                                                                     <  SRR19973071.769675/150‑1 (MQ=60)
                                                                              GGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGAGTTTTTTCATTGACGATATGCG‑AAAATAATGtatttgcctgcgccagacgccgcacagtcgtgcattcggagtggtgcgcgggtttttatcaccag                                                                     >  SRR19973071.2987143/1‑85 (MQ=60)
                                                                                       GATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACG                                                            <  SRR19973071.2377073/150‑1 (MQ=60)
                                                                                       GATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACG                                                            <  SRR19973071.2376093/150‑1 (MQ=60)
                                                                                            CGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTG                                                       <  SRR19973071.3612958/150‑1 (MQ=60)
                                                                                            CGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAATTTAATCTATTCG‑AAAAAAATTAATTTCCCAGATCCACAACCAGCACAATATTGCATTAGACGTGCTGCGCCGTTATTTATCACCAGCAGGAAACCCCATG                                                       >  SRR19973071.4320745/1‑150 (MQ=44)
                                                                                             GAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGT                                                      <  SRR19973071.4613490/150‑1 (MQ=60)
                                                                                              AACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTC                                                     <  SRR19973071.4661203/150‑1 (MQ=60)
                                                                                              AACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTC                                                     <  SRR19973071.4661415/150‑1 (MQ=60)
                                                                                               ACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCC                                                    <  SRR19973071.986804/150‑1 (MQ=60)
                                                                                                CCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCA                                                   <  SRR19973071.2316241/150‑1 (MQ=60)
                                                                                                 CGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCAT                                                  <  SRR19973071.4518227/150‑1 (MQ=60)
                                                                                                  GTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATAATTTAATATAATCG‑AAAATAATTTAATGACAAGATCCAGACGCAGCACAATGTTGAATTAGTCGTGCTGCGCCTTTATTTATAAACAGCAGGAAACGCCTTGTCCATA                                                 >  SRR19973071.1178881/1‑150 (MQ=38)
                                                                                                       CTGATTATCGACTTGATCCACACAGTCGATTTTTTTAATTTTCTATATTCG‑AAAATAATTAAATGCCATGAGCCAGACGCACCACAATCTTGCATTAGACGTGCTGCGCCGTTATTTATAACCAGCAGGAAACGCCTTGTCCAAAGAAGC                                            >  SRR19973071.2588847/1‑150 (MQ=60)
                                                                                                       CTGATTATCGACTTGATCCACACAGTCGAGTTTTATAATTGTCGCTATTCG‑AAAATAAAGAAATTCCCTGTGCCAGACGCCGCCCAATCTTGCAATCGACGTGCTGGGCCGTTATTTATAACCAGCACGAAACCCCTTGTCCAAAGACGC                                            >  SRR19973071.2588315/1‑150 (MQ=45)
                                                                                                          ATTATCGACTTGATCCACACAGTCGATTTTTATCAATGTAGCAATTCG‑AAAAAAATGAATTGCCAAGATCCCGACGCAGCACAATCTTTCATTAGACGTGCGGCGCCGTTATTTATAACCAGCAGGAAACcccatttacaaataaccccc                                         >  SRR19973071.1596325/1‑130 (MQ=50)
                                                                                                           TTATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCT                                        <  SRR19973071.3322109/150‑1 (MQ=60)
                                                                                                            TATCGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTT                                       <  SRR19973071.3311503/150‑1 (MQ=60)
                                                                                                             ATCGACTTGATCCACACAGTCGATTTTTATAATTGACTCGAAGCG‑AAAATAATTAAATGCCATGTGCCAGACGCAGCACATTCTTTCAGTCGACGTGCTTCGCCGTTATTTATCACCACCAGGAAACGCCGTGTCCATAGACGCCCCATC                                      >  SRR19973071.4357895/1‑150 (MQ=60)
                                                                                                               CGACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCA                                    <  SRR19973071.4788096/150‑1 (MQ=60)
                                                                                                                GACTTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCAC                                   <  SRR19973071.4155935/150‑1 (MQ=60)
                                                                                                                  CTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACAT                                 <  SRR19973071.4583427/150‑1 (MQ=60)
                                                                                                                   TTGATCCACACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATG                                <  SRR19973071.1589866/150‑1 (MQ=60)
                                                                                                                            ACAGTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAG                       <  SRR19973071.58384/150‑1 (MQ=60)
                                                                                                                               GTCGATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAA                    >  SRR19973071.1004220/1‑150 (MQ=60)
                                                                                                                               GTCGATTTTTATCATTGACGCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAA                    <  SRR19973071.4654521/150‑1 (MQ=60)
                                                                                                                                  GATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCT                 >  SRR19973071.620897/1‑150 (MQ=60)
                                                                                                                                   ATTTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTC                >  SRR19973071.2753270/1‑150 (MQ=60)
                                                                                                                                     TTTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTA              <  SRR19973071.2921824/150‑1 (MQ=60)
                                                                                                                                      TTTATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTAT             <  SRR19973071.4216830/150‑1 (MQ=60)
                                                                                                                                        TATCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTC           <  SRR19973071.649042/150‑1 (MQ=60)
                                                                                                                                          TCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCA         >  SRR19973071.4715246/1‑150 (MQ=60)
                                                                                                                                          TCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCA         >  SRR19973071.2316722/1‑150 (MQ=60)
                                                                                                                                          TCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCA         >  SRR19973071.4715172/1‑150 (MQ=60)
                                                                                                                                          TCATTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCA         >  SRR19973071.2316324/1‑150 (MQ=60)
                                                                                                                                             GTGACTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTAGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCAGTG      >  SRR19973071.1351050/1‑150 (MQ=60)
                                                                                                                                                 CTCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGACACCTCTATTCCAGTGACAC  >  SRR19973071.1545222/1‑150 (MQ=60)
                                                                                                                                                  |                                                                                                                                                     
AGATGCCTTCATCGCCTGAGCTTAGAAACGGGTACACAGGCTGGATTTACGGCCTGCCATCAACTTTATTTGAAGCATGGTTTCGCTGATTGCGAACCGTTTGCTGATTATCGACTTGATCCACACAGTCGATTTTTATCATTGACGCTATGCG‑AAAATAATGAATTGCCATGAGCCAGACGCAGCACATTCTTGCATTCGACGTGCTGCGCCTTTATTTATCACCAGCAGGAAACGCCTTGTCCATAGACGCCCCTTCCACATGCGCCACAAGAAACCTCTATTCCAGTGACAC  >  CP009273/2004387‑2004681

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