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 3,475,701 T→C S352P (TCG→CCG)  yheS → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,475,7010TC100.0% 200.2 / NA 60S352P (TCG→CCG) yheSABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/30);  total (30/30)

GATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTT  >  NZ_CP009273/3475617‑3475786
                                                                                    |                                                                                     
gATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACa                                                                                  <  1:1084582/90‑1 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 >  1:632232/1‑90 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 >  1:794303/1‑90 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 <  2:1081658/90‑1 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 >  1:283401/1‑90 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 >  2:1506543/1‑90 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 >  1:1461977/1‑90 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAt                                                                                 >  1:1058309/1‑90 (MQ=255)
 aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGGAAACCGACAt                                                                                 >  1:921911/1‑90 (MQ=255)
      attattCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATc                                                                            <  2:591262/90‑1 (MQ=255)
        tattCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCaa                                                                          <  2:641230/90‑1 (MQ=255)
        tattCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCaa                                                                          <  2:486225/90‑1 (MQ=255)
        tattCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCaa                                                                          <  2:1511903/90‑1 (MQ=255)
              cgactcgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGtt                                                                    >  1:1163293/1‑90 (MQ=255)
               gactcgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTa                                                                   >  2:697814/1‑90 (MQ=255)
                  tcgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGcc                                                                <  2:903267/90‑1 (MQ=255)
                   cgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCaa                                                                          <  2:1157711/79‑1 (MQ=255)
                   cgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCaa                                                                          >  1:1157711/1‑79 (MQ=255)
                    gaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCgg                                                              >  2:987995/1‑90 (MQ=255)
                      ttAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTg                                                            <  1:87622/90‑1 (MQ=255)
                          aCTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACt                                                        <  2:1163293/90‑1 (MQ=255)
                          aCTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACt                                                        <  1:987995/90‑1 (MQ=255)
                          aCTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACt                                                        <  2:1058596/90‑1 (MQ=255)
                                  tGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAg                                                >  1:556488/1‑90 (MQ=255)
                                   ggTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGt                                               >  2:977011/1‑90 (MQ=255)
                                    gTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTc                                              <  2:209154/90‑1 (MQ=255)
                                    gTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTc                                              <  2:579088/90‑1 (MQ=255)
                                     tGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCa                                             >  2:1571104/1‑90 (MQ=255)
                                       cccGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGc                                           <  1:977011/90‑1 (MQ=255)
                                        ccGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCg                                          >  1:1523168/1‑90 (MQ=255)
                                         cGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCgg                                         >  1:137572/1‑90 (MQ=255)
                                          ggCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGt                                        >  1:170051/1‑90 (MQ=255)
                                                    ttGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTc                              <  1:1198053/90‑1 (MQ=255)
                                                    ttGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTc                              <  2:1523168/90‑1 (MQ=255)
                                                    ttGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTc                              <  2:1266144/90‑1 (MQ=255)
                                                        tCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGc                          >  2:842368/1‑90 (MQ=255)
                                                         cTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCg                         >  2:287892/1‑90 (MQ=255)
                                                           gTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGaa                       >  1:493465/1‑90 (MQ=255)
                                                            ttAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGaaa                      <  1:202521/90‑1 (MQ=255)
                                                             tAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAg                     <  2:1175298/90‑1 (MQ=255)
                                                             tAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAg                     <  2:357196/90‑1 (MQ=255)
                                                             tAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAg                     <  1:1506543/90‑1 (MQ=255)
                                                                gCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGa                  >  1:399294/1‑90 (MQ=255)
                                                                gCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGGGAAATTGGTCTGGCGAAAGGGa                  >  1:921829/1‑90 (MQ=255)
                                                                 ccGCAATGGCGCGGGTAAACCGACATTAATCAAAc                                                                        >  1:207413/1‑35 (MQ=38)
                                                                 ccGCAATGGCGCGGGTAAACCGACATTAATCAAAc                                                                        <  2:207413/35‑1 (MQ=38)
                                                                 ccGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGt                                                             <  1:307228/46‑1 (MQ=255)
                                                                 ccGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGt                                                             >  2:307228/1‑46 (MQ=255)
                                                                 ccGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGAt                 <  2:24275/90‑1 (MQ=255)
                                                                 ccGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGAt                 >  1:1498729/1‑90 (MQ=255)
                                                                   gCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCa               >  1:1274739/1‑90 (MQ=255)
                                                                      aTGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAAc            <  1:441843/90‑1 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGAt                 <  2:376789/82‑1 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGAt                 >  1:376789/1‑82 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCg         >  1:622685/1‑90 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCg         >  2:1070674/1‑90 (MQ=255)
                                                                         gcgcGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCg         >  1:458560/1‑90 (MQ=255)
                                                                           gcgGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGc       <  2:921829/90‑1 (MQ=255)
                                                                                tAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACtt  <  1:468915/90‑1 (MQ=255)
                                                                                tAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACtt  <  2:1461977/90‑1 (MQ=255)
                                                                                    |                                                                                     
GATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTT  >  NZ_CP009273/3475617‑3475786

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
CGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTTCGCCCAGCATCA  >  NZ_CP009273/3475604‑3475798
                                                                                                 |                                                                                                 
CGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCG                                                                                                 <  SRR3722078.751511/100‑1 (MQ=60)
    GGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAT                                                                                             >  SRR3722078.1074218/1‑100 (MQ=60)
    GGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAT                                                                                             >  SRR3722078.804713/1‑100 (MQ=60)
    GGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAT                                                                                             >  SRR3722078.1485832/1‑100 (MQ=60)
    GGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAT                                                                                             >  SRR3722078.639503/1‑100 (MQ=60)
    GGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACAT                                                                                             >  SRR3722078.286863/1‑100 (MQ=60)
    GGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGGAAACCGACAT                                                                                             >  SRR3722078.934523/1‑100 (MQ=60)
         TGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATC                                                                                        <  SRR3722078.1251755/100‑1 (MQ=60)
             GATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAAC                                                                                    <  SRR3722078.1101001/100‑1 (MQ=60)
                 GCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTT                                                                                >  SRR3722078.1181338/1‑100 (MQ=60)
                      ATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCC                                                                           >  SRR3722078.1175620/1‑100 (MQ=60)
                                  atcgtcggcagcgtcagatgtgtataagagacaGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGC                                                               <  SRR3722078.310978/67‑1 (MQ=60)
                                   TTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCC                                                              <  SRR3722078.88656/100‑1 (MQ=60)
                                     AAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAG                                                            >  SRR3722078.563047/1‑100 (MQ=60)
                                       ACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTC                                                          <  SRR3722078.1002316/100‑1 (MQ=60)
                                           AACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCG                                                      >  SRR3722078.1547999/1‑100 (MQ=60)
                                            ACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGG                                                     >  SRR3722078.139170/1‑100 (MQ=60)
                                             CCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGT                                                    >  SRR3722078.172033/1‑100 (MQ=60)
                                                    CCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTG                                             <  SRR3722078.991061/100‑1 (MQ=60)
                                                              GTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAA                                   >  SRR3722078.499311/1‑100 (MQ=60)
                                                                 TTGGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGG                                <  SRR3722078.1216845/100‑1 (MQ=60)
                                                                   GGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGA                              >  SRR3722078.404133/1‑100 (MQ=60)
                                                                   GGTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGGGAAATTGGTCTGGCGAAAGGGA                              >  SRR3722078.934440/1‑100 (MQ=60)
                                                                    GTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGctgtctcttatacacatctgacgctgccgacgatggcggtcgtgt                             >  SRR3722078.209921/1‑55 (MQ=60)
                                                                    GTCTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGAT                             >  SRR3722078.1523173/1‑100 (MQ=60)
                                                                      CTGTTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCA                           >  SRR3722078.1294988/1‑100 (MQ=60)
                                                                         TTAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAAC                        <  SRR3722078.204961/100‑1 (MQ=60)
                                                                          TAGGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACT                       <  SRR3722078.1531109/100‑1 (MQ=60)
                                                                            GGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCG                     >  SRR3722078.464036/1‑100 (MQ=60)
                                                                            GGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCG                     >  SRR3722078.381390/1‑100 (MQ=60)
                                                                            GGCCGCAATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCG                     >  SRR3722078.629861/1‑100 (MQ=60)
                                                                                   ATGGCGCGGGTAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTT              <  SRR3722078.447140/100‑1 (MQ=60)
                                                                                             TAAACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTTCGCCCAGCAT    <  SRR3722078.474493/100‑1 (MQ=60)
                                                                                               AACCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTTCGCCCAGCATCA  >  SRR3722078.551978/1‑100 (MQ=60)
                                                                                                 |                                                                                                 
CGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGGCTACTTCGCCCAGCATCA  >  NZ_CP009273/3475604‑3475798

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