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 485,143 A→G D1051G (GAT→GGT)  mscK → mechanosensitive channel MscK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273485,1430AG100.0% 221.0 / NA 59D1051G (GAT→GGT) mscKmechanosensitive channel MscK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/34);  total (25/34)

TTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGC  >  NZ_CP009273/485055‑485228
                                                                                        |                                                                                     
ttACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTc                                                                                     >  1:68186/1‑91 (MQ=255)
ttACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTc                                                                                     >  2:703684/1‑91 (MQ=255)
ttACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTc                                                                                     >  2:700266/1‑91 (MQ=255)
 tACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCa                                                                                    >  1:983045/1‑91 (MQ=255)
      aaGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGc                                                                               <  1:399574/91‑1 (MQ=255)
      aaGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTCGGGTCATGAGc                                                                               <  2:1124199/91‑1 (MQ=255)
       aGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCt                                                                              >  1:558009/1‑91 (MQ=255)
       aGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCt                                                                              >  1:163621/1‑91 (MQ=255)
        gggcgcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTg                                                                             <  1:1084436/87‑1 (MQ=255)
        ggcggcGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTg                                                                             <  2:877903/91‑1 (MQ=255)
        ggcggcGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCCGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTg                                                                             <  1:1051697/91‑1 (MQ=255)
         gcggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGc                                                                            >  1:684671/1‑91 (MQ=255)
          cggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCg                                                                           <  1:327484/91‑1 (MQ=255)
          cggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCg                                                                           <  2:216232/91‑1 (MQ=255)
            gcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTc                                                                         <  1:703684/91‑1 (MQ=255)
            gcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTc                                                                         >  2:415929/1‑91 (MQ=255)
                 tGAGCACCCAAGGGTTATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTAt                                                                    <  2:354366/91‑1 (MQ=255)
                 tGAGCACCCAAGGGTGATGCGCGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTAt                                                                    <  2:504986/91‑1 (MQ=255)
                 tGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTAt                                                                    <  2:1030385/91‑1 (MQ=255)
                 tGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTAt                                                                    <  2:260231/91‑1 (MQ=255)
                 tGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTAt                                                                    <  2:140983/91‑1 (MQ=255)
                     cACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGc                                                                <  2:629024/91‑1 (MQ=255)
                         cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGa                                                            <  1:356339/91‑1 (MQ=255)
                         cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGa                                                            >  1:315447/1‑91 (MQ=255)
                            gggTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACt                                                         >  1:702048/1‑91 (MQ=255)
                              gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGc                                                       >  1:879645/1‑91 (MQ=255)
                                 aTGCACGACCCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTg                                                    >  2:1111869/1‑91 (MQ=255)
                                 aTGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTg                                                                             <  1:595612/66‑1 (MQ=255)
                                 aTGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTg                                                                             >  2:595612/1‑66 (MQ=255)
                                 aTGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTg                                                    >  1:780703/1‑91 (MQ=255)
                                 aTGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTg                                                    >  2:740529/1‑91 (MQ=255)
                                    cACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGAcc                                                 <  2:1052196/91‑1 (MQ=255)
                                    cACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGAcc                                                 <  1:1116174/91‑1 (MQ=255)
                                      cGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGt                                               >  2:631285/1‑91 (MQ=255)
                                      cGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGt                                               >  2:315809/1‑91 (MQ=255)
                                        aaCCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAg                                             >  1:629615/1‑91 (MQ=255)
                                             aTGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTa                                        >  1:613114/1‑91 (MQ=255)
                                               gCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACt                                      <  2:999928/91‑1 (MQ=255)
                                                  ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTc                                   >  1:19434/1‑91 (MQ=255)
                                                    aaGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                 <  2:702048/91‑1 (MQ=255)
                                                         ttcttTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGc                            <  2:1104998/91‑1 (MQ=255)
                                                         ttcttTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGc                            <  2:251781/91‑1 (MQ=255)
                                                            tttACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGa                         <  2:1046266/91‑1 (MQ=255)
                                                            tttACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGa                         <  2:405060/91‑1 (MQ=255)
                                                            tttACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGa                         <  2:95353/91‑1 (MQ=255)
                                                            tttACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGa                         <  2:710488/91‑1 (MQ=255)
                                                              tACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                       <  1:1111869/91‑1 (MQ=255)
                                                              tACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                       <  2:80398/91‑1 (MQ=255)
                                                              tACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                       <  2:613114/91‑1 (MQ=255)
                                                              tACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                       <  1:137057/91‑1 (MQ=255)
                                                                  gCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTa                   >  1:519804/1‑91 (MQ=255)
                                                                    aTTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACt                 <  1:155314/91‑1 (MQ=255)
                                                                    aTTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACt                 <  1:25248/91‑1 (MQ=255)
                                                                    aTTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACt                 >  1:516504/1‑91 (MQ=255)
                                                                         gTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGa            <  2:684671/91‑1 (MQ=255)
                                                                              aGCTCGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTGCTATCGATCAGc       <  2:77109/91‑1 (MQ=255)
                                                                              aGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc       >  1:201831/1‑91 (MQ=255)
                                                                              aGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc       >  1:369259/1‑91 (MQ=255)
                                                                               gCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCt      >  2:804405/1‑91 (MQ=255)
                                                                                   gTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGc  >  2:223061/1‑91 (MQ=255)
                                                                                        |                                                                                     
TTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGC  >  NZ_CP009273/485055‑485228

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ 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
GTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATC  >  NZ_CP009273/485043‑485243
                                                                                                    |                                                                                                    
GTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGG                                                                                                      >  SRR3722125.738369/1‑101 (MQ=60)
  GCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTC                                                                                                    >  SRR3722125.69442/1‑101 (MQ=60)
   CGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCA                                                                                                   >  SRR3722125.1005000/1‑101 (MQ=60)
       AAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAG                                                                                               >  SRR3722125.337836/1‑101 (MQ=60)
         GTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCT                                                                                             >  SRR3722125.569952/1‑101 (MQ=60)
         GTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCT                                                                                             >  SRR3722125.166756/1‑101 (MQ=60)
           GTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGC                                                                                           >  SRR3722125.699444/1‑101 (MQ=60)
                  AAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTA                                                                                    <  SRR3722125.407970/101‑1 (MQ=60)
                    gtgtataagagacagCCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATG                                                                                  <  SRR3722125.608423/86‑1 (MQ=60)
                    gggcGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATG                                                                                  <  SRR3722125.1108856/97‑1 (MQ=60)
                    cgagatcatgcctagtctcgtgggctcggagatgtgtataagagacagCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATG                                                                                  <  SRR3722125.617135/53‑1 (MQ=60)
                    GGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCCGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATG                                                                                  <  SRR3722125.1075319/101‑1 (MQ=60)
                      CGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTG                                                                                <  SRR3722125.334244/101‑1 (MQ=60)
                        GCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCG                                                                              <  SRR3722125.718844/101‑1 (MQ=60)
                           ACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGA                                                                           >  SRR3722125.321929/1‑101 (MQ=60)
                              GAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACT                                                                        >  SRR3722125.717169/1‑101 (MQ=60)
                                GCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGC                                                                      >  SRR3722125.899100/1‑101 (MQ=60)
                                   CCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTG                                                                   >  SRR3722125.797706/1‑101 (MQ=60)
                                     CAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGAC                                                                 <  SRR3722125.363799/101‑1 (MQ=60)
                                          GTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAG                                                            >  SRR3722125.643180/1‑101 (MQ=60)
                                               GCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTA                                                       >  SRR3722125.626337/1‑101 (MQ=60)
                                                CACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTAC                                                      <  SRR3722125.1141368/101‑1 (MQ=60)
                                                    AACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTC                                                  >  SRR3722125.19790/1‑101 (MQ=60)
                                                                    CTTCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTA                                  >  SRR3722125.530861/1‑101 (MQ=60)
                                                                      TCTTTACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACT                                >  SRR3722125.527473/1‑101 (MQ=60)
                                                                          TACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCG                            <  SRR3722125.139675/101‑1 (MQ=60)
                                                                          TACGGCATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCG                            <  SRR3722125.1136948/101‑1 (MQ=60)
                                                                                ATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGC                      >  SRR3722125.376997/1‑101 (MQ=60)
                                                                                ATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGC                      <  SRR3722125.25690/101‑1 (MQ=60)
                                                                                ATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGC                      >  SRR3722125.205812/1‑101 (MQ=60)
                                                                                ATTTGGTGCCAGCACGTTGGGTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGC                      <  SRR3722125.158277/101‑1 (MQ=60)
                                                                                                    GTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGCACCGGACTATCGATCAGCTGTGCCGTGAAAACGACATC  >  SRR3722125.195711/1‑101 (MQ=60)
                                                                                                    GTCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATC  >  SRR3722125.1098270/1‑101 (MQ=60)
                                                                                                    |                                                                                                    
GTGCGTAAAGTGTTACTGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACGACATC  >  NZ_CP009273/485043‑485243

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