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 3,898,828 (T)5→4 coding (130/1878 nt) bglF ← PTS beta‑glucoside transporter subunit IIABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,898,8240T.98.2% 224.5 / ‑2.7 54coding (134/1878 nt)bglFPTS beta‑glucoside transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (35/18);  minor base C (0/1);  total (35/19)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG  >  NZ_CP009273/3898735‑3898907
                                                                                         |                                                                                   
tAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCt                                                                                     >  1:523073/1‑90 (MQ=255)
tAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCt                                                                                     <  2:697427/90‑1 (MQ=255)
tAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCt                                                                                     <  2:602433/90‑1 (MQ=255)
 aaCCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCtt                                                                                    >  2:528750/1‑90 (MQ=255)
    cGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttc                                                                                 <  2:216208/90‑2 (MQ=255)
                aTCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCt                                                                    >  2:353688/1‑90 (MQ=255)
                 tCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCtt                                                                   >  1:11512/1‑90 (MQ=255)
                  cGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTg                                                                  <  1:506890/90‑1 (MQ=255)
                      cacaTGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCgct                                                              >  1:441162/1‑90 (MQ=255)
                              tACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTc                                                      >  1:367721/1‑90 (MQ=255)
                               aCCTATGACCACCTGAATCTGGCCACCGCTTGCCCCCCCCATAACAAAACCGGGGGCCTTTTCAGTACCTCTTCCTGCGCTTTGCTTTCa                                                     <  1:20612/90‑1 (MQ=255)
                                ccTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCAt                                                    <  1:480866/90‑1 (MQ=255)
                                  tATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc                                                  >  2:429500/1‑90 (MQ=255)
                                  tATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc                                                  >  1:349103/1‑90 (MQ=255)
                                  tATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATcc                                                  >  2:208391/1‑90 (MQ=255)
                                    tGACCACCTTAAACTGGCCACCGCTTTCCCCCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtt                                                <  1:212078/90‑1 (MQ=255)
                                    tGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcag                                                                               <  2:552769/60‑4 (MQ=255)
                                    tGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCttttcag                                                                               >  1:552769/1‑57 (MQ=255)
                                          ccTGAAACTGGCCACCGCTTTCCCCCCCCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaattt                                          <  1:560205/90‑1 (MQ=255)
                                          ccTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaattt                                          >  2:557186/1‑90 (MQ=255)
                                            tGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCtttaatttaa                                        <  1:755173/90‑1 (MQ=255)
                                                   ggCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTaa                                 >  2:644449/1‑90 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAATACGTAAt                                >  2:389861/1‑90 (MQ=255)
                                                    gCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                <  2:214108/90‑1 (MQ=255)
                                                       aCCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCCTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATcgc                             >  1:371908/1‑90 (MQ=255)
                                                             ttCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcg                       >  2:501000/1‑90 (MQ=255)
                                                             ttCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcg                       >  2:758048/1‑90 (MQ=255)
                                                              tCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcgc                      <  2:335678/90‑1 (MQ=255)
                                                                caccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                    <  2:198672/90‑1 (MQ=255)
                                                                caccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCaa                    <  2:687617/90‑1 (MQ=255)
                                                                  ccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATg                  >  1:642571/1‑90 (MQ=255)
                                                                  ccaccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATg                  >  1:445148/1‑90 (MQ=255)
                                                                     ccaTAATAATACTGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAt               >  1:38717/1‑90 (MQ=255)
                                                                     ccaTAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAt               >  1:282836/1‑90 (MQ=255)
                                                                       ataataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                <  1:478179/71‑1 (MQ=255)
                                                                       ataataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAt                                >  2:478179/1‑71 (MQ=255)
                                                                       ataataataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCa             >  2:126309/1‑90 (MQ=255)
                                                                             ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca       <  2:99797/90‑1 (MQ=255)
                                                                             ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca       <  2:65102/90‑1 (MQ=255)
                                                                             ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca       <  1:398296/90‑1 (MQ=255)
                                                                             ataCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTca       >  1:436501/1‑90 (MQ=255)
                                                                              taCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcac      <  2:164225/90‑1 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca     >  1:106900/1‑90 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca     >  1:461601/1‑90 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca     >  2:481398/1‑90 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcaca     >  1:648199/1‑90 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcac      <  2:721822/89‑1 (MQ=255)
                                                                               aCCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTcac      >  1:721822/1‑89 (MQ=255)
                                                                                 cGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAAt   >  1:570104/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAAc                                     <  2:147713/55‑1 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAAc                                     >  1:147713/1‑55 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcgc                      >  1:309936/1‑70 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTgcgc                      <  2:309936/70‑1 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:213286/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  <  2:461601/90‑1 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  1:85262/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:33521/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  1:27466/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  1:225466/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:711233/1‑90 (MQ=255)
                                                                                  gggggTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATg  >  2:302990/1‑90 (MQ=255)
                                                                                         |                                                                                   
TAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG  >  NZ_CP009273/3898735‑3898907

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ 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
TGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT  >  NZ_CP009273/3898732‑3898925
                                                                                                |                                                                                                 
TGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGT                                                                                               <  SRR3722073.380615/100‑1 (MQ=60)
  TTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTAC                                                                                             <  SRR3722073.209065/100‑1 (MQ=60)
          CAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTT                                                                                     >  SRR3722073.11806/1‑100 (MQ=60)
               AGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCT                                                                                >  SRR3722073.450592/1‑100 (MQ=60)
                     CGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTT                                                                          <  SRR3722073.517588/100‑1 (MQ=60)
                       GCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTC                                                                        >  SRR3722073.375527/1‑100 (MQ=60)
                           CATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCC                                                                    >  SRR3722073.356497/1‑100 (MQ=60)
                             TGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTctgtctcttatacacatctg                                                                  >  SRR3722073.564320/1‑80 (MQ=60)
                                  ACCTATGACCACCTGAATCTGGCCACCGCTTGCCCCCCCCATAACAAAACCGGGGGCCTTTTCAGTACCTCTTCCTGCGCTTTGCTTTCATCCTGTAATT                                                             <  SRR3722073.21136/100‑1 (MQ=37)
                                   CCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTT                                                            <  SRR3722073.491090/100‑1 (MQ=60)
                                       TGACCACCTTAAACTGGCCACCGCTTTCCCCCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAA                                                        <  SRR3722073.216443/100‑1 (MQ=60)
                                             CCTGAAACTGGCCACCGCTTTCCCCCCCCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAAT                                                  <  SRR3722073.571895/100‑1 (MQ=60)
                                               TGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCG                                                <  SRR3722073.770445/100‑1 (MQ=60)
                                                GAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCCTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGC                                               >  SRR3722073.379792/1‑100 (MQ=60)
                                                 AAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTctgtctcttatacac                                              <  SRR3722073.688280/100‑16 (MQ=60)
                                                       aagagacagTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGC                                        <  SRR3722073.488354/91‑1 (MQ=60)
                                                           CCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATG                                    >  SRR3722073.655823/1‑100 (MQ=60)
                                                           CCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATG                                    >  SRR3722073.454658/1‑100 (MQ=60)
                                                              CTTTCCACCACCATAATAATACTGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAT                                 >  SRR3722073.39663/1‑100 (MQ=60)
                                                              CTTTCCACCACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCAT                                 >  SRR3722073.288824/1‑100 (MQ=60)
                                                                      CACCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCA                         >  SRR3722073.445838/1‑100 (MQ=60)
                                                                        CCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACA                       >  SRR3722073.736405/1‑100 (MQ=60)
                                                                        CCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACA                       >  SRR3722073.661543/1‑100 (MQ=60)
                                                                        CCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACA                       >  SRR3722073.471449/1‑100 (MQ=60)
                                                                        CCATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACA                       >  SRR3722073.109210/1‑100 (MQ=60)
                                                                          ATAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAAT                     >  SRR3722073.582019/1‑100 (MQ=60)
                                                                           TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTctgtctcttatacacatctgacgct                    >  SRR3722073.150756/1‑75 (MQ=60)
                                                                           TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGctgtctctta                    >  SRR3722073.316502/1‑90 (MQ=60)
                                                                           TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                    >  SRR3722073.28165/1‑100 (MQ=60)
                                                                           TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                    >  SRR3722073.230169/1‑100 (MQ=60)
                                                                           TAATAATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATG                    >  SRR3722073.87184/1‑100 (MQ=60)
                                                                                ATACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATC               <  SRR3722073.406780/100‑1 (MQ=60)
                                                                                 TACCGGGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCT              >  SRR3722073.19009/1‑100 (MQ=60)
                                                                                      GGGGTCTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCC         >  SRR3722073.532905/1‑100 (MQ=60)
                                                                                              TTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCctgtctcttatac  >  SRR3722073.731826/1‑87 (MQ=60)
                                                                                              TTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT  >  SRR3722073.74651/1‑100 (MQ=60)
                                                                                                |                                                                                                 
TGTTAACCGCCAGGAAGACATCGGCCACATGGTTACCTATGACCACCTGAAACTGGCCACCGCTTTCCACCACCATAATAATACCGGGGGTCTTTTTCAGTACCTCTGCTTGCGCTTTGCTTTCATCCTTTAATTTAAAACGTAATCGCGTTGCGCAATGCATCAGACTCACAATGTTATCTGCGCCCCCGACT  >  NZ_CP009273/3898732‑3898925

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: