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 F5 I97 R1 275 44.0 1733530 96.5% 1672856 122.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 834,220 A→G D148D (GAT→GAC ybiX ← PKHD‑type hydroxylase YbiX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273834,2200AG100.0% 119.6 / NA 35D148D (GAT→GACybiXPKHD‑type hydroxylase YbiX
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/17);  total (18/17)

TTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTG  >  NZ_CP009273/834086‑834350
                                                                                                                                      |                                                                                                                                  
ttgtCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCg                                                                                                                                 <  2:302535/138‑1 (MQ=255)
   tCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTg                                                                                                                             >  1:2243/1‑139 (MQ=255)
        ttCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCgg                                                                                                                         >  1:728576/1‑138 (MQ=255)
        ttCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCgg                                                                                                                         >  1:64507/1‑138 (MQ=255)
            aaCAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGtt                                                                                                                    <  2:2243/139‑1 (MQ=255)
                   tGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTAccc                                                                                                              >  1:75293/1‑138 (MQ=255)
                      cgcgcTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtg                                                                                                          <  2:781976/139‑1 (MQ=255)
                      cgcgcTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtg                                                                                                          >  2:307581/1‑139 (MQ=255)
                                 tcatcGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAg                                                                                                <  2:147339/138‑1 (MQ=255)
                                   atcGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGgt                                                                                              <  1:694581/138‑1 (MQ=255)
                                     cgcgGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGgtgt                                                                                            >  1:116662/1‑138 (MQ=255)
                                               cGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGAc                                                                                  >  2:393770/1‑138 (MQ=255)
                                                    ggATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTc                                                                            >  1:294716/1‑139 (MQ=255)
                                                     gATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTcg                                                                           >  1:301812/1‑139 (MQ=255)
                                                       tCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTcgcc                                                                         <  2:274668/139‑1 (MQ=255)
                                                         cacaTAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTcgccg                                                                        <  2:339189/138‑1 (MQ=255)
                                                               aaTGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCcgtcg                                                                  >  2:364503/1‑138 (MQ=255)
                                                                atgatgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTa                                                                <  2:294716/139‑1 (MQ=255)
                                                                     gCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtg                                                                                                          <  1:266019/92‑1 (MQ=255)
                                                                     gCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtg                                                                                                          >  2:266019/1‑92 (MQ=255)
                                                                     gCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCGGCCGTCGTAGCttt                                                           >  1:433562/1‑139 (MQ=255)
                                                                     gCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCt                                                             >  1:375706/1‑137 (MQ=255)
                                                                     gCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCt                                                             <  2:375706/137‑1 (MQ=255)
                                                                         cTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTtgtg                                                        >  1:670713/1‑138 (MQ=255)
                                                                          tCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGa                                                      >  2:565370/1‑139 (MQ=255)
                                                                                   gcgGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTaa                                              >  2:499894/1‑138 (MQ=255)
                                                                                    cgGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTaaa                                             >  1:602078/1‑138 (MQ=255)
                                                                                           cGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCg                                      <  1:448069/138‑1 (MQ=255)
                                                                                           cGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCg                                      <  2:64507/138‑1 (MQ=255)
                                                                                           cGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCg                                      <  2:728576/138‑1 (MQ=255)
                                                                                                 tCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGa                               <  2:670713/139‑1 (MQ=255)
                                                                                                      cAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGt                          <  1:41740/139‑1 (MQ=255)
                                                                                                                  ctgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATc               >  1:844877/1‑138 (MQ=255)
                                                                                                                      aGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGcc          <  1:364503/139‑1 (MQ=255)
                                                                                                                              aacaCGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTtgtg  <  1:499894/139‑1 (MQ=255)
                                                                                                                                   gaggTCGCCTGCCGGGAGTTTTACCCGATGTTGACCGAAGGTGTCATTAACGAc                                                                                  >  1:425249/1‑54 (MQ=255)
                                                                                                                                   gaggTCGCCTGCCGGGAGTTTTACCCGATGTTGACCGAAGGTGTCATTAACGAc                                                                                  <  2:425249/54‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
TTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTG  >  NZ_CP009273/834086‑834350

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGCGACTGAATATTGTTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGC  >  NZ_CP009273/834071‑834368
                                                                                                                                                     |                                                                                                                                                    
AGCGACTGAATATTGTTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGT                                                                                                                                                     <  SRR3722120.643652/151‑1 (MQ=60)
       GAATATTGTTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGC                                                                                                                                              >  SRR3722120.2275/1‑151 (MQ=60)
            TTGTTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGG                                                                                                                                          >  SRR3722120.738448/1‑150 (MQ=60)
            TTGTTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGG                                                                                                                                          >  SRR3722120.65345/1‑150 (MQ=60)
                       TTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCG                                                                                                                               >  SRR3722120.76282/1‑150 (MQ=60)
                                     tagtctcgtgggctcggagatgtgtataagagacagACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGG                                                                                                                 <  SRR3722120.269839/114‑1 (MQ=60)
                                         CTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTC                                                                                                             >  SRR3722120.118199/1‑150 (MQ=60)
                                                 CATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGA                                                                                                     <  SRR3722120.704063/150‑1 (MQ=60)
                                                        GATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCG                                                                                             >  SRR3722120.298930/1‑151 (MQ=60)
                                                         ATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGC                                                                                            >  SRR3722120.306135/1‑151 (MQ=60)
                                                                         ACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCGGCCGTCGTAGCTTTG                                                                            >  SRR3722120.439913/1‑151 (MQ=60)
                                                                         ACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTG                                                                            >  SRR3722120.381245/1‑151 (MQ=60)
                                                                             AAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGG                                                                         >  SRR3722120.679902/1‑150 (MQ=60)
                                                                                        CTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaca                                                              >  SRR3722120.89242/1‑71 (MQ=60)
                                                                                        CTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAA                                                              >  SRR3722120.610440/1‑150 (MQ=60)
                                                                                                         ACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGG                                             <  SRR3722120.454610/150‑1 (MQ=60)
                                                                                                                    GCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATC                                 <  SRR3722120.42243/151‑1 (MQ=60)
                                                                                                                      AATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCC                                >  SRR3722120.856234/1‑150 (MQ=60)
                                                                                                                                    GAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGAT                 <  SRR3722120.369808/151‑1 (MQ=60)
                                                                                                                                       GGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGACCGAAGGTGTCATTAACGACCAGTTCGCCGCCtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaa               >  SRR3722120.431456/1‑77 (MQ=60)
                                                                                                                                            CAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtat         >  SRR3722120.315277/1‑100 (MQ=60)
                                                                                                                                            CAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGT         <  SRR3722120.507065/151‑1 (MQ=60)
                                                                                                                                                   AGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGC  <  SRR3722120.573276/151‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                                    
AGCGACTGAATATTGTTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGC  >  NZ_CP009273/834071‑834368

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