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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,980,420 A→G G5G (GGA→GGG ftnB → non‑heme ferritin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,980,4200AG100.0% 208.9 / NA 62G5G (GGA→GGGftnBnon‑heme ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/22);  total (40/22)

ATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACA  >  NZ_CP009273/1980333‑1980506
                                                                                       |                                                                                      
attCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGAt                                                                                      >  1:110466/1‑90 (MQ=255)
attCTGCAAACAGCAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGAt                                                                                      >  1:541595/1‑90 (MQ=255)
 ttCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATg                                                                                     <  1:92432/90‑1 (MQ=255)
  tCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGc                                                                                    <  2:109645/90‑1 (MQ=255)
    tGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCtt                                                                                  >  2:719207/1‑90 (MQ=255)
     gCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTtc                                                                                 >  1:1345219/1‑90 (MQ=255)
          cAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCaaa                                                                            >  2:678349/1‑90 (MQ=255)
          cAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCaaa                                                                            >  2:618721/1‑90 (MQ=255)
            gTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAaact                                                                          >  2:646215/1‑90 (MQ=255)
                    ggTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACtctc                                                                  >  1:1080300/1‑90 (MQ=255)
                    ggTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACtctc                                                                  >  2:1500318/1‑90 (MQ=255)
                    ggTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACtctc                                                                  >  1:904840/1‑90 (MQ=255)
                      tgtTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCaa                                                                >  1:957864/1‑90 (MQ=255)
                       gtTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCaaa                                                               >  2:597379/1‑90 (MQ=255)
                         tttGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATg                                                             <  2:945402/90‑1 (MQ=255)
                               ttATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACcgc                                                       >  1:651709/1‑90 (MQ=255)
                                      gCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGtttt                                                >  1:761943/1‑90 (MQ=255)
                                      gCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGtttt                                                >  2:1413422/1‑90 (MQ=255)
                                          ctctCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGc                                            <  1:880023/90‑1 (MQ=255)
                                           tctcCACTTCTGGATATAAGGGTATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCa                                           >  1:897364/1‑90 (MQ=255)
                                           tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCa                                           >  2:674306/1‑90 (MQ=255)
                                           tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCa                                           >  2:329633/1‑90 (MQ=255)
                                           tctcCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCa                                           >  2:645033/1‑90 (MQ=255)
                                              ccACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATcc                                        <  1:143818/90‑1 (MQ=255)
                                               cACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCa                                       <  2:432757/90‑1 (MQ=255)
                                                aCTTCTGGATATAAGGGTATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCaa                                      <  2:897364/90‑1 (MQ=255)
                                                aCTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCaa                                      <  2:1540304/90‑1 (MQ=255)
                                                   tCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtct                                   >  2:449822/1‑90 (MQ=255)
                                                    cTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                                  <  1:265630/90‑1 (MQ=255)
                                                      ggATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTa                                >  2:1207971/1‑90 (MQ=255)
                                                           taAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTc                           >  1:853866/1‑90 (MQ=255)
                                                            aaGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCa                          <  2:718522/90‑1 (MQ=255)
                                                             aGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCAc                         >  2:431116/1‑90 (MQ=255)
                                                               gATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCt                       <  1:719207/90‑1 (MQ=255)
                                                                atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                      >  2:104657/1‑90 (MQ=255)
                                                                atatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATAAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTg                      >  2:1564227/1‑90 (MQ=255)
                                                                 tatTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGa                     >  2:349775/1‑90 (MQ=255)
                                                                   ttAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGa                     <  1:426809/88‑1 (MQ=255)
                                                                   ttAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGa                     >  2:426809/1‑88 (MQ=255)
                                                                      ggTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAc                <  2:761943/90‑1 (MQ=255)
                                                                      ggTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAc                >  1:361920/1‑90 (MQ=255)
                                                                       gTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACt               >  2:584242/1‑90 (MQ=255)
                                                                        tATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTg              >  2:1084304/1‑90 (MQ=255)
                                                                         aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             <  1:1122709/90‑1 (MQ=255)
                                                                         aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             <  2:350892/90‑1 (MQ=255)
                                                                         aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             >  1:881187/1‑90 (MQ=255)
                                                                         aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             >  1:563259/1‑90 (MQ=255)
                                                                         aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             >  1:409971/1‑90 (MQ=255)
                                                                         aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             <  2:1345219/90‑1 (MQ=255)
                                                                         aTGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTgg             <  2:1576905/90‑1 (MQ=255)
                                                                          tGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgt            >  1:1570824/1‑90 (MQ=255)
                                                                          tGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgt            >  1:903423/1‑90 (MQ=255)
                                                                          tGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgt            >  2:769259/1‑90 (MQ=255)
                                                                          tGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgt            >  2:375277/1‑90 (MQ=255)
                                                                           ggCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGgtg           <  2:1098094/90‑1 (MQ=255)
                                                                            gCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                                  >  2:1455352/1‑66 (MQ=255)
                                                                            gCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAAtctc                                  <  1:1455352/66‑1 (MQ=255)
                                                                               aCCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATcc                                        >  2:460484/1‑57 (MQ=255)
                                                                               aCCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATcc                                        <  1:460484/57‑1 (MQ=255)
                                                                                 cGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGa     >  1:166951/1‑90 (MQ=255)
                                                                                    tGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACa  <  2:881187/90‑1 (MQ=255)
                                                                                    tGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACa  <  1:597379/90‑1 (MQ=255)
                                                                                       |                                                                                      
ATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACA  >  NZ_CP009273/1980333‑1980506

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACG  >  NZ_CP009273/1980322‑1980517
                                                                                                  |                                                                                                 
AGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGA                                                                                                  <  SRR3722078.882700/100‑1 (MQ=60)
AGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGA                                                                                                  <  SRR3722078.712272/100‑1 (MQ=60)
 GAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGAT                                                                                                 >  SRR3722078.111770/1‑100 (MQ=60)
 GAAAGGAATTATTCTGCAAACAGCAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGAT                                                                                                 >  SRR3722078.547991/1‑100 (MQ=60)
      GAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTC                                                                                            >  SRR3722078.1366811/1‑100 (MQ=60)
          TATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAA                                                                                        <  SRR3722078.399321/100‑1 (MQ=60)
            TTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAAC                                                                                      <  SRR3722078.93517/100‑1 (MQ=60)
                     CAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTC                                                                             >  SRR3722078.1096646/1‑100 (MQ=60)
                     CAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTC                                                                             >  SRR3722078.916988/1‑100 (MQ=60)
                       GTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAA                                                                           >  SRR3722078.971422/1‑100 (MQ=60)
                                GTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGC                                                                  >  SRR3722078.659178/1‑100 (MQ=60)
                                       GATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTT                                                           >  SRR3722078.771737/1‑100 (MQ=60)
                                            ATCTTGCACCTCTCCACTTCTGGATATAAGGGTATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCA                                                      >  SRR3722078.909415/1‑100 (MQ=60)
                                                     CTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTC                                             <  SRR3722078.891857/100‑1 (MQ=60)
                                                         cgtcagatgtgtataagagacagAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACC                                         <  SRR3722078.465986/77‑1 (MQ=60)
                                                         CCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACC                                         <  SRR3722078.145485/100‑1 (MQ=60)
                                                            CTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTC                                      >  SRR3722078.865268/1‑100 (MQ=60)
                                                               tgtataagagacagATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACC                                   <  SRR3722078.1479082/86‑1 (MQ=60)
                                                               CTGGATATAAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACC                                   <  SRR3722078.268875/100‑1 (MQ=60)
                                                                       AAGGATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAAC                           >  SRR3722078.366333/1‑100 (MQ=60)
                                                                          GATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACctgtctcttatacacatctgacgctgccgacgatggcggtcgtgtagatctcggtggtcgcc                        >  SRR3722078.904865/1‑38 (MQ=60)
                                                                          GATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGG                        >  SRR3722078.414927/1‑100 (MQ=60)
                                                                          GATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGG                        >  SRR3722078.569870/1‑100 (MQ=60)
                                                                          GATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGG                        >  SRR3722078.893036/1‑100 (MQ=60)
                                                                          GATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGG                        <  SRR3722078.728098/100‑1 (MQ=60)
                                                                           ATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGT                       >  SRR3722078.915539/1‑100 (MQ=60)
                                                                           ATATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGT                       >  SRR3722078.1596383/1‑100 (MQ=60)
                                                                            TATTAGGTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTG                      <  SRR3722078.431947/100‑1 (MQ=60)
                                                                                  GTATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGA                >  SRR3722078.168892/1‑100 (MQ=60)
                                                                                    ATGGCAACCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAAC              <  SRR3722078.1139852/100‑1 (MQ=60)
                                                                                           CCGCTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCT       >  SRR3722078.127029/1‑100 (MQ=60)
                                                                                           CCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCT       >  SRR3722078.1169533/1‑100 (MQ=60)
                                                                                              CTGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTATGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACctgtctctta    >  SRR3722078.107895/1‑90 (MQ=60)
                                                                                               TGGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAAC   <  SRR3722078.604302/100‑1 (MQ=60)
                                                                                                GGGATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACG  <  SRR3722078.105900/100‑1 (MQ=60)
                                                                                                  |                                                                                                 
AGAAAGGAATTATTCTGCAAACAGTAATTATGGTGTTTTGATTTATCTTGCACCTCTCCACTTCTGGATATAAGGATATTAGGTATGGCAACCGCTGGAATGCTTCTCAAACTCAACTCTCAAATGAACCGCGAGTTTTACGCATCCAATCTCTACCTTCACCTGAGTAACTGGTGTTCTGAACAGAGTCTGAACG  >  NZ_CP009273/1980322‑1980517

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