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 858,772 T→C I66V (ATT→GTT)  ybiY ← glycyl‑radical enzyme activating protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273858,7720TC100.0% 197.5 / NA 53I66V (ATT→GTT) ybiYglycyl‑radical enzyme activating protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/27);  total (26/27)

GTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCG  >  NZ_CP009273/858685‑858854
                                                                                       |                                                                                   
gtgtCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCAc                                                                                 >  2:670773/1‑91 (MQ=255)
gtgtCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCAc                                                                                 >  2:555593/1‑91 (MQ=255)
gtgtCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCAc                                                                                 >  2:673359/1‑91 (MQ=255)
gtgtCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCAc                                                                                 >  2:257990/1‑91 (MQ=255)
gtgtCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCAc                                                                                 >  1:761443/1‑91 (MQ=255)
     ggACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCg                                                                            >  2:308239/1‑91 (MQ=255)
         agcagTCGTTTAACGCCGTCAGATGCCCGGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGgcgc                                                                        <  1:916035/91‑1 (MQ=255)
         agcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGgcgc                                                                        <  1:780240/91‑1 (MQ=255)
         agcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGgcgc                                                                        <  2:491206/91‑1 (MQ=255)
         agcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGgcgc                                                                        <  2:1096000/91‑1 (MQ=255)
         agcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGgcgc                                                                        <  2:707359/91‑1 (MQ=255)
          gcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGcgcg                                                                       <  1:888927/91‑1 (MQ=255)
          gcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGcgcg                                                                       <  1:42060/91‑1 (MQ=255)
          gcagTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGcgcg                                                                       >  2:1095734/1‑91 (MQ=255)
          gcagTCCGTTAACGCCGTCAGCTGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGcgcg                                                                       <  2:172490/91‑1 (MQ=255)
            agTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGc                                                                     >  1:19023/1‑91 (MQ=255)
                   tAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAgcg                                                              <  1:127931/91‑1 (MQ=255)
                    aaCGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAgcgc                                                             >  1:505266/1‑91 (MQ=255)
                    aaCGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAgcgc                                                             >  1:544951/1‑91 (MQ=255)
                    aaCGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAgcgc                                                             <  1:402672/91‑1 (MQ=255)
                      cGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGcac                                                           <  1:569928/91‑1 (MQ=255)
                      cGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGcac                                                           <  2:835539/91‑1 (MQ=255)
                      cGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGcac                                                           <  2:761443/91‑1 (MQ=255)
                       gCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGcaca                                                          <  1:40754/91‑1 (MQ=255)
                         cGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGc                                                        >  2:411409/1‑91 (MQ=255)
                         cGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGc                                                        >  2:1026520/1‑91 (MQ=255)
                               aTGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAg                                                  <  1:308239/91‑1 (MQ=255)
                                  cTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCt                                               >  1:115559/1‑91 (MQ=255)
                                  cTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCt                                               >  2:243326/1‑91 (MQ=255)
                                    ccGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGc                                                                     <  1:145984/67‑1 (MQ=255)
                                    ccGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGc                                                                     >  2:145984/1‑67 (MQ=255)
                                    ccGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTc                                             <  1:673359/91‑1 (MQ=255)
                                    ccGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTc                                             <  1:488693/91‑1 (MQ=255)
                                           tAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGacac                                      >  1:1015999/1‑91 (MQ=255)
                                           tAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGacac                                      <  2:218647/91‑1 (MQ=255)
                                           tAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGacac                                      >  1:1068653/1‑91 (MQ=255)
                                           tAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGacac                                      >  1:843324/1‑91 (MQ=255)
                                           tAACTTTTCCCGATGAATAACCAACCCATTACGCGCCGCCCCATGCCCTCCCGGCGCGGCCTTAGAGCCCAGCTCGCAGCCTTGCAGacac                                      <  2:923377/91‑1 (MQ=255)
                                                tttCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCg                                 <  2:1098445/91‑1 (MQ=255)
                                                tttCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCg                                 <  1:1095734/91‑1 (MQ=255)
                                                     cGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGt                            <  2:70309/91‑1 (MQ=255)
                                                      gATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTc                           >  1:641366/1‑91 (MQ=255)
                                                           aTAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGcc                                                <  2:332451/65‑1 (MQ=255)
                                                           aTAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGcc                                                >  1:332451/1‑65 (MQ=255)
                                                             aaGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATAcagc                    <  2:1068653/91‑1 (MQ=255)
                                                             aaGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATAcagc                    <  2:1022101/91‑1 (MQ=255)
                                                                 aaaCCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGAt                >  1:1016818/1‑91 (MQ=255)
                                                                   aCCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGc                                                                     <  1:835530/36‑1 (MQ=255)
                                                                   aCCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGc                                                                     >  2:835530/1‑36 (MQ=38)
                                                                     cATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCtg            >  2:1000721/1‑91 (MQ=255)
                                                                      aTTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCtgc           >  1:336943/1‑91 (MQ=255)
                                                                        tCAGCGCGCGCTCAACCCCTGGCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACTCAGTCGTGCGTCATAAAGCAGATCCtgcgt         <  2:801342/91‑1 (MQ=255)
                                                                             gcgcgcTCAACCACTTCCGGCGCGGCCTTCGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCggg    >  2:426176/1‑91 (MQ=255)
                                                                              cgcgcTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCg                                 >  2:730828/1‑61 (MQ=255)
                                                                              cgcgcTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCg                                 <  1:730828/61‑1 (MQ=255)
                                                                               gcgcTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGgcg  >  1:664645/1‑91 (MQ=255)
                                                                                       |                                                                                   
GTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCG  >  NZ_CP009273/858685‑858854

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ 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
CGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTG  >  NZ_CP009273/858673‑858871
                                                                                                   |                                                                                                   
CGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACC                                                                                                    <  SRR3722125.971823/101‑1 (MQ=60)
CGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACC                                                                                                    >  SRR3722125.302931/1‑101 (MQ=60)
  GTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCAC                                                                                                  >  SRR3722125.777936/1‑101 (MQ=60)
        GCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGG                                                                                            <  SRR3722125.323901/101‑1 (MQ=60)
              GTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGC                                                                                      >  SRR3722125.19371/1‑101 (MQ=60)
                     AGCAGTCGTTTAACGCCGTCAGATGCCCGGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCG                                                                               <  SRR3722125.936349/101‑1 (MQ=60)
                     AGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCG                                                                               <  SRR3722125.797230/101‑1 (MQ=60)
                     AGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCG                                                                               <  SRR3722125.513209/101‑1 (MQ=60)
                      GCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGC                                                                              <  SRR3722125.42811/101‑1 (MQ=60)
                      GCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGC                                                                              <  SRR3722125.908596/101‑1 (MQ=60)
                      GCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGC                                                                              >  SRR3722125.556589/1‑101 (MQ=60)
                      GCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGC                                                                              >  SRR3722125.515975/1‑101 (MQ=60)
                        tgtataagagacagGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCAC                                                                            <  SRR3722125.148762/87‑1 (MQ=60)
                        gatcatgcctagtctcgtgggctcggagatgtgtataagagacagGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCAC                                                                            <  SRR3722125.853894/56‑1 (MQ=60)
                               TAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGC                                                                     <  SRR3722125.130338/101‑1 (MQ=60)
                                AACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCA                                                                    <  SRR3722125.411130/101‑1 (MQ=60)
                                  CGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGC                                                                  <  SRR3722125.582115/101‑1 (MQ=60)
                                   GCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCC                                                                 <  SRR3722125.41479/101‑1 (MQ=60)
                                    CCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCT                                                                >  SRR3722125.117732/1‑101 (MQ=60)
                                           ATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGAC                                                         <  SRR3722125.314567/101‑1 (MQ=60)
                                             GCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACAC                                                       >  SRR3722125.861877/1‑101 (MQ=60)
                                             GCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACAC                                                       >  SRR3722125.1092743/1‑101 (MQ=60)
                                             GCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACAC                                                       >  SRR3722125.1038730/1‑101 (MQ=60)
                                                CCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGT                                                    <  SRR3722125.687862/101‑1 (MQ=60)
                                                CCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGT                                                    <  SRR3722125.499014/101‑1 (MQ=60)
                                                        AACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTC                                            >  SRR3722125.655192/1‑101 (MQ=60)
                                                            gagatgtgtataagagacagCCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATAC                                        <  SRR3722125.746565/81‑1 (MQ=60)
                                                            TTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATAC                                        <  SRR3722125.1120399/101‑1 (MQ=60)
                                                             TTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACctgtctcttatacaca                                       >  SRR3722125.339334/1‑85 (MQ=60)
                                                                   ATGAATAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGAT                                 >  SRR3722125.1039580/1‑101 (MQ=60)
                                                                        TAAGCAAACCATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGC                            >  SRR3722125.343949/1‑101 (MQ=60)
                                                                                 CATTCAGCGCGCGCTCAACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCG                   >  SRR3722125.678962/1‑101 (MQ=60)
                                                                                                  ACCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTG  >  SRR3722125.273765/1‑101 (MQ=60)
                                                                                                   |                                                                                                   
CGGTTAATGCCTGTGTCGGACAGCAGTCGGTTAACGCCGTCAGATGCTCCGGGGTTAACTTTTCCCGATGAATAAGCAAACCATTCAGCGCGCGCTCAATCACTTCCGGCGCGGCCTTAGCGCACAGCTCGCAGCCTTCCAGACACAGTCGTGCGTCATACAGCAGATCCTGCGTGCGGGCGCGGCTTTCCGGGTTCTG  >  NZ_CP009273/858673‑858871

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