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 928,200 T→C R50R (CGT→CGC lrp → leucine‑responsive transcriptional regulator Lrp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273928,2000TC100.0% 182.1 / NA 50R50R (CGT→CGClrpleucine‑responsive transcriptional regulator Lrp
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/30);  total (20/30)

GTTGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGTCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGT  >  NZ_CP009273/928116‑928283
                                                                                    |                                                                                   
gTTGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTg                                                                               >  2:551017/1‑91 (MQ=255)
gTTGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTg                                                                               >  1:608765/1‑91 (MQ=255)
 ttGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTgg                                                                              >  1:242528/1‑91 (MQ=255)
  tGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGa                                                                             <  1:551017/91‑1 (MQ=255)
  tGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGa                                                                             >  1:739629/1‑91 (MQ=255)
         ggATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACaa                                                                      <  1:228694/91‑1 (MQ=255)
         agATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACaa                                                                      <  2:757651/90‑1 (MQ=255)
           aTGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAgg                                                                    >  1:1110589/1‑91 (MQ=255)
           aTGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAgg                                                                    >  1:608513/1‑91 (MQ=255)
                 gTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTAt                                                              >  1:749725/1‑91 (MQ=255)
                   aattCTAACGTCGAGCTGTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTc                                                            <  2:357018/89‑1 (MQ=255)
                   aTTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTc                                                            <  2:749725/91‑1 (MQ=255)
                   aTTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTc                                                            <  2:608513/91‑1 (MQ=255)
                   aTTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTc                                                            <  1:822489/91‑1 (MQ=255)
                   aTTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTc                                                            <  1:1080345/91‑1 (MQ=255)
                   aTTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTc                                                            <  1:680434/91‑1 (MQ=255)
                    tttCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCa                                                           >  1:81494/1‑91 (MQ=255)
                      tCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAgg                                                         <  1:647611/91‑1 (MQ=255)
                        tAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGa                                                                             <  1:225898/69‑1 (MQ=255)
                        tAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGa                                                                             >  2:225898/1‑69 (MQ=255)
                        tAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACaa                                                                      >  2:1031968/1‑76 (MQ=255)
                        tAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACaa                                                                      <  1:1031968/76‑1 (MQ=255)
                             tCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATAc                                                  <  2:407906/91‑1 (MQ=255)
                             tCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATAc                                                  <  1:249759/91‑1 (MQ=255)
                              cGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTc                                                            <  2:447744/80‑1 (MQ=255)
                              cGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTc                                                            >  1:447744/1‑80 (MQ=255)
                                  cTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGgcgc                                             <  2:887999/91‑1 (MQ=255)
                                      cTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCgctgct                                         <  2:619949/91‑1 (MQ=255)
                                              gtgGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAAccccc                                 >  2:784271/1‑91 (MQ=255)
                                              gtgGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAAccccc                                 >  1:982178/1‑91 (MQ=255)
                                               tgGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCa                                <  1:967688/91‑1 (MQ=255)
                                                gggACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCAt                               <  1:784271/91‑1 (MQ=255)
                                                     tttCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCt                          <  1:762644/91‑1 (MQ=255)
                                                     tttCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCt                          >  1:436139/1‑91 (MQ=255)
                                                      ttCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAGGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAAc                                     >  1:865510/1‑79 (MQ=255)
                                                      ttCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAGGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAAc                                     <  2:865510/79‑1 (MQ=255)
                                                              cGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATc                 <  1:831478/91‑1 (MQ=255)
                                                              cGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATc                 <  2:368910/91‑1 (MQ=255)
                                                              cGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATc                 <  1:763179/91‑1 (MQ=255)
                                                              cGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATc                 <  1:965286/91‑1 (MQ=255)
                                                                ccGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCAc               >  2:84658/1‑91 (MQ=255)
                                                                ccGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCAc               >  2:808390/1‑91 (MQ=255)
                                                                 cGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACt              <  2:1110589/91‑1 (MQ=255)
                                                                  gTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGc                                                       >  1:690797/1‑49 (MQ=255)
                                                                  gTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGc                                                       <  2:690797/49‑1 (MQ=255)
                                                                  gTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACtt             >  1:561335/1‑91 (MQ=255)
                                                                        tGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTa       <  2:886431/91‑1 (MQ=255)
                                                                        tGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTa       <  1:101859/91‑1 (MQ=255)
                                                                          aGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTAtt     <  1:307753/91‑1 (MQ=255)
                                                                             gtgtGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGt  >  2:1023489/1‑91 (MQ=255)
                                                                                    |                                                                                   
GTTGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGTCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGT  >  NZ_CP009273/928116‑928283

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
TTCTTAATGAGTTGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGTCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATT  >  NZ_CP009273/928106‑928290
                                                                                              |                                                                                          
TTCTTAATGAGTTGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTG                                                                                      >  SRR3722125.621875/1‑101 (MQ=60)
 TCTTAATGAGTTGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGG                                                                                     >  SRR3722125.247389/1‑101 (MQ=60)
  CTTAATGAGTTGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGA                                                                                    >  SRR3722125.755587/1‑101 (MQ=60)
    cagaagacggcatacgatatcatgcctagtctcgtgggctcggagatgtgtataagagacagCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAA                                                                                  <  SRR3722125.754285/39‑1 (MQ=60)
    TAATGAGTTGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGGAA                                                                                  <  SRR3722125.1149629/101‑1 (MQ=60)
           TTGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGG                                                                           >  SRR3722125.1135606/1‑101 (MQ=60)
           TTGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGG                                                                           >  SRR3722125.621619/1‑101 (MQ=60)
            tataagagacaGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGG                                                                          <  SRR3722125.230393/90‑1 (MQ=60)
            TGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGG                                                                          <  SRR3722125.562817/101‑1 (MQ=60)
                 AAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTAT                                                                     >  SRR3722125.765908/1‑101 (MQ=60)
                   gacaGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTC                                                                   <  SRR3722125.1055105/97‑1 (MQ=60)
                   GGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTC                                                                   <  SRR3722125.233255/101‑1 (MQ=60)
                    GATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCA                                                                  >  SRR3722125.83013/1‑101 (MQ=60)
                             ATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATAC                                                         <  SRR3722125.1104648/101‑1 (MQ=60)
                             ATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATAC                                                         <  SRR3722125.840522/101‑1 (MQ=60)
                             ATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATAC                                                         <  SRR3722125.695112/101‑1 (MQ=60)
                              TTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACC                                                        >  SRR3722125.457195/1‑101 (MQ=60)
                                TCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGC                                                      <  SRR3722125.661574/101‑1 (MQ=60)
                                       TCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTT                                               <  SRR3722125.254774/101‑1 (MQ=60)
                                              TTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCC                                        >  SRR3722125.1004118/1‑101 (MQ=60)
                                                     CGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCT                                 >  SRR3722125.445336/1‑101 (MQ=60)
                                                      GTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAGGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCct                                >  SRR3722125.884606/1‑99 (MQ=60)
                                                         TGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGAT                             <  SRR3722125.989274/101‑1 (MQ=60)
                                                          GGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATG                            <  SRR3722125.801354/101‑1 (MQ=60)
                                                               TTTCCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCA                       <  SRR3722125.779171/101‑1 (MQ=60)
                                                                  CCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCctgtctcttatacacatctccgagcccacgag                    >  SRR3722125.705702/1‑69 (MQ=60)
                                                                  CCCCAACGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTT                    >  SRR3722125.573349/1‑101 (MQ=60)
                                                                        CGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTA              <  SRR3722125.779716/101‑1 (MQ=60)
                                                                        CGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTA              <  SRR3722125.849748/101‑1 (MQ=60)
                                                                        CGCCGTGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTA              <  SRR3722125.986805/101‑1 (MQ=60)
                                                                             TGCCTTGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCctgtctcttatacacatctccgagcccacgagactaggcatgatctcgtatgccgtct         >  SRR3722125.645745/1‑43 (MQ=60)
                                                                                  TGAGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGA    <  SRR3722125.103769/101‑1 (MQ=60)
                                                                                    AGCGTGTGCGCCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATT  <  SRR3722125.314065/101‑1 (MQ=60)
                                                                                              |                                                                                          
TTCTTAATGAGTTGCAAAAGGATGGGCGTATTTCTAACGTCGAGCTTTCTAAACGTGTGGGACTTTCCCCAACGCCGTGCCTTGAGCGTGTGCGTCGGCTGGAAAGACAAGGGTTTATTCAGGGCTATACGGCGCTGCTTAACCCCCATTATCTGGATGCATCACTTCTGGTATTCGTTGAGATT  >  NZ_CP009273/928106‑928290

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: