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 F4 I96 R1 227 23.3 882352 98.0% 864704 126.9

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% 146.6 / NA 42D148D (GAT→GACybiXPKHD‑type hydroxylase YbiX
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/22);  total (20/22)

TCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGAT  >  NZ_CP009273/834089‑834353
                                                                                                                                   |                                                                                                                                     
tCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTg                                                                                                                                >  1:7459/1‑139 (MQ=255)
tCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTg                                                                                                                                >  1:274270/1‑139 (MQ=255)
     ttCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCgg                                                                                                                            >  1:197852/1‑138 (MQ=255)
     ttCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCgg                                                                                                                            >  1:54348/1‑138 (MQ=255)
          aCAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGttt                                                                                                                      >  2:435860/1‑139 (MQ=255)
               aTGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTAccc                                                                                                                 >  1:203089/1‑139 (MQ=255)
                  gcgcgcTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAt                                                                                                              <  2:197852/139‑1 (MQ=255)
                   cgcgcTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGAtg                                                                                                             <  2:201463/139‑1 (MQ=255)
                        tttttATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTc                                                                                                        >  2:392403/1‑139 (MQ=255)
                        tttttATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTc                                                                                                        >  1:388849/1‑139 (MQ=255)
                              tcatcGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAg                                                                                                   <  2:54348/138‑1 (MQ=255)
                                  cgcgGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGgtgt                                                                                               >  2:330511/1‑138 (MQ=255)
                                    cgGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCAt                                                                                            >  2:433904/1‑139 (MQ=255)
                                            cGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGAc                                                                                     >  2:146469/1‑138 (MQ=255)
                                                      cacaTAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTcgccg                                                                           <  1:392403/138‑1 (MQ=255)
                                                        caTAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTcgc                                                                             >  2:248582/1‑134 (MQ=255)
                                                        caTAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTcgc                                                                             <  1:248582/134‑1 (MQ=255)
                                                             atgatgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTa                                                                   <  2:272565/139‑1 (MQ=255)
                                                             atgatgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTa                                                                   <  1:157280/139‑1 (MQ=255)
                                                             atgatgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTa                                                                   <  1:433904/139‑1 (MQ=255)
                                                                atgCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCt                                                                <  2:203089/139‑1 (MQ=255)
                                                                  gCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCttt                                                              >  2:245545/1‑139 (MQ=255)
                                                                         gTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGAt                                                        <  2:292302/138‑1 (MQ=255)
                                                                                 cgGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTaaa                                                >  2:6476/1‑138 (MQ=255)
                                                                                        cGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCg                                         <  2:7459/138‑1 (MQ=255)
                                                                                         gggTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTg                                       <  1:146469/139‑1 (MQ=255)
                                                                                              tCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAg                                   >  2:29531/1‑138 (MQ=255)
                                                                                                 cGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAg                               >  1:345411/1‑139 (MQ=255)
                                                                                                  gCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAgg                              <  1:245545/139‑1 (MQ=255)
                                                                                                   cAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGt                             <  1:5555/139‑1 (MQ=255)
                                                                                                   cAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGt                             <  2:404685/139‑1 (MQ=255)
                                                                                                                tgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTa                                                                   <  2:439830/88‑1 (MQ=255)
                                                                                                                tgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTa                                                                   >  1:439830/1‑88 (MQ=255)
                                                                                                                tgGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCATGGCAGAAAGGTCAGTACGCATc                  >  1:352001/1‑137 (MQ=255)
                                                                                                                  gAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAg               <  1:330511/138‑1 (MQ=255)
                                                                                                                      gacacaacaCGCGGTCGCCTGCCGGGAATTTTACCCGATGTTGTCCGAAGACGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTTTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGtt          <  1:401624/136‑1 (MQ=255)
                                                                                                                      gATACATCACGAAGTCGCCTGCCGTGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGtt          <  2:303757/139‑1 (MQ=255)
                                                                                                                             cacGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCg            <  2:375870/130‑1 (MQ=255)
                                                                                                                             cacGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCg            >  1:375870/1‑130 (MQ=255)
                                                                                                                             cacGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGa   <  2:388849/139‑1 (MQ=255)
                                                                                                                             cacGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGa   <  2:274270/139‑1 (MQ=255)
                                                                                                                              acGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGAt  >  1:83607/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                     
TCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGAT  >  NZ_CP009273/834089‑834353

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 35 ≤ 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
GAATATTGTTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGC  >  NZ_CP009273/834078‑834368
                                                                                                                                              |                                                                                                                                                    
GAATATTGTTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGC                                                                                                                                              >  SRR3722119.277114/1‑151 (MQ=60)
GAATATTGTTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGC                                                                                                                                              >  SRR3722119.7515/1‑151 (MQ=60)
     TTGTTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGG                                                                                                                                          >  SRR3722119.199964/1‑150 (MQ=60)
     TTGTTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGG                                                                                                                                          >  SRR3722119.54848/1‑150 (MQ=60)
               GTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCG                                                                                                                               >  SRR3722119.205262/1‑151 (MQ=60)
                        GCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCC                                                                                                                      >  SRR3722119.392749/1‑151 (MQ=60)
                                                             GATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAG                                                                                 <  SRR3722119.251163/151‑1 (MQ=60)
                                                                CCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCT                                                                               <  SRR3722119.396323/150‑1 (MQ=60)
                                                                       AATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGAT                                                                       <  SRR3722119.158888/151‑1 (MQ=60)
                                                                       AATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGAT                                                                       <  SRR3722119.438244/151‑1 (MQ=60)
                                                                          aagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagCAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTGCGCCGCCGTCGTAGCTTTGTGGATCGC                                                                    <  SRR3722119.205655/92‑1 (MQ=60)
                                                                                                 TACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGG                                             >  SRR3722119.348920/1‑151 (MQ=60)
                                                                                                   CGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTC                                           <  SRR3722119.147975/151‑1 (MQ=60)
                                                                                                            CGCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCAT                                  <  SRR3722119.248096/151‑1 (MQ=60)
                                                                                                             GCAATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATC                                 <  SRR3722119.5603/151‑1 (MQ=60)
                                                                                                                ATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCtgtctcttatacacatctccgagcccacgagactaaggcg                              >  SRR3722119.444221/1‑111 (MQ=60)
                                                                                                                ATGCAGGCTGCTGGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCATGGCAGAAAGGTCAGTACGCATCC                                >  SRR3722119.355578/1‑149 (MQ=60)
                                                                                                                            GGAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGG                   <  SRR3722119.333874/150‑1 (MQ=60)
                                                                                                                             GAGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGAT                 >  SRR3722119.379660/1‑151 (MQ=60)
                                                                                                                              AGGGATACAACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATG                >  SRR3722119.84413/1‑151 (MQ=60)
                                                                                                                                GGACACAACACGCGGTCGCCTGCCGGGAATTTTACCCGATGTTGTCCGAAGACGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTTTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGC              <  SRR3722119.405633/151‑1 (MQ=60)
                                                                                                                                       ACACGAGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaacataagatacacgagagcgggt       >  SRR3722119.40081/1‑52 (MQ=60)
                                                                                                                                            AGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGC  <  SRR3722119.440219/151‑1 (MQ=60)
                                                                                                                                            AGGTCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGC  <  SRR3722119.332861/151‑1 (MQ=60)
                                                                                                                                              |                                                                                                                                                    
GAATATTGTTGTCCAGTTCAAACAGCATGGCGCGCTTTTTATCATCGCGGATCATCGACTGGATCCACATAAATGATGCCACTCGTACGCCGCGGGTTACGGGTGTCACGCAATGCAGGCTGCTGGAGGGATACAACACGAGATCGCCTGCCGGGAGTTTTACCCGATGTTGTCCGAAGGTGTCATTAACGACCAGTTCGCCGCCGTCGTAGCTTTGTGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGC  >  NZ_CP009273/834078‑834368

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