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 I195 R1 222 27.1 1551108 96.0% 1489063 86.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,494,539 T→C F222S (TTC→TCC)  nanX → MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,494,5390TC100.0% 176.6 / NA 51F222S (TTC→TCC) nanXMFS transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/26);  total (25/26)

AGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA  >  NZ_CP009273/4494452‑4494624
                                                                                       |                                                                                     
aGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTccc                                                                                     <  1:111160/90‑1 (MQ=255)
     tGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTc                                                                                <  1:110558/90‑1 (MQ=255)
       gATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTccc                                                                                     <  2:242475/83‑1 (MQ=255)
       gATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTccc                                                                                     >  1:242475/1‑83 (MQ=255)
           gAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGtttt                                                                          >  2:292390/1‑90 (MQ=255)
           gAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGtttt                                                                          >  2:655618/1‑90 (MQ=255)
                   atatAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGCTCGTTTCCCTCGTCTGTTTTTGTCTAtt                                                                  <  1:236598/90‑1 (MQ=255)
                     ataAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTg                                                                >  2:110216/1‑90 (MQ=255)
                     ataAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTg                                                                >  1:341379/1‑90 (MQ=255)
                     ataAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTg                                                                >  1:222667/1‑90 (MQ=255)
                       aaGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGt                                                              >  1:241723/1‑90 (MQ=255)
                       aaGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGt                                                              >  1:297401/1‑90 (MQ=255)
                        aGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTg                                                             >  1:40951/1‑90 (MQ=255)
                           aTAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCaa                                                          >  2:656094/1‑90 (MQ=255)
                            tAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCaaa                                                         >  1:334465/1‑90 (MQ=255)
                            tAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCaaa                                                         >  2:260377/1‑90 (MQ=255)
                             aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                        <  2:632867/90‑1 (MQ=255)
                             aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                        >  2:453445/1‑90 (MQ=255)
                             aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                        <  2:41344/90‑1 (MQ=255)
                             aaaTCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAc                                                        >  2:257269/1‑90 (MQ=255)
                                   aCATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCg                                                  >  1:484631/1‑90 (MQ=255)
                                          tgtctgtcTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACg                                           <  2:241723/90‑1 (MQ=255)
                                          tgtctgtcTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACg                                           <  1:100452/90‑1 (MQ=255)
                                             ctgtctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGac                                        >  2:314992/1‑90 (MQ=255)
                                              tgtctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGact                                       <  1:338589/90‑1 (MQ=255)
                                                tctTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGactac                                     >  2:599799/1‑90 (MQ=255)
                                                   tCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTActtc                                  <  1:644276/90‑1 (MQ=255)
                                                             acaTCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCt                        >  1:55112/1‑90 (MQ=255)
                                                             acaTCTTTCAATCTCTATGATCGTTCCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCt                        >  2:743616/1‑90 (MQ=255)
                                                               aTCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTgg                      <  2:341379/90‑1 (MQ=255)
                                                                tCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGc                     >  1:567008/1‑90 (MQ=255)
                                                                 cTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCa                    <  2:324079/90‑1 (MQ=255)
                                                                     cAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATa                >  1:73365/1‑90 (MQ=255)
                                                                      aaTCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATaa               <  2:297401/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCTCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  1:453053/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  1:292390/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  2:210724/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  1:592294/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  1:580433/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  1:655618/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  1:556541/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  2:629860/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  2:567008/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  1:453445/90‑1 (MQ=255)
                                                                        tctcTATGATCGTTCCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATg             <  1:743616/90‑1 (MQ=255)
                                                                          tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa           <  1:153436/90‑1 (MQ=255)
                                                                          tctATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGa           <  1:599799/90‑1 (MQ=255)
                                                                             ttGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGtt        <  1:197126/89‑1 (MQ=255)
                                                                             aTGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGtt        <  1:314992/90‑1 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa  >  1:457631/1‑90 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa  >  1:479498/1‑90 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa  >  2:703370/1‑90 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACa  >  2:107846/1‑90 (MQ=255)
                                                                                   gTTTCCCTCGTCTGTTTTTGTCTATTTGGGGCAAACTGGCCGATAAACGGGCTACTTCCTTCCTACCTGGCATATAATGGAGTTAATACa  >  1:179321/1‑90 (MQ=255)
                                                                                       |                                                                                     
AGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA  >  NZ_CP009273/4494452‑4494624

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCATTTCAAC  >  NZ_CP009273/4494443‑4494638
                                                                                                |                                                                                                   
CAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCT                                                                                                  <  SRR3722073.333971/100‑1 (MQ=60)
      GTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTG                                                                                            >  SRR3722073.247561/1‑100 (MQ=60)
         AGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTT                                                                                         <  SRR3722073.113566/100‑1 (MQ=60)
              TGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTC                                                                                    <  SRR3722073.112953/100‑1 (MQ=60)
                    GAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTG                                                                              >  SRR3722073.348627/1‑100 (MQ=60)
                    GAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTG                                                                              >  SRR3722073.227286/1‑100 (MQ=60)
                      AGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGT                                                                            >  SRR3722073.303695/1‑100 (MQ=60)
                      AGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGT                                                                            >  SRR3722073.246793/1‑100 (MQ=60)
                       GATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTG                                                                           >  SRR3722073.41944/1‑100 (MQ=60)
                           AATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAA                                                                       >  SRR3722073.341560/1‑100 (MQ=60)
                            ATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGCTCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAAC                                                                      <  SRR3722073.241547/100‑1 (MQ=60)
                                  GGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCG                                                                >  SRR3722073.494927/1‑100 (MQ=60)
                                                   TGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCC                                               <  SRR3722073.102643/100‑1 (MQ=60)
                                                       TGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCC                                           <  SRR3722073.345778/100‑1 (MQ=60)
                                                            TCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCT                                      <  SRR3722073.657556/100‑1 (MQ=60)
                                                            TCAGAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCT                                      >  SRR3722073.56409/1‑100 (MQ=60)
                                                               GAAAACCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGC                                   >  SRR3722073.578841/1‑100 (MQ=60)
                                                                    CCACATCTTTCAATCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATA                              >  SRR3722073.75057/1‑100 (MQ=60)
                                                                                 TCTCTATGATCGTTTCTCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC                 <  SRR3722073.462736/100‑1 (MQ=60)
                                                                                 TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC                 <  SRR3722073.298585/100‑1 (MQ=60)
                                                                                 TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC                 <  SRR3722073.669064/100‑1 (MQ=60)
                                                                                 TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC                 <  SRR3722073.604660/100‑1 (MQ=60)
                                                                                 TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC                 <  SRR3722073.592547/100‑1 (MQ=60)
                                                                                 TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC                 <  SRR3722073.463136/100‑1 (MQ=60)
                                                                                 TCTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC                 <  SRR3722073.568161/100‑1 (MQ=60)
                                                                                 TCTCTATGATCGTTCCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATAC                 <  SRR3722073.758686/100‑1 (MQ=60)
                                                                                  CTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA                >  SRR3722073.489700/1‑100 (MQ=60)
                                                                                  CTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACA                >  SRR3722073.467399/1‑100 (MQ=60)
                                                                                  CTCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGGGCAAACTGGCCGATAAACGGGCTACTTCCTTCCTACCTGGCATATAATGGAGTTAATACA                >  SRR3722073.182788/1‑100 (MQ=60)
                                                                                   TCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAG               <  SRR3722073.156554/100‑1 (MQ=60)
                                                                                   TCTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAG               <  SRR3722073.612311/100‑1 (MQ=60)
                                                                                    CTATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGT              >  SRR3722073.118380/1‑100 (MQ=60)
                                                                                      TTGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGG            <  SRR3722073.201075/100‑1 (MQ=60)
                                                                                      ATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGG            >  SRR3722073.790345/1‑100 (MQ=60)
                                                                                      ATGATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGG            <  SRR3722073.321663/100‑1 (MQ=60)
                                                                                         ATCGTTTCCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCA         >  SRR3722073.241525/1‑100 (MQ=60)
                                                                                                CCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCATTTCAAC  <  SRR3722073.110178/100‑1 (MQ=60)
                                                                                                |                                                                                                   
CAGAAAGTCAGGAGTGGATTGAAGATAAATATAAGGATAAATCAACATTTTTGTCTGTCTTCAGAAAACCACATCTTTCAATCTCTATGATCGTTTTCCTCGTCTGTTTTTGTCTATTTGGTGCAAACTGGCCGATAAACGGACTACTTCCTTCCTACCTGGCAGATAATGGAGTTAATACAGTGGTCATTTCAAC  >  NZ_CP009273/4494443‑4494638

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