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 F2 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 444,519 G→A R527C (CGT→TGT)  cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273444,5190GA100.0% 245.4 / NA 69R527C (CGT→TGT) cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (45/24);  total (45/24)

ACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGG  >  NZ_CP009273/444437‑444602
                                                                                  |                                                                                   
aCGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTgg                                                                              >  2:89852/1‑90 (MQ=255)
   gCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTggtcg                                                                           >  2:601494/1‑90 (MQ=255)
     aaaGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGTCACAGTTCTggtcgcg                                                                         >  1:1180091/1‑90 (MQ=255)
      aaGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTggtcgcgg                                                                        <  2:800214/90‑1 (MQ=255)
        gTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTggtcgcggtc                                                                      >  1:978771/1‑90 (MQ=255)
          tataGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTggtcgcggtcgc                                                                    >  2:733369/1‑90 (MQ=255)
                aCGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAAt                                                                >  2:938919/1‑88 (MQ=255)
                aCGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAAt                                                                <  1:938919/88‑1 (MQ=255)
                aCGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAg                                                              >  1:376391/1‑90 (MQ=255)
                aCGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAg                                                              >  1:925899/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  2:990173/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  2:647297/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  2:181540/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  2:1080431/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  1:712908/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  1:1046376/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  1:408753/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  1:1011106/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGGCGCGGTCGCGAATAGaa                                                            >  1:594847/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  1:959403/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  1:876121/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGGCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  1:1009261/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGGCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  1:920437/1‑90 (MQ=255)
                  ggAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGGCACAGTTCTGGTCGCGGTCGCGAATAGaa                                                            >  1:586544/1‑90 (MQ=255)
                        ggggAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTAc                                                      <  2:1180091/90‑1 (MQ=255)
                           gAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATc                                                   <  2:1199458/90‑1 (MQ=255)
                           gAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATc                                                   <  1:990173/90‑1 (MQ=255)
                               agGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGaa                                               >  2:75268/1‑90 (MQ=255)
                                   tGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAAc                                           >  1:218731/1‑90 (MQ=255)
                                   tGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAAc                                           >  2:321483/1‑90 (MQ=255)
                                       cACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAgg                                       <  1:222239/90‑1 (MQ=255)
                                         cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGt                                                        >  2:479689/1‑71 (MQ=255)
                                         cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGt                                                        <  1:479689/71‑1 (MQ=255)
                                         cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                     >  1:265101/1‑90 (MQ=255)
                                         cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                     <  1:842152/90‑1 (MQ=255)
                                         cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                     <  1:915164/90‑1 (MQ=255)
                                         cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                     <  2:674829/90‑1 (MQ=255)
                                         cTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                     <  2:941767/90‑1 (MQ=255)
                                              gCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGaa                                >  1:257158/1‑90 (MQ=255)
                                              gCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGaa                                >  1:846156/1‑90 (MQ=255)
                                              gCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGaa                                >  2:584839/1‑90 (MQ=255)
                                                  aCGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAcc                            >  1:245988/1‑90 (MQ=255)
                                                  aCGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAcc                            <  2:1011106/90‑1 (MQ=255)
                                                  aCGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAcc                            <  2:277315/90‑1 (MQ=255)
                                                   cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                           <  2:1128014/90‑1 (MQ=255)
                                                   cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                           >  1:722057/1‑90 (MQ=255)
                                                   cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                           >  1:810573/1‑90 (MQ=255)
                                                   cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                           >  2:732673/1‑90 (MQ=255)
                                                   cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                           >  2:1076318/1‑90 (MQ=255)
                                                   cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                           >  2:844605/1‑90 (MQ=255)
                                                   cGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAccc                           >  1:89600/1‑90 (MQ=255)
                                                     gCCACCCCACGGGCCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAg                         >  1:734737/1‑90 (MQ=255)
                                                       cACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAgcg                       <  2:495209/90‑1 (MQ=255)
                                                       cACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAgcg                       <  1:473343/90‑1 (MQ=255)
                                                        aCCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAgcgc                      >  2:300680/1‑90 (MQ=255)
                                                        aCCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAgcgc                      >  1:1052050/1‑90 (MQ=255)
                                                          cccACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCaa                    <  2:1046376/90‑1 (MQ=255)
                                                          cccACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCaa                    <  1:800011/90‑1 (MQ=255)
                                                            cACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATc                  <  1:1189517/90‑1 (MQ=255)
                                                            cACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATc                  <  2:933649/90‑1 (MQ=255)
                                                            cACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATc                  <  1:733369/90‑1 (MQ=255)
                                                               gggTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGt               >  1:222932/1‑90 (MQ=255)
                                                               gggTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGt               >  1:925054/1‑90 (MQ=255)
                                                                   cGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                     <  1:1061673/64‑1 (MQ=255)
                                                                   cGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCa                                     >  2:1061673/1‑64 (MQ=255)
                                                                   cGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTg           >  2:176409/1‑90 (MQ=255)
                                                                      cAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCAc        <  1:433833/90‑1 (MQ=255)
                                                                           aGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCgctg   <  2:408753/90‑1 (MQ=255)
                                                                            ggTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCgctgg  >  2:1166996/1‑90 (MQ=255)
                                                                                  |                                                                                   
ACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGG  >  NZ_CP009273/444437‑444602

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCATCAGCAT  >  NZ_CP009273/444422‑444618
                                                                                                 |                                                                                                   
TGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGT                                                                                                   <  SRR3722089.486022/100‑1 (MQ=60)
  AACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTC                                                                                                 <  SRR3722089.657480/100‑1 (MQ=60)
          GCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGTCACAGTTCTGGTCGCG                                                                                         >  SRR3722089.1198070/1‑100 (MQ=60)
             CTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTC                                                                                      >  SRR3722089.993246/1‑100 (MQ=60)
                     AAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAG                                                                              >  SRR3722089.380995/1‑100 (MQ=60)
                     AAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAG                                                                              >  SRR3722089.939406/1‑100 (MQ=60)
                       GTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAA                                                                            >  SRR3722089.1026167/1‑100 (MQ=60)
                       GTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAA                                                                            >  SRR3722089.1062060/1‑100 (MQ=60)
                       GTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAA                                                                            >  SRR3722089.722184/1‑100 (MQ=60)
                       GTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAA                                                                            >  SRR3722089.413701/1‑100 (MQ=60)
                       GTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGGCGCGGTCGCGAATAGAA                                                                            >  SRR3722089.602104/1‑100 (MQ=60)
                       GTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAA                                                                            >  SRR3722089.973531/1‑100 (MQ=60)
                       GTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAA                                                                            >  SRR3722089.888724/1‑100 (MQ=60)
                       GTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGGCACAGTTCTGGTCGCGGTCGCGAATAGAA                                                                            >  SRR3722089.1024289/1‑100 (MQ=60)
                       GTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGGCACAGTTCTGGTCGCGGTCGCGAATAGAA                                                                            >  SRR3722089.933867/1‑100 (MQ=60)
                       GTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGGCGCCAGTCAGGGCACAGTTCTGGTCGCGGTCGCGAATAGAA                                                                            >  SRR3722089.593608/1‑100 (MQ=60)
                             GAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTAC                                                                      <  SRR3722089.952661/100‑1 (MQ=60)
                                     aagagacAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAAT                                                              <  SRR3722089.485393/93‑1 (MQ=60)
                                        GGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAAC                                                           >  SRR3722089.221503/1‑100 (MQ=60)
                                          GAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGA                                                         <  SRR3722089.1004869/100‑1 (MQ=60)
                                              AGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCA                                                     >  SRR3722089.268491/1‑100 (MQ=60)
                                                   GCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAA                                                >  SRR3722089.260460/1‑100 (MQ=60)
                                                   GCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAA                                                >  SRR3722089.858186/1‑100 (MQ=60)
                                                      CACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATAC                                             <  SRR3722089.225059/100‑1 (MQ=60)
                                                       ACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACC                                            >  SRR3722089.249124/1‑100 (MQ=60)
                                                        tgtgtataagagacagCCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCC                                           <  SRR3722089.1077605/84‑1 (MQ=60)
                                                        CTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCC                                           <  SRR3722089.928495/100‑1 (MQ=60)
                                                        CTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCC                                           >  SRR3722089.731452/1‑100 (MQ=60)
                                                        CTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCC                                           >  SRR3722089.90703/1‑100 (MQ=60)
                                                        CTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCC                                           <  SRR3722089.854092/100‑1 (MQ=60)
                                                        CTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCC                                           >  SRR3722089.821883/1‑100 (MQ=60)
                                                          CCAGCGTACGGCCACCCCACGGGCCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAG                                         >  SRR3722089.744323/1‑100 (MQ=60)
                                                             GCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGC                                      >  SRR3722089.1067821/1‑100 (MQ=60)
                                                                    GCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGT                               >  SRR3722089.225758/1‑100 (MQ=60)
                                                                    GCCACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGT                               >  SRR3722089.938549/1‑100 (MQ=60)
                                                                      CACCCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAC                             <  SRR3722089.478988/100‑1 (MQ=60)
                                                                         CCCACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGC                          <  SRR3722089.811070/100‑1 (MQ=60)
                                                                           CACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCAC                        <  SRR3722089.742931/100‑1 (MQ=60)
                                                                           CACGGGTCGCCAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCAC                        <  SRR3722089.1207608/100‑1 (MQ=60)
                                                                                     CAGTCAGGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGC              <  SRR3722089.439053/100‑1 (MQ=60)
                                                                                           GGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCAT        >  SRR3722089.130085/1‑100 (MQ=60)
                                                                                           GGTCACAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCAT        >  SRR3722089.701914/1‑100 (MQ=60)
                                                                                                CAGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCATCAGCA   >  SRR3722089.505248/1‑100 (MQ=60)
                                                                                                 AGTTCTGGTCGCGGTCGCGAATAGAAACGTACCTCTGAACAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCCCCGCTGGCTGCAATCATCAGCAT  <  SRR3722089.595853/100‑1 (MQ=60)
                                                                                                 AGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCATCAGCAT  <  SRR3722089.178612/100‑1 (MQ=60)
                                                                                                 AGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCATCAGCAT  <  SRR3722089.183807/100‑1 (MQ=60)
                                                                                                 |                                                                                                   
TGAACGTGCGGCACTACGGCAAAGTTATAGAACGGAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCATCAGCAT  >  NZ_CP009273/444422‑444618

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