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 543,592 C→T T486I (ACC→ATC)  fdrA → acyl‑CoA synthetase FdrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273543,5920CT94.6% 193.1 / 1.1 56T486I (ACC→ATC) fdrAacyl‑CoA synthetase FdrA
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (24/29);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGC  >  NZ_CP009273/543457‑543718
                                                                                                                                       |                                                                                                                              
cgcgTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCtt                                                                                                                             >  1:60651/1‑139 (MQ=255)
cgcgTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCtt                                                                                                                             >  1:519700/1‑139 (MQ=255)
        gATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTggcg                                                                                                                      <  2:217481/138‑1 (MQ=255)
         aTAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTggcggc                                                                                                                    <  2:956875/139‑1 (MQ=255)
         aTAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTggcggc                                                                                                                    <  1:1277914/139‑1 (MQ=255)
             tCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCg                                                                                                                >  1:1072257/1‑139 (MQ=255)
             tCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCg                                                                                                                >  1:1051893/1‑139 (MQ=255)
             tCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCg                                                                                                                >  1:987764/1‑139 (MQ=255)
                   aCTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAAtt                                                                                                          >  2:1004025/1‑139 (MQ=255)
                       tATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTggcg                                                                                                                      <  2:131830/123‑1 (MQ=255)
                       tATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTggcg                                                                                                                      >  1:131830/1‑123 (MQ=255)
                       tATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATc                                                                                                      >  1:590603/1‑139 (MQ=255)
                       tATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATc                                                                                                      >  2:325996/1‑139 (MQ=255)
                       tATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATc                                                                                                      >  2:682293/1‑139 (MQ=255)
                       tATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATc                                                                                                      >  1:1154853/1‑139 (MQ=255)
                            cATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGctc                                                                                                 <  2:82323/139‑1 (MQ=255)
                            cATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGctc                                                                                                 >  1:451540/1‑139 (MQ=255)
                               tGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCg                                                                                              >  1:1101765/1‑139 (MQ=255)
                               tGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCATGCTGGCGGCAGCGTTAATTCATCCGCTCTCg                                                                                              >  1:356541/1‑139 (MQ=255)
                                 ccACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGcc                                                                                            <  2:590603/139‑1 (MQ=255)
                                 ccACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGcc                                                                                            >  2:837762/1‑139 (MQ=255)
                                      gTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCg                                                                                       <  2:1072257/139‑1 (MQ=255)
                                      gTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCg                                                                                       <  2:519700/139‑1 (MQ=255)
                                         aCAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGca                                                                                     <  2:990245/138‑1 (MQ=255)
                                         aCAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGca                                                                                     <  2:356541/138‑1 (MQ=255)
                                          cAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGcaca                                                                                   <  1:930505/139‑1 (MQ=255)
                                           aGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAg                                                                                  >  2:110730/1‑139 (MQ=255)
                                           aGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAg                                                                                  <  2:1242365/139‑1 (MQ=255)
                                            ggCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGc                                                                                 <  1:1244313/139‑1 (MQ=255)
                                                   aaCGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCg                                                                                       <  2:402786/126‑1 (MQ=255)
                                                   aaCGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCg                                                                                       >  1:402786/1‑126 (MQ=255)
                                                   aaCGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAacaca                                                                           <  2:1101765/138‑1 (MQ=255)
                                                   aaCGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAacaca                                                                           <  2:638937/138‑1 (MQ=255)
                                                         aCCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGt                                                                     <  2:451540/138‑1 (MQ=255)
                                                                      cTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGa                                                           <  2:1113128/135‑1 (MQ=255)
                                                                      cTCGCAGCAAATCGCCACGCTGGAAGATGCGGAGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGa                                                           >  1:1113128/1‑135 (MQ=255)
                                                                          cagcaAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCg                                                   <  1:99920/139‑1 (MQ=255)
                                                                                aTCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGa                                             >  1:223228/1‑139 (MQ=255)
                                                                                       cGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCAt                                                                  <  2:1058342/111‑1 (MQ=255)
                                                                                       cGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCAt                                                                  >  1:1058342/1‑111 (MQ=255)
                                                                                                aTGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACg                             <  2:1173589/139‑1 (MQ=255)
                                                                                                aTGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACg                             <  1:110730/139‑1 (MQ=255)
                                                                                                aTGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACg                             <  1:30025/139‑1 (MQ=255)
                                                                                                     gggATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCt                        <  1:837762/139‑1 (MQ=255)
                                                                                                          tGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCAt                                                                                                       >  2:317284/1‑55 (MQ=255)
                                                                                                          tGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCAt                                                                                                       <  1:317284/55‑1 (MQ=255)
                                                                                                            cGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTgcgc                  <  1:682293/138‑1 (MQ=255)
                                                                                                            cGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTgcgc                  <  2:313761/138‑1 (MQ=255)
                                                                                                              gTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGca                                                                                     <  2:58436/69‑1 (MQ=255)
                                                                                                              gTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGca                                                                                     >  1:58436/1‑69 (MQ=255)
                                                                                                                  tGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACg                                                      >  1:739270/1‑96 (MQ=255)
                                                                                                                  tGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACg                                                      <  2:739270/96‑1 (MQ=255)
                                                                                                                      ttCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGca                                                                                     <  1:1150112/61‑1 (MQ=255)
                                                                                                                      ttCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGca                                                                                     >  2:1150112/1‑61 (MQ=255)
                                                                                                                         gCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAgc    >  2:493951/1‑139 (MQ=255)
                                                                                                                             ccGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAgcgc  >  2:1091181/1‑137 (MQ=255)
                                                                                                                                       |                                                                                                                              
CGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGC  >  NZ_CP009273/543457‑543718

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ 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
AGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTT  >  NZ_CP009273/543446‑543731
                                                                                                                                                  |                                                                                                                                           
AGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTG                                                                                                                                         >  SRR3722129.61456/1‑151 (MQ=60)
AGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTG                                                                                                                                         >  SRR3722129.527292/1‑151 (MQ=60)
             CGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGT                                                                                                                            >  SRR3722129.1066659/1‑151 (MQ=60)
             CGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGT                                                                                                                            >  SRR3722129.1087286/1‑151 (MQ=60)
             CGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGT                                                                                                                            >  SRR3722129.1001724/1‑151 (MQ=60)
                   GATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTC                                                                                                                      <  SRR3722129.1295658/151‑1 (MQ=60)
                       ATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATctgtct                                                                                                                  >  SRR3722129.133608/1‑145 (MQ=60)
                       ATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCC                                                                                                                  >  SRR3722129.599332/1‑151 (MQ=60)
                       ATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCC                                                                                                                  >  SRR3722129.1170971/1‑151 (MQ=60)
                            CCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCT                                                                                                             >  SRR3722129.458070/1‑151 (MQ=60)
                               CTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGC                                                                                                          >  SRR3722129.1117128/1‑151 (MQ=60)
                               CTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCATGCTGGCGGCAGCGTTAATTCATCCGCTCTCGC                                                                                                          >  SRR3722129.361688/1‑151 (MQ=60)
                                 tttttttcaagcagaagacggcatacgagattcgccttagtctcgtgagctcggagatgtgtataagagacagGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCC                                                                                                         <  SRR3722129.321814/77‑1 (MQ=60)
                                                   tcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACAC                                                                                       <  SRR3722129.1166149/83‑1 (MQ=60)
                                                   GACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACctg                                                                                      >  SRR3722129.408632/1‑148 (MQ=60)
                                                    ACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACAC                                                                                     <  SRR3722129.943706/151‑1 (MQ=60)
                                                     aagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACC                                                                                    <  SRR3722129.922452/86‑1 (MQ=60)
                                                      AGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCG                                                                                   <  SRR3722129.1261654/151‑1 (MQ=60)
                                                                      CCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGAGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACG                                                                   >  SRR3722129.1128662/1‑151 (MQ=60)
                                                                                GCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGAT                                                         >  SRR3722129.226327/1‑151 (MQ=60)
                                                                                    GCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAAC                                                     <  SRR3722129.101289/151‑1 (MQ=60)
                                                                                       GCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACctgtctcttatacacatc                                                  >  SRR3722129.1073215/1‑133 (MQ=60)
                                                                                                          GATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGC                               <  SRR3722129.30400/151‑1 (MQ=60)
                                                                                                          GATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGC                               <  SRR3722129.112239/151‑1 (MQ=60)
                                                                                                              CGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttc                           >  SRR3722129.59205/1‑91 (MQ=60)
                                                                                                               GGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAG                          <  SRR3722129.849786/151‑1 (MQ=60)
                                                                                                                  GATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaagacataaaatcagt                           >  SRR3722129.382805/1‑54 (MQ=60)
                                                                                                                  GATTGCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTctgtctcttatacacatctccgagcccacg                          >  SRR3722129.750089/1‑118 (MQ=60)
                                                                                                                      GCGGTCGTGAGTTCGCTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAA                    <  SRR3722129.692365/150‑1 (MQ=60)
                                                                                                                                     CTACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCACTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGG    >  SRR3722129.1135048/1‑151 (MQ=60)
                                                                                                                                       ACCGGAAGCCATCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTT  >  SRR3722129.1097212/1‑151 (MQ=60)
                                                                                                                                                  |                                                                                                                                           
AGCCTGTGCCGCGCGTTTAGATAATCAACCACTGTATGCCATTGCCACGGTGACAGGCACTGAACGTGACCCGCAATGCCGCTCGCAGCAAATCGCCACGCTGGAAGATGCGGGGATTGCGGTCGTGAGTTCGCTACCGGAAGCCACCTTGCTGGCGGCAGCGTTAATTCATCCGCTCTCGCCTGCCGCACAGCAACACACACCGTCATTACTGGAAAACGTCGCCGTGATTAACATCGGATTACGCAGCTTTGCGCTGGAGCTACAAAGCGCCAGCAAACCGGTT  >  NZ_CP009273/543446‑543731

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