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 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 386,045 T→C pseudogene (339/1527 nt) BW25113_RS01915 → outer membrane autotransporter barrel domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273386,0450TC100.0% 218.0 / NA 60pseudogene (339/1527 nt)BW25113_RS01915outer membrane autotransporter barrel domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/26);  total (34/26)

GCACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACACT  >  NZ_CP009273/385960‑386129
                                                                                     |                                                                                    
gCACGATCTCAACCTATTATTTGAACCATGATTATGCGGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGa                                                                                  >  2:268204/1‑90 (MQ=255)
gCACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGa                                                                                  >  2:93995/1‑90 (MQ=255)
gCACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGa                                                                                  >  1:172480/1‑90 (MQ=255)
gCACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGa                                                                                  >  1:287191/1‑90 (MQ=255)
gCACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATCCAGTGATTCACGGCTCGa                                                                                  >  2:377733/1‑90 (MQ=255)
  aCGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAtt                                                                                <  1:252359/90‑1 (MQ=255)
  aCGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAtt                                                                                >  2:282112/1‑90 (MQ=255)
  aCGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAtt                                                                                <  1:520091/90‑1 (MQ=255)
              tattatTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctg                                                                    >  2:106302/1‑90 (MQ=255)
              tattatTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctg                                                                    <  2:112512/90‑1 (MQ=255)
               attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTc                                                                           <  1:404396/82‑1 (MQ=255)
               attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                   >  1:910646/1‑90 (MQ=255)
               attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                   >  2:334739/1‑90 (MQ=255)
               attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                   >  2:435675/1‑90 (MQ=255)
               attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                   >  2:4500/1‑90 (MQ=255)
               attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTAtgctgc                                                                   >  1:274025/1‑90 (MQ=255)
               attattTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGAATATTTc                                                                           >  2:404396/1‑82 (MQ=255)
                    ttGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGc                                                              <  1:457158/90‑1 (MQ=255)
                    ttGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGc                                                              <  1:618066/90‑1 (MQ=255)
                    ttGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGc                                                              <  1:141923/90‑1 (MQ=255)
                     tGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCg                                                             >  1:714933/1‑90 (MQ=255)
                     tGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCg                                                             >  1:785614/1‑90 (MQ=255)
                      gAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGGTTACTTCTATGCTTCCTGGCgg                                                            >  2:510775/1‑90 (MQ=255)
                      gAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCgg                                                            <  1:772697/90‑1 (MQ=255)
                       aaCCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGt                                                           >  1:341364/1‑90 (MQ=255)
                        aCCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGtt                                                          <  1:642527/90‑1 (MQ=255)
                            tGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTAATTTTATGCTGCCTGGCGGttatta                                                      >  2:112115/1‑90 (MQ=255)
                                tATGCGGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGAt                                                  <  1:268204/90‑1 (MQ=255)
                                 aTGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATc                                                 <  2:20189/90‑1 (MQ=255)
                                      gACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTTTTTTTATCGtttt                                            >  2:25278/1‑90 (MQ=255)
                                      gACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGtttt                                            >  1:617496/1‑90 (MQ=255)
                                      gACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGtttt                                            >  2:87893/1‑90 (MQ=255)
                                        cAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGa                                          >  1:356820/1‑90 (MQ=255)
                                          gTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATg                                        >  2:841912/1‑90 (MQ=255)
                                            aCTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCa                                      <  1:106302/90‑1 (MQ=255)
                                               gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATTCTCCCTGGAGGTTATTTAGATCGTTTTTATGCaatt                                   >  2:854417/1‑87 (MQ=255)
                                               gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                   <  2:910646/90‑1 (MQ=255)
                                               gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                   <  1:745766/90‑1 (MQ=255)
                                               gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                   <  1:518807/90‑1 (MQ=255)
                                               gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                   <  2:785614/90‑1 (MQ=255)
                                               gCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGGTATTATGATCGTTTTGATGCAGAt                                   >  2:7111/1‑90 (MQ=255)
                                                  aaTCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTTTTTTGATCGTTTTGATGCAGATTGt                                >  2:595584/1‑90 (MQ=255)
                                                    tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTTTGATCGTTTTGATGCAGATGGtaa                              >  2:560891/1‑90 (MQ=255)
                                                    tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                   <  1:804619/85‑1 (MQ=255)
                                                    tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGAt                                   >  2:804619/1‑85 (MQ=255)
                                                    tCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGtaa                              <  1:556822/90‑1 (MQ=255)
                                                      aGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGGTATTATGATCGTTTTGATGGAGATGGgaata                            >  1:818984/1‑90 (MQ=255)
                                                                 aGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGGAATAATCTggggggg                 >  2:875178/1‑89 (MQ=255)
                                                                   tAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGata               <  1:103585/90‑1 (MQ=255)
                                                                     aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg             <  1:32543/90‑1 (MQ=255)
                                                                     aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg             <  2:743978/90‑1 (MQ=255)
                                                                     aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg             <  2:223654/90‑1 (MQ=255)
                                                                     aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg             <  2:215874/90‑1 (MQ=255)
                                                                     aTTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATg             <  1:875178/90‑1 (MQ=255)
                                                                          gTGATTCACGGCTCGATTACTTCTTTGCTGCCTGGGGGTTATTTTGATCGGTTTGATGCAGATGGTAATAATCTGGGTGGATATGAtttt        >  2:40812/1‑90 (MQ=255)
                                                                          gTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTTATGCAGATGGTAATAATCTGGGGGGATATGAtttt        >  2:911517/1‑90 (MQ=255)
                                                                            gATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGGGGATATGATTTTTa      >  1:77139/1‑90 (MQ=255)
                                                                             aTTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTac     <  2:39567/90‑1 (MQ=255)
                                                                              ttCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTaca    <  2:483579/90‑1 (MQ=255)
                                                                                cACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACACt  >  1:308407/1‑90 (MQ=255)
                                                                                     |                                                                                    
GCACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACACT  >  NZ_CP009273/385960‑386129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ 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
ACCGTTAACGGCACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACACTGATGCGGTTGTTGAT  >  NZ_CP009273/385950‑386144
                                                                                               |                                                                                                   
ACCGTTAACGGCACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGA                                                                                                 >  SRR3722072.175466/1‑100 (MQ=60)
ACCGTTAACGGCACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGA                                                                                                 >  SRR3722072.292139/1‑100 (MQ=60)
            ACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTGTGC                                                                                     <  SRR3722072.529447/100‑1 (MQ=60)
            ACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGC                                                                                     <  SRR3722072.256631/100‑1 (MQ=60)
               ATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGC                                                                                  >  SRR3722072.926303/1‑100 (MQ=60)
               ATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGC                                                                                  >  SRR3722072.278714/1‑100 (MQ=60)
                 CTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCT                                                                                <  SRR3722072.411690/100‑1 (MQ=60)
                     ACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCG                                                                            >  SRR3722072.799378/1‑100 (MQ=60)
                     ACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCG                                                                            >  SRR3722072.727665/1‑100 (MQ=60)
                      CCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCctgtctctta                                                                           >  SRR3722072.14959/1‑90 (MQ=60)
                       CTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGT                                                                          >  SRR3722072.347377/1‑100 (MQ=60)
                              TTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATG                                                                   <  SRR3722072.144460/100‑1 (MQ=60)
                              TTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATG                                                                   <  SRR3722072.629091/100‑1 (MQ=60)
                              TTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATG                                                                   <  SRR3722072.465525/100‑1 (MQ=60)
                                GAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGAT                                                                 <  SRR3722072.786272/100‑1 (MQ=60)
                                  ACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCG                                                               <  SRR3722072.653984/100‑1 (MQ=60)
                                      TGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTT                                                           >  SRR3722072.628509/1‑100 (MQ=60)
                                        ATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGA                                                         >  SRR3722072.363133/1‑100 (MQ=60)
                                          TATGCGGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATG                                                       <  SRR3722072.272768/100‑1 (MQ=60)
                                                      ACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATA                                           <  SRR3722072.108260/100‑1 (MQ=60)
                                                      ACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGGTATTATGATCGTTTTGATGGAGATGGGAATA                                           >  SRR3722072.833246/1‑100 (MQ=60)
                                                         GCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATC                                        <  SRR3722072.528142/100‑1 (MQ=60)
                                                         GCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATC                                        <  SRR3722072.758954/100‑1 (MQ=60)
                                                         GCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATC                                        <  SRR3722072.818661/100‑1 (MQ=60)
                                                              TCAGCTTGATATCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGT                                   <  SRR3722072.566829/100‑1 (MQ=60)
                                                                         TCAGTAATTCAGTGATTCACGGCTCGATTACTTCTATGCctgtctcttatacacatctccgagcccacgagacggccgattttctcgtatgccggcttct                        >  SRR3722072.383175/1‑39 (MQ=60)
                                                                            GTAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGGGGATATGATTTTTA                     >  SRR3722072.78571/1‑100 (MQ=60)
                                                                             TAATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTAC                    <  SRR3722072.105495/100‑1 (MQ=60)
                                                                               ATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACAC                  <  SRR3722072.33122/100‑1 (MQ=60)
                                                                               ATTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACAC                  <  SRR3722072.890250/100‑1 (MQ=60)
                                                                                TTCAGTGATTCACGGCTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACACT                 >  SRR3722072.313757/1‑100 (MQ=60)
                                                                                               CTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACACTGATGCGGTTGTTGAT  <  SRR3722072.570954/100‑1 (MQ=60)
                                                                                               |                                                                                                   
ACCGTTAACGGCACGATCTCAACCTATTATTTGAACCATGATTATGCAGACAGTACTGCTAATCAGCTTGATATCAGTAATTCAGTGATTCACGGTTCGATTACTTCTATGCTGCCTGGCGGTTATTATGATCGTTTTGATGCAGATGGTAATAATCTGGGTGGATATGATTTTTACACTGATGCGGTTGTTGAT  >  NZ_CP009273/385950‑386144

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