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 F5 I102 R1 278 39.1 2251254 95.8% 2156701 87.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 834,220 A→G D148D (GAT→GAC ybiX ← PKHD‑type hydroxylase YbiX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273834,2200AG100.0% 198.8 / NA 56D148D (GAT→GACybiXPKHD‑type hydroxylase YbiX
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/29);  total (27/29)

TCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTA  >  NZ_CP009273/834132‑834305
                                                                                        |                                                                                     
tcGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTc                                                                                     <  1:936957/91‑1 (MQ=255)
     tGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTg                                                                                >  2:681068/1‑91 (MQ=255)
      ggATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGc                                                                               >  2:783678/1‑91 (MQ=255)
               tAAATGATGCCACTCGTACGCCGCGGGTTATGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGttt                                                                      <  1:622543/91‑1 (MQ=255)
               tAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGttt                                                                      <  2:98816/91‑1 (MQ=255)
               tAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGttt                                                                      <  2:562051/91‑1 (MQ=255)
               tAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGttt                                                                      <  1:783678/91‑1 (MQ=255)
                 aaTGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTa                                                                    >  1:829358/1‑91 (MQ=255)
                     atgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCg                                                                >  1:787352/1‑91 (MQ=255)
                      tgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGa                                                               >  2:628387/1‑91 (MQ=255)
                      tgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGa                                                               >  1:286685/1‑91 (MQ=255)
                      tgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGa                                                               >  2:181182/1‑91 (MQ=255)
                       gCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAt                                                              >  1:673470/1‑91 (MQ=255)
                        ccACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtg                                                             <  2:544830/91‑1 (MQ=255)
                         cACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtgt                                                            <  2:20744/91‑1 (MQ=255)
                               tACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCg                                                      <  1:961542/91‑1 (MQ=255)
                                aCGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGa                                                     <  2:829358/91‑1 (MQ=255)
                                     gcgGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtgt                                                            <  2:877733/79‑1 (MQ=255)
                                     gcgGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtgt                                                            >  1:877733/1‑79 (MQ=255)
                                     gcgGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGgtg                                                >  2:1058896/1‑91 (MQ=255)
                                     gcgGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGgtg                                                >  2:436759/1‑91 (MQ=255)
                                       gggTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTc                                              >  2:753632/1‑91 (MQ=255)
                                       gggTTACGAGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTc                                              >  1:209204/1‑91 (MQ=255)
                                               ggTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtgt                                                            <  1:775561/69‑1 (MQ=255)
                                               ggTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtgt                                                            >  2:775561/1‑69 (MQ=255)
                                                 tgtCACCGAATGCAGGCTGGTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGAcc                                    <  2:520774/91‑1 (MQ=255)
                                                  gtCACGCAATGCAGGCTGCTGGAGGGATACAGCACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCa                                   >  1:679892/1‑91 (MQ=255)
                                                  gtCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCa                                   >  1:548909/1‑91 (MQ=255)
                                                  gtCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCa                                   >  2:205522/1‑91 (MQ=255)
                                                  gtCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCa                                   >  1:555471/1‑91 (MQ=255)
                                                     aCGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGtt                                >  2:10474/1‑91 (MQ=255)
                                                      cGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTc                               <  1:177146/91‑1 (MQ=255)
                                                          aTGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGa                                                                          >  2:510643/1‑44 (MQ=255)
                                                          aTGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGa                                                                          <  1:510643/44‑1 (MQ=255)
                                                           tGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTcgccgc                          >  1:500594/1‑91 (MQ=255)
                                                           tGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTcgccgc                          <  2:382503/91‑1 (MQ=255)
                                                                 ctgctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCcgtcgt                    <  1:1012868/91‑1 (MQ=255)
                                                                   gctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAg                  <  1:205522/91‑1 (MQ=255)
                                                                   gctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAg                  <  2:22455/91‑1 (MQ=255)
                                                                   gctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAg                  <  1:595233/91‑1 (MQ=255)
                                                                    ctgGAGGGATATAACACGAGGCCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGc                 >  1:1014968/1‑91 (MQ=255)
                                                                    ctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGa                                                               >  2:803762/1‑45 (MQ=255)
                                                                    ctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGa                                                               <  1:803762/45‑1 (MQ=255)
                                                                    ctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGc                 >  1:333664/1‑91 (MQ=255)
                                                                    ctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGc                 >  1:87969/1‑91 (MQ=255)
                                                                           gATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTa                                          >  2:1083909/1‑59 (MQ=255)
                                                                           gATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTa                                          <  1:1083909/59‑1 (MQ=255)
                                                                              ataacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATc       <  2:1014968/89‑1 (MQ=255)
                                                                              acaacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATc       <  2:209204/91‑1 (MQ=255)
                                                                              acaacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATc       <  1:346240/91‑1 (MQ=255)
                                                                              acaacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATc       <  2:566050/91‑1 (MQ=255)
                                                                              acaacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATc       <  2:568392/91‑1 (MQ=255)
                                                                              acaacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATc       <  1:181182/91‑1 (MQ=255)
                                                                              acaacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATc       <  1:761345/91‑1 (MQ=255)
                                                                              acaacaCGAGGTCCCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATc       <  2:720628/91‑1 (MQ=255)
                                                                                  cacGAGGTCGCCGGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCtt   <  1:917628/91‑1 (MQ=255)
                                                                                   acGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTa  >  2:737437/1‑91 (MQ=255)
                                                                                        |                                                                                     
TCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTA  >  NZ_CP009273/834132‑834305

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAG  >  NZ_CP009273/834128‑834320
                                                                                            |                                                                                                    
ATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGC                                                                                              <  SRR3722125.700177/101‑1 (MQ=60)
    TCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGG                                                                                          <  SRR3722125.957758/101‑1 (MQ=60)
           GATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTA                                                                                   >  SRR3722125.847575/1‑101 (MQ=60)
               ctagtctcgtgggctcggagatgtgtataagagacaGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCG                                                                               <  SRR3722125.521476/65‑1 (MQ=60)
               CACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCG                                                                               >  SRR3722125.804518/1‑101 (MQ=60)
                ACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGA                                                                              >  SRR3722125.292526/1‑101 (MQ=60)
                 CATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAT                                                                             >  SRR3722125.687978/1‑101 (MQ=60)
                   TAAATGATGCCACTCGTACGCCGCGGGTTATGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGT                                                                           <  SRR3722125.635968/101‑1 (MQ=60)
                   TAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGT                                                                           <  SRR3722125.800743/101‑1 (MQ=60)
                          tagtctcgtgggctcggagatgtgtataagagacagATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGA                                                                    <  SRR3722125.821339/65‑1 (MQ=60)
                             tataagagacagGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGG                                                                 <  SRR3722125.792432/89‑1 (MQ=60)
                               CTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGct                                                               >  SRR3722125.897145/1‑99 (MQ=60)
                                 CGTACGCCGCGGGTTACGAGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTC                                                             >  SRR3722125.213336/1‑101 (MQ=60)
                                   TACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCAT                                                           <  SRR3722125.982974/101‑1 (MQ=60)
                                            GGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAGCACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCA                                                  >  SRR3722125.694557/1‑101 (MQ=60)
                                            GGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCA                                                  >  SRR3722125.560662/1‑101 (MQ=60)
                                            GGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCA                                                  >  SRR3722125.567356/1‑101 (MQ=60)
                                               cggagatgtgtataagagacaGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTT                                               <  SRR3722125.1108319/80‑1 (MQ=60)
                                                     TGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGC                                         >  SRR3722125.511191/1‑101 (MQ=60)
                                                          CGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCG                                    <  SRR3722125.180588/101‑1 (MQ=60)
                                                              ATGCAGGCTGCTGGAGGGATATAACACGAGGCCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGC                                >  SRR3722125.1037668/1‑101 (MQ=60)
                                                              ATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGC                                >  SRR3722125.340579/1‑101 (MQ=60)
                                                              ATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGC                                >  SRR3722125.844389/1‑101 (MQ=60)
                                                              ATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGC                                >  SRR3722125.89610/1‑101 (MQ=60)
                                                                     CTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGG                         <  SRR3722125.1035527/101‑1 (MQ=60)
                                                                       GCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGAT                       <  SRR3722125.608035/101‑1 (MQ=60)
                                                                       GCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGAT                       <  SRR3722125.209574/101‑1 (MQ=60)
                                                                               GATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAA               >  SRR3722125.576465/1‑101 (MQ=60)
                                                                                  ACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAAC            <  SRR3722125.777835/101‑1 (MQ=60)
                                                                                  ACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAAC            <  SRR3722125.353431/101‑1 (MQ=60)
                                                                                  ACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAAC            <  SRR3722125.184702/101‑1 (MQ=60)
                                                                                  ACAACACGAGGGCGCCTGCCGGGAGTTTTGCCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAAC            <  SRR3722125.432352/101‑1 (MQ=60)
                                                                                      CACGAGGTCGCCGGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCG        <  SRR3722125.937986/101‑1 (MQ=60)
                                                                                            GTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAG  <  SRR3722125.2169/101‑1 (MQ=60)
                                                                                            |                                                                                                    
ATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAG  >  NZ_CP009273/834128‑834320

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