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
A2 F5 I175 R1 328 64.7 2582366 96.2% 2484236 120.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,049,216 A→G S807P (TCG→CCG)  torS ← TMAO reductase system sensor histidine kinase/response regulator TorS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,049,2160AG100.0% 144.2 / NA 41S807P (TCG→CCG) torSTMAO reductase system sensor histidine kinase/response regulator TorS
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/21);  total (20/21)

CGCGGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGATTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAACTGGTACGCTGGCGCAGCGTTTCGTCAATGACATG  >  NZ_CP009273/1049093‑1049348
                                                                                                                           |                                                                                                                                    
cgcgGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTTTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGAc                                                                                                                       >  1:145403/1‑139 (MQ=255)
cgcgGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGAc                                                                                                                       >  1:335430/1‑139 (MQ=255)
 gcgGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACt                                                                                                                      <  1:1126036/139‑1 (MQ=255)
   gggCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTg                                                                                                                    <  2:829893/139‑1 (MQ=255)
   gggCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTg                                                                                                                    <  1:281471/139‑1 (MQ=255)
   gggCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTg                                                                                                                    <  1:1262034/139‑1 (MQ=255)
     gCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTTTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCa                                                                                                                  <  2:145403/139‑1 (MQ=255)
                   aTCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGtt                                                                                                                                    <  2:439432/107‑1 (MQ=255)
                   aTCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGtt                                                                                                                                    >  1:439432/1‑107 (MQ=255)
                      gagaAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTgcgc                                                                                                  <  2:511807/138‑1 (MQ=255)
                            cGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAAc                                                                                           <  1:501991/139‑1 (MQ=255)
                                  ggCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGt                                                                                                      >  2:501483/1‑122 (MQ=255)
                                  ggCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGt                                                                                                      <  1:501483/122‑1 (MQ=255)
                                    cATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAg                                                                                   >  2:293479/1‑139 (MQ=255)
                                            gtAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCAc                                                                           >  1:574387/1‑139 (MQ=255)
                                             tAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCAcg                                                                          >  1:415481/1‑139 (MQ=255)
                                                   ataCCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCAc                                                                    <  1:269390/139‑1 (MQ=255)
                                                   ataCCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCAc                                                                    <  2:574387/139‑1 (MQ=255)
                                                   ataCCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCAc                                                                    <  1:33279/139‑1 (MQ=255)
                                                    taCCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGAt                                                                                                                                 <  2:871323/77‑1 (MQ=255)
                                                    taCCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGAt                                                                                                                                 >  1:871323/1‑77 (MQ=255)
                                                    taCCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCAcc                                                                   >  2:240581/1‑139 (MQ=255)
                                                     aCCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTg                                                                                                                    >  1:194842/1‑89 (MQ=255)
                                                     aCCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTg                                                                                                                    <  2:194842/89‑1 (MQ=255)
                                                                         ctcCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCcgcg                                               <  1:240581/138‑1 (MQ=255)
                                                                                cccATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATg                                       <  1:1007941/139‑1 (MQ=255)
                                                                                        cTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTgcgc                                                                                                  >  2:1205901/1‑72 (MQ=255)
                                                                                        cTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTgcgc                                                                                                  <  1:1205901/72‑1 (MQ=255)
                                                                                          gAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTCGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGACCTGGTACg                             >  2:603308/1‑139 (MQ=255)
                                                                                           aGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCAc                                                                    <  2:165265/99‑1 (MQ=255)
                                                                                           aGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCAc                                                                    >  1:165265/1‑99 (MQ=255)
                                                                                           aGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAACTGGTACGc                            >  2:353095/1‑139 (MQ=255)
                                                                                           aGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAACTGGTACGc                            >  2:878766/1‑139 (MQ=255)
                                                                                            gCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTGAAGCGGCACCGGATTAGGGATAATCCCGCGGAATAATGAAGTGGCACGCt                           >  2:960743/1‑139 (MQ=255)
                                                                                                cTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAg                                                                                   >  1:327825/1‑79 (MQ=255)
                                                                                                cTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAg                                                                                   <  2:327825/79‑1 (MQ=255)
                                                                                                          ttGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAACTGGTACGCTGGCGCAGCGTTTCg             >  2:1059702/1‑139 (MQ=255)
                                                                                                             cGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCAc                                                                    >  2:1214377/1‑81 (MQ=255)
                                                                                                             cGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCAc                                                                    <  1:1214377/81‑1 (MQ=255)
                                                                                                                  cATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAACTGGTACGCTGGCGCAGCGTTTCGTCAATGa      >  1:1244601/1‑138 (MQ=255)
                                                                                                                      cAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGAGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAACTGGTACGCTGGCGCAGCGTTTCGTCAATGACATg  <  2:440565/138‑1 (MQ=255)
                                                                                                                           |                                                                                                                                    
CGCGGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGATTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAACTGGTACGCTGGCGCAGCGTTTCGTCAATGACATG  >  NZ_CP009273/1049093‑1049348

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ 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
GCTTTATTTTCTCGCTATCCTGACTGGCGCGGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGATTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAACTGGTACGCTGGCGCAGCGTTTCGTCAATGACATGGGCGCTAAAGCCAATCA  >  NZ_CP009273/1049066‑1049365
                                                                                                                                                      |                                                                                                                                                     
ataagagacagTCGCTATCCTGACTGGCGCGGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGG                                                                                                                                                       <  SRR3722129.662629/140‑1 (MQ=60)
agaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagaCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGG                                                                                                                                                       <  SRR3722129.387164/93‑1 (MQ=60)
  tttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGT                                                                                                                                                      <  SRR3722129.978319/77‑1 (MQ=60)
     ATTTTCTCGCTATCCTGACTGGCGCGGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGAT                                                                                                                                                  <  SRR3722129.33902/151‑1 (MQ=60)
                ATCCTGACTGGCGCGGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTTTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACT                                                                                                                                       >  SRR3722129.147434/1‑151 (MQ=60)
                ATCCTGACTGGCGCGGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACT                                                                                                                                       >  SRR3722129.340267/1‑151 (MQ=60)
                           CGCGGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAG                                                                                                                            <  SRR3722129.1141731/151‑1 (MQ=60)
                             CGGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGAT                                                                                                                          <  SRR3722129.285513/151‑1 (MQ=60)
                             CGGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGAT                                                                                                                          <  SRR3722129.1279575/151‑1 (MQ=60)
                                   ATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGctgtctcttatacacatctccg                                                                                                                    >  SRR3722129.445807/1‑129 (MQ=60)
                                           gagacaGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAAC                                                                                                            <  SRR3722129.508736/145‑1 (MQ=60)
                                               gacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGAC                                                                                                        <  SRR3722129.1222721/94‑1 (MQ=60)
                                                      GCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCAC                                                                                                 <  SRR3722129.509257/151‑1 (MQ=60)
                                                            GGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACG                                                                                           >  SRR3722129.582885/1‑151 (MQ=60)
                                                             GGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGC                                                                                          >  SRR3722129.421537/1‑151 (MQ=60)
                                                                    TGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatg                                                                                   >  SRR3722129.883761/1‑99 (MQ=60)
                                                                     GTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATctgtctcttatacacatctccgagcccacgagactaaggc                                                                                  >  SRR3722129.197559/1‑111 (MQ=60)
                                                                             cgagattcgccttagtctcgtgggctcggagatgtgtataagagacagTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGA                                                                          <  SRR3722129.1231300/103‑1 (MQ=60)
                                                                             AATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGA                                                                          <  SRR3722129.273202/151‑1 (MQ=60)
                                                                             AATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGA                                                                          <  SRR3722129.33688/151‑1 (MQ=60)
                                                                                                agattcgccttagtctcgtgggctcggagatgtgtataagagacagATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGA                                                       <  SRR3722129.431358/105‑1 (MQ=60)
                                                                                                   TCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAAC                                                     <  SRR3722129.243999/150‑1 (MQ=60)
                                                                                                          CCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAACTGGTACGC                                             <  SRR3722129.1022177/151‑1 (MQ=60)
                                                                                                           CCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGActgtctcttatacacatctccgagcccacg                                            >  SRR3722129.167554/1‑121 (MQ=60)
                                                                                                                TAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgt                                        >  SRR3722129.332560/1‑101 (MQ=60)
                                                                                                                                  GAGTTGCGATACATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAACTGGTACGCTGGCGCAGCGTTTCGTCAATGAC                      >  SRR3722129.1261947/1‑150 (MQ=60)
                                                                                                                                             CATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaagaaaaaatgaaaaac          >  SRR3722129.329026/1‑58 (MQ=60)
                                                                                                                                             CATCCAGCGGTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAG                                                                                                                >  SRR3722129.821637/1‑49 (MQ=60)
                                                                                                                                                     GATTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaaacaacatcaaaaccaaaac  >  SRR3722129.254822/1‑50 (MQ=60)
                                                                                                                                                      GTTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAACTGGTACGCTGGCGCAGCGTTTCGTCAATGACATGGGCGCTAAAGCCAATCA  >  SRR3722129.227032/1‑150 (MQ=60)
                                                                                                                                                      |                                                                                                                                                     
GCTTTATTTTCTCGCTATCCTGACTGGCGCGGGCAATGTCGATTTCATCGAGAAGCGGCAGGGCATGTTGTGTAAATAATACCAGCCATTCGTGGATCTTCTCCGTCCCCATTAACTGAGCATCTTCATTGAGTTGCGATACATCCAGCGATTGATCATTATTGACTTGCAGTTGGAGATAGTGCGCCAGTAACTGACCGAGCACTTCACGCGGCACCGGTTTAGGGATAATCCCGCGGAATAATGAACTGGTACGCTGGCGCAGCGTTTCGTCAATGACATGGGCGCTAAAGCCAATCA  >  NZ_CP009273/1049066‑1049365

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