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 F3 I149 R1 8 65.0 3813846 96.1% 3665106 85.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,267,454 C→T R20C (CGC→TGC)  soxR → redox‑sensitive transcriptional activator SoxR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,267,4540CT98.4% 237.6 / ‑3.9 62R20C (CGC→TGC) soxRredox‑sensitive transcriptional activator SoxR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (24/37);  total (24/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

TCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCG  >  NZ_CP009273/4267370‑4267539
                                                                                    |                                                                                     
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  >  1:1061041/1‑90 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  >  2:661340/1‑90 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  <  2:1831952/90‑1 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  >  1:1362925/1‑90 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  <  2:1128150/90‑1 (MQ=255)
tCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  >  1:1510968/1‑90 (MQ=255)
       aaCTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTgg                                                                           <  2:1535764/90‑1 (MQ=255)
        aCTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGc                                                                          >  1:807032/1‑90 (MQ=255)
         cTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCg                                                                         <  1:1877216/90‑1 (MQ=255)
           tGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGt                                                                       >  2:213664/1‑90 (MQ=255)
            gAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTa                                                                      <  1:1596196/90‑1 (MQ=255)
              ggTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATc                                                                    <  2:966009/90‑1 (MQ=255)
                 aaaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgc                                                                 >  1:768171/1‑90 (MQ=255)
                  aaGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgcg                                                                >  2:579902/1‑90 (MQ=255)
                   aGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgcgc                                                               <  2:1061041/90‑1 (MQ=255)
                   aGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGgcgc                                                               <  1:1479828/90‑1 (MQ=255)
                      gATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGc                                                            >  2:482309/1‑90 (MQ=255)
                        tttATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  >  2:772234/1‑66 (MQ=255)
                        tttATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGc                                                                                  <  1:772234/66‑1 (MQ=255)
                        tttATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCAt                                                          >  2:561134/1‑90 (MQ=255)
                               aaaaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATg                                                   <  1:723555/90‑1 (MQ=255)
                                aaaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGa                                                  <  1:1192497/90‑1 (MQ=255)
                                 aaGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGaa                                                 <  1:215850/90‑1 (MQ=255)
                                         aCCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAgg                                         <  1:1390967/90‑1 (MQ=255)
                                         aCCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAgg                                         <  1:212868/90‑1 (MQ=255)
                                           cccGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                       <  2:1891118/90‑1 (MQ=255)
                                           cccGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                       <  2:733992/90‑1 (MQ=255)
                                           cccGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGt                                       <  2:603027/90‑1 (MQ=255)
                                               cATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGAt                                   <  1:1773663/90‑1 (MQ=255)
                                                aTTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGAtt                                  <  1:482309/90‑1 (MQ=255)
                                                   aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                               <  2:415898/90‑1 (MQ=255)
                                                   aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                               <  2:1095572/90‑1 (MQ=255)
                                                   aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                               <  2:55347/90‑1 (MQ=255)
                                                   aaaGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTAcc                               >  2:633684/1‑90 (MQ=255)
                                                     aGCGCTGCTACCCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                             <  2:202943/90‑1 (MQ=255)
                                                     aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                             <  1:1070996/90‑1 (MQ=255)
                                                     aGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAg                             <  2:1629233/90‑1 (MQ=255)
                                                       cgcTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTa                           <  2:347622/90‑1 (MQ=255)
                                                        gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCAtt                                                         >  1:19534/1‑59 (MQ=255)
                                                        gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCAtt                                                         <  2:19534/59‑1 (MQ=255)
                                                        gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTATGAAAGTAAAGGGGTGATTACCAGTAt                          >  2:1065872/1‑90 (MQ=255)
                                                        gctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTAt                          >  2:1278567/1‑90 (MQ=255)
                                                         ctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCg                                                                         <  1:1043357/42‑1 (MQ=255)
                                                         ctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCg                                                                         >  2:1043357/1‑42 (MQ=255)
                                                         ctgctAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGATTACCAGTATc                         >  2:1259976/1‑90 (MQ=255)
                                                                   ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCCTTTCTATGAAAGTAAAGGGGTGATTACCAGTATCCGTAACAGCg               >  1:876652/1‑90 (MQ=255)
                                                                   ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTTTGAAAGTAAAGGGTTGGTTACCAGTATCCGTAACAGCg               >  2:1832309/1‑90 (MQ=255)
                                                                   ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTATGAAAGGAAAGGGGTGATTACCAGTATCCGTAACAGCg               >  1:763475/1‑90 (MQ=255)
                                                                   ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCg               >  1:258136/1‑90 (MQ=255)
                                                                   ccGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGGTTACCAGTATCCGTAACAGCg               >  1:19462/1‑90 (MQ=255)
                                                                    cGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCgg              <  1:1763336/90‑1 (MQ=255)
                                                                    cGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCgg              <  1:599227/90‑1 (MQ=255)
                                                                     ggCGAAGTGGCGAAATGCAGCGGTGTGGCGGGATCGGCGCCGCATTTCTCTGAAAGGAAAGGGGTGGTTTCCAGTATCCGTAACAGCggg             >  1:1418614/1‑89 (MQ=255)
                                                                             ggCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGaa                                                 <  2:1506337/46‑1 (MQ=255)
                                                                             ggCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGaa                                                 >  1:1506337/1‑46 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc   <  2:480271/90‑1 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc   <  2:1198779/90‑1 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc   <  1:27322/90‑1 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc   <  2:276548/90‑1 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc   <  2:864470/90‑1 (MQ=255)
                                                                               cGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggc   <  2:91278/90‑1 (MQ=255)
                                                                                gAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAgcggcg  <  1:1278567/90‑1 (MQ=255)
                                                                                    |                                                                                     
TCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCG  >  NZ_CP009273/4267370‑4267539

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ 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
GTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATATAAACGTGAT  >  NZ_CP009273/4267360‑4267552
                                                                                              |                                                                                                  
GTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGC                                                                                               >  SRR3722018.1077648/1‑100 (MQ=60)
GTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGC                                                                                               >  SRR3722018.1936650/1‑100 (MQ=60)
GTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGC                                                                                               >  SRR3722018.1537558/1‑100 (MQ=60)
GTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGC                                                                                               >  SRR3722018.1386376/1‑100 (MQ=60)
  ATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGG                                                                                             <  SRR3722018.206016/100‑1 (MQ=60)
  ATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGG                                                                                             <  SRR3722018.1873221/100‑1 (MQ=60)
        CCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGC                                                                                       >  SRR3722018.818528/1‑100 (MQ=60)
          tgtataagagacagGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGG                                                                                     <  SRR3722018.783339/86‑1 (MQ=60)
                 AACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGC                                                                              >  SRR3722018.779227/1‑100 (MQ=60)
                   actcctcgtcggcagcgtcagatgtgtataagagacaGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGC                                                                            <  SRR3722018.1059675/63‑1 (MQ=60)
                   CTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGC                                                                            <  SRR3722018.1910671/100‑1 (MQ=60)
                      GAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGC                                                                         <  SRR3722018.1624444/100‑1 (MQ=60)
                             AGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTA                                                                  <  SRR3722018.1505728/100‑1 (MQ=60)
                                         AAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGG                                                      <  SRR3722018.734064/100‑1 (MQ=60)
                                          AAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGG                                                     <  SRR3722018.1212045/100‑1 (MQ=60)
                                           AAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGT                                                    <  SRR3722018.219009/100‑1 (MQ=60)
                                                   ACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACC                                            <  SRR3722018.1414968/100‑1 (MQ=60)
                                                   ACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACC                                            <  SRR3722018.215985/100‑1 (MQ=60)
                                                        GCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTATGAAAactgtctcttatacacatctga                                       >  SRR3722018.19825/1‑78 (MQ=60)
                                                         CATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATC                                      <  SRR3722018.1805323/100‑1 (MQ=60)
                                                          ATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCC                                     <  SRR3722018.489569/100‑1 (MQ=60)
                                                               AGCGCTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAAC                                <  SRR3722018.1087778/100‑1 (MQ=60)
                                                                   CTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCCTTTCTATGAAAGTAAAGGGGTGATTACCAGTATCCGTAACAGCG                            >  SRR3722018.889639/1‑100 (MQ=60)
                                                                   CTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTTTATGAAAGGAAAGGGGTGATTACCAGTATCCGTAACAGCG                            >  SRR3722018.774488/1‑100 (MQ=60)
                                                                   CTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCG                            >  SRR3722018.261990/1‑100 (MQ=60)
                                                                   CTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGGTGGTTACCAGTATCCGTAACAGCG                            >  SRR3722018.19752/1‑100 (MQ=60)
                                                                   CTGCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGAGGCGGgagcgacgccgcctttccgagaaaggaaagggggaagtgcccgcctacggccacgcg                            >  SRR3722018.526420/1‑43 (MQ=60)
                                                                     GCTAACCCCCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGGATCGGCGCCGCATTTCTCTGAAAGGAAAGGGGTGGTTTCCAGTATCCGTAACAGCGGG                          >  SRR3722018.1443174/1‑100 (MQ=60)
                                                                             CCGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAagggctgtctcttttacacatctgacgctgccgacgag                  >  SRR3722018.1532809/1‑62 (MQ=60)
                                                                              CGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGG                 <  SRR3722018.608134/100‑1 (MQ=60)
                                                                              CGGCGAAGTGGCGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGG                 <  SRR3722018.1794844/100‑1 (MQ=60)
                                                                                         CGAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATATAAACG      <  SRR3722018.27723/100‑1 (MQ=60)
                                                                                          GAAATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATATAAACGT     <  SRR3722018.1300174/100‑1 (MQ=60)
                                                                                             ATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATTAAAGGAAAGGGGTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATATAAACGTGAT  >  SRR3722018.1308775/1‑100 (MQ=60)
                                                                                             ATGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATATAAACGTGAT  >  SRR3722018.484966/1‑100 (MQ=60)
                                                                                              |                                                                                                  
GTATAATTCCTCAAGTTAACTTGAGGTAAAGCGATTTATGGAAAAGAAATTACCCCGCATTAAAGCGCTGCTAACCCCCGGCGAAGTGGCGAAACGCAGCGGTGTGGCGGTATCGGCGCTGCATTTCTATGAAAGTAAAGGGTTGATTACCAGTATCCGTAACAGCGGCAATCAGCGGCGATATAAACGTGAT  >  NZ_CP009273/4267360‑4267552

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