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 475,405 Δ1 bp intergenic (‑31/+141) maa ← / ← hha maltose O‑acetyltransferase/HHA domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273475,4040A.100.0% 262.2 / NA 57intergenic (‑30/+142)maa/hhamaltose O‑acetyltransferase/HHA domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (22/35);  total (22/35)

ACGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCAACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTATGTCTTC  >  NZ_CP009273/475314‑475482
                                                                                          |                                                                              
aCGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCa                                                                                <  1:9979/91‑1 (MQ=255)
 cGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCac                                                                               <  2:738666/91‑2 (MQ=255)
 cGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCac                                                                               <  2:674969/91‑2 (MQ=255)
 cGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCac                                                                               <  2:491753/91‑2 (MQ=255)
 cGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCac                                                                               <  1:218719/91‑2 (MQ=255)
  gTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCaca                                                                              <  2:749912/91‑3 (MQ=255)
   tctcATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCacac                                                                             <  1:100460/91‑4 (MQ=255)
   tctcATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCacac                                                                             <  1:979879/91‑4 (MQ=255)
      cATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGt                                                                         >  1:755494/1‑91 (MQ=255)
         cTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAg                                                                      <  2:770735/91‑1 (MQ=255)
              gAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGaa                                                                 <  2:580271/91‑1 (MQ=255)
                   aTACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTc                                                            <  1:1113026/91‑1 (MQ=255)
                   aTACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTc                                                            >  1:396392/1‑91 (MQ=255)
                     aCAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAg                                                          <  1:832416/91‑1 (MQ=255)
                     aCAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAg                                                          <  1:133320/91‑1 (MQ=255)
                           cACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCac                                                                               >  1:169135/1‑64 (MQ=255)
                           cACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCac                                                                               <  2:169135/65‑2 (MQ=255)
                           cACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAaga                                                    >  1:150499/1‑91 (MQ=255)
                             ccAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTa                                                                        <  1:213597/69‑1 (MQ=255)
                             ccAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTa                                                                        >  2:213597/1‑69 (MQ=255)
                             ccAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAgaga                                                  >  1:271564/1‑91 (MQ=255)
                              cAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGaa                                                 <  1:54863/91‑1 (MQ=255)
                                gCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACg                                               >  2:16687/1‑91 (MQ=255)
                                  aaTCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGtt                                             >  1:182794/1‑91 (MQ=255)
                                   atcatcTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTc                                            >  2:707857/1‑91 (MQ=255)
                                   atcatcTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTc                                            <  1:558975/91‑1 (MQ=255)
                                   atcatcTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTc                                            <  1:469694/91‑1 (MQ=255)
                                   atcatcTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTc                                            <  2:787124/91‑1 (MQ=255)
                                   atcatcTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTc                                            <  2:829807/91‑1 (MQ=255)
                                   atcatcTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTc                                            <  2:966116/91‑1 (MQ=255)
                                        cttttcttttTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCAtt                                       <  1:668805/91‑1 (MQ=255)
                                          tttctttttctGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTAc                                     >  1:698309/1‑91 (MQ=255)
                                             ctttttctGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATc                                  <  1:707857/91‑1 (MQ=255)
                                             ctttttctGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATc                                  <  2:22308/91‑1 (MQ=255)
                                                 ttctGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCtgtg                              <  2:811463/91‑1 (MQ=255)
                                                 ttctGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCtgtg                              <  1:1065400/91‑1 (MQ=255)
                                                      tgCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATa                         <  2:804297/91‑1 (MQ=255)
                                                       gCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTc                                            <  1:65861/71‑1 (MQ=255)
                                                       gCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTc                                            >  2:65861/1‑71 (MQ=255)
                                                        cTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATaaa                       <  1:1100010/91‑1 (MQ=255)
                                                         tCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAg                      <  2:541667/91‑1 (MQ=255)
                                                          cATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGa                     <  2:107637/91‑1 (MQ=255)
                                                          cATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGa                     <  2:912921/91‑1 (MQ=255)
                                                          cATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGa                     >  1:498608/1‑91 (MQ=255)
                                                          cATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGa                     <  2:720363/91‑1 (MQ=255)
                                                          cATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGa                     <  2:707443/91‑1 (MQ=255)
                                                          cATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGa                     >  2:29022/1‑91 (MQ=255)
                                                           aTTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGAt                    <  1:26920/91‑1 (MQ=255)
                                                            ttACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATc                   <  2:734150/91‑1 (MQ=255)
                                                            ttACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATc                   <  2:1095555/91‑1 (MQ=255)
                                                               cTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATcaca                >  2:441583/1‑91 (MQ=255)
                                                                tGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACAt               <  2:320447/91‑1 (MQ=255)
                                                                tGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACAt               <  1:26929/91‑1 (MQ=255)
                                                                 gACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATc                                  >  1:231302/1‑71 (MQ=255)
                                                                 gACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATc                                  <  2:231302/71‑1 (MQ=255)
                                                                 gACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATa              >  1:867281/1‑91 (MQ=255)
                                                                  aCCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAg             <  2:498608/91‑1 (MQ=255)
                                                                   ccTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAgg            <  2:1055609/91‑1 (MQ=255)
                                                                    cTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAg                                                          >  2:921956/1‑44 (MQ=38)
                                                                    cTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAg                                                          <  1:921956/44‑1 (MQ=38)
                                                                    cTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAggg           >  2:58070/1‑91 (MQ=255)
                                                                      ccGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTa         >  2:270289/1‑91 (MQ=255)
                                                                      ccGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTa         >  2:240851/1‑91 (MQ=255)
                                                                        ggTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTATg       >  2:1038102/1‑91 (MQ=255)
                                                                             cTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATc                                  <  1:720529/59‑1 (MQ=255)
                                                                             cTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATc                                  >  2:720529/1‑59 (MQ=255)
                                                                             cTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTATGTCTTc  >  1:726614/1‑91 (MQ=255)
                                                                                          |                                                                              
ACGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCAACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTATGTCTTC  >  NZ_CP009273/475314‑475482

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGCGAGATAACGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCAACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTATGTCTTC  >  NZ_CP009273/475305‑475482
                                                                                                    |                                                                             
CGCGAGATAACGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCAC                                                                               <  SRR3722125.105467/101‑1 (MQ=60)
  ggctcggagatgtgtataagagacaGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACAC                                                                            <  SRR3722125.796152/76‑1 (MQ=60)
     GATAACGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGT                                                                         >  SRR3722125.771836/1‑101 (MQ=60)
       TAACGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAA                                                                       <  SRR3722125.246473/101‑1 (MQ=60)
         ACGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGA                                                                     <  SRR3722125.10170/101‑1 (MQ=60)
          CGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGAC                                                                    <  SRR3722125.223037/101‑1 (MQ=60)
            TCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGA                                                                  <  SRR3722125.102341/101‑1 (MQ=60)
            TCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGA                                                                  <  SRR3722125.1001753/101‑1 (MQ=60)
                tataagagacaGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAG                                                              <  SRR3722125.217816/90‑1 (MQ=60)
                  CTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTC                                                            >  SRR3722125.404711/1‑101 (MQ=60)
                          CGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACctgtctcttatacaca                                                    >  SRR3722125.172397/1‑85 (MQ=60)
                          CGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGA                                                    >  SRR3722125.153357/1‑101 (MQ=60)
                            ATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGA                                                  >  SRR3722125.277039/1‑101 (MQ=60)
                            ATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGA                                                  <  SRR3722125.1138139/101‑1 (MQ=60)
                              ctagtctcgtgggctcggagatgtgtataagagacagCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAAC                                                <  SRR3722125.942420/64‑1 (MQ=60)
                              ACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAAC                                                <  SRR3722125.135838/101‑1 (MQ=60)
                              ACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAAC                                                <  SRR3722125.850709/101‑1 (MQ=60)
                                 ACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTT                                             >  SRR3722125.186346/1‑101 (MQ=60)
                                       CAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATT                                       <  SRR3722125.55873/101‑1 (MQ=60)
                                         GCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTAC                                     >  SRR3722125.713353/1‑101 (MQ=60)
                                            taagagacagCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATC                                  <  SRR3722125.67070/91‑1 (MQ=60)
                                            ATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATC                                  <  SRR3722125.570936/101‑1 (MQ=60)
                                            ATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATC                                  <  SRR3722125.479608/101‑1 (MQ=60)
                                                 CTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGT                             <  SRR3722125.683199/101‑1 (MQ=60)
                                                      cggagatgtgtataagagacagCCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAA                        <  SRR3722125.736068/79‑1 (MQ=60)
                                                      CTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCtcttcactgacttttcgtc                        <  SRR3722125.961845/101‑20 (MQ=60)
                                                      CTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAA                        <  SRR3722125.723122/101‑1 (MQ=60)
                                                         TTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGA                     >  SRR3722125.509152/1‑101 (MQ=60)
                                                          TTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGAT                    <  SRR3722125.1089387/101‑1 (MQ=60)
                                                                GCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAActgtctctta              >  SRR3722125.235927/1‑91 (MQ=60)
                                                                GCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATA              >  SRR3722125.886420/1‑101 (MQ=60)
                                                                 CTCATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAG             <  SRR3722125.1124786/101‑1 (MQ=60)
                                                                    ATTACTGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGT          <  SRR3722125.27393/101‑1 (MQ=60)
                                                                         TGACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTATGTC     <  SRR3722125.27402/101‑1 (MQ=60)
                                                                          GACCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGctgtctcttatacacatctccgagcccacgagactaggcatgatctc    >  SRR3722125.843182/1‑54 (MQ=60)
                                                                            CCTCCGGTTGCTGAAAATGCGGCACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTATGTCTTC  >  SRR3722125.742274/1‑101 (MQ=60)
                                                                                                    |                                                                             
CGCGAGATAACGTCTCATCTGCCGAGCGATACAACTCACCAGCAATCATCTTTTCTTTTTCTGTGCTCATTACTGACCTCCGGTTGCTGAAAATGCGGCAACACAGTAAGACGAAAAGTCAGTGAAGAGAACGTTCCCATTACATCTGTGTGATAAAGATCACATAGGGTATGTCTTC  >  NZ_CP009273/475305‑475482

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