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 929,816 T→C G379G (GGT→GGC ftsK → DNA translocase FtsK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273929,8160TC93.8% 200.6 / 4.2 65G379G (GGT→GGCftsKDNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (29/32);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

CGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTAT  >  NZ_CP009273/929683‑929942
                                                                                                                                     |                                                                                                                              
cgcCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTAcc                                                                                                                            >  2:582953/1‑138 (MQ=255)
    ggTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCAcag                                                                                                                       <  2:676560/139‑1 (MQ=255)
         aaCCTGTGGCTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTAGCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTc                                                                                                                  <  2:560397/139‑1 (MQ=255)
                     aGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCa                                                                                                      <  1:995923/139‑1 (MQ=255)
                      gACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCAcag                                                                                                                       <  2:1268267/121‑1 (MQ=255)
                      gACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCAcag                                                                                                                       >  1:1268267/1‑121 (MQ=255)
                          ccgccTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTg                                                                                                 <  1:1217461/139‑1 (MQ=255)
                              cTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGt                                                                                             <  1:1164090/139‑1 (MQ=255)
                                           tGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCcaca                                                                                                                        <  1:1033991/99‑1 (MQ=255)
                                           tGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCcaca                                                                                                                        >  2:1033991/1‑99 (MQ=255)
                                             aTGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCa                                                                                               <  1:749876/122‑1 (MQ=255)
                                             aTGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCa                                                                                               >  2:749876/1‑122 (MQ=255)
                                              tGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCg                                                                             <  1:1005447/139‑1 (MQ=255)
                                                  ccACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTg                                                                          >  2:746943/1‑138 (MQ=255)
                                                   cACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGca                                                                        <  1:1280398/139‑1 (MQ=255)
                                                    aCCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTcaca                                                                                                               <  1:1098070/99‑1 (MQ=255)
                                                    aCCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTcaca                                                                                                               >  2:1098070/1‑99 (MQ=255)
                                                    aCCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTc                                                                                                                  >  2:757857/1‑96 (MQ=255)
                                                    aCCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTc                                                                                                                  <  1:757857/96‑1 (MQ=255)
                                                    aCCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGcaa                                                                       >  1:440772/1‑139 (MQ=255)
                                                      cTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGcaaca                                                                     <  1:746943/139‑1 (MQ=255)
                                                        gcgcAACCTACAGTAGCGTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCa                                                                                                      <  2:1099575/104‑1 (MQ=255)
                                                        gcgcAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCa                                                                                                      >  1:1099575/1‑104 (MQ=255)
                                                         cgcAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCaacaac                                                                   >  2:145697/1‑138 (MQ=255)
                                                           cAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAg                                                                >  1:561819/1‑139 (MQ=255)
                                                              ccTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTa                                                              >  2:1069745/1‑138 (MQ=255)
                                                                       gCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCgcag                                                     <  1:555477/138‑1 (MQ=255)
                                                                        ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCcaca                                                                                                                        <  2:200579/70‑1 (MQ=255)
                                                                        ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCcaca                                                                                                                        >  1:200579/1‑70 (MQ=255)
                                                                        ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCAcag                                                                                                                       <  1:708055/71‑1 (MQ=255)
                                                                        ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCAcag                                                                                                                       >  2:708055/1‑71 (MQ=255)
                                                                        ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTAc                                                             >  2:831248/1‑129 (MQ=255)
                                                                        ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTAc                                                             <  1:831248/129‑1 (MQ=255)
                                                                        ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCgcag                                                     >  2:555477/1‑137 (MQ=255)
                                                                        ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCgcag                                                     >  2:799854/1‑137 (MQ=255)
                                                                        ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCgcag                                                     >  2:987821/1‑137 (MQ=255)
                                                                        ccTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCgcag                                                     >  1:122260/1‑137 (MQ=255)
                                                                         cTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGcagcag                                                  <  1:1049364/139‑1 (MQ=255)
                                                                         cTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGcagcag                                                  <  1:736951/139‑1 (MQ=255)
                                                                          tGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCagcagc                                                 >  1:1173384/1‑139 (MQ=255)
                                                                          tGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCagcagc                                                 >  1:261933/1‑139 (MQ=255)
                                                                            gCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCagcagc                                                 <  1:582953/137‑1 (MQ=255)
                                                                            gCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCg                                               <  2:921134/139‑1 (MQ=255)
                                                                            gCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCg                                               >  1:650834/1‑139 (MQ=255)
                                                                                         ggTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTg                                  <  2:613978/139‑1 (MQ=255)
                                                                                         ggTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGAAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTg                                  <  2:425610/139‑1 (MQ=255)
                                                                                           tCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCa                                <  1:799854/139‑1 (MQ=255)
                                                                                            ccACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATg                                                                                  <  2:620411/88‑1 (MQ=255)
                                                                                            ccACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATg                                                                                  >  1:620411/1‑88 (MQ=255)
                                                                                              acaAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCt                             >  2:459123/1‑139 (MQ=255)
                                                                                              acaAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCt                             >  1:120746/1‑139 (MQ=255)
                                                                                                     ggAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGaacaac                      >  1:1242749/1‑139 (MQ=255)
                                                                                                            cGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGACGTATTATGCACCTGCAGCTGAACAACCTgcgc                >  1:977992/1‑138 (MQ=255)
                                                                                                             ggTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCaa              <  2:654829/139‑1 (MQ=255)
                                                                                                              gTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTg                                                                                            <  1:4467/60‑1 (MQ=255)
                                                                                                              gTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTg                                                                                            >  2:4467/1‑60 (MQ=255)
                                                                                                                  ttGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGc          <  2:650834/138‑1 (MQ=255)
                                                                                                                        cTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGtatta   <  2:122260/139‑1 (MQ=255)
                                                                                                                           cACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGtattat  <  1:145697/137‑1 (MQ=255)
                                                                                                                             ccGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTg                                                                                                 <  2:757174/40‑1 (MQ=255)
                                                                                                                             ccGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTg                                                                                                 >  1:757174/1‑40 (MQ=255)
                                                                                                                                gAAGGTTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCg                                                                             <  1:890602/57‑1 (MQ=255)
                                                                                                                                gAAGGTTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCg                                                                             >  2:890602/1‑57 (MQ=255)
                                                                                                                                gAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGc          <  2:635177/124‑1 (MQ=255)
                                                                                                                                gAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGc          >  1:635177/1‑124 (MQ=255)
                                                                                                                                     |                                                                                                                              
CGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTAT  >  NZ_CP009273/929683‑929942

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAAC  >  NZ_CP009273/929671‑929961
                                                                                                                                                 |                                                                                                                                                 
AAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCC                                                                                                                                              <  SRR3722129.681470/151‑1 (MQ=60)
  AGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCAC                                                                                                                                            <  SRR3722129.934070/151‑1 (MQ=60)
              cgtgggctcggagatgtgtataagagacagGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAAT                                                                                                                                <  SRR3722129.1048551/121‑1 (MQ=60)
               agacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATA                                                                                                                               <  SRR3722129.718465/93‑1 (MQ=60)
                       cgtgggctcggagatgtgtataagagacagGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGC                                                                                                                       <  SRR3722129.1113398/121‑1 (MQ=60)
                       cgtgggctcggagatgtgtataagagacagGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGC                                                                                                                          <  SRR3722129.768901/118‑1 (MQ=60)
                       CCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATActgtctct                                                                                                                       >  SRR3722129.1285885/1‑143 (MQ=60)
                                CAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGC                                                                                                              <  SRR3722129.1009970/151‑1 (MQ=60)
                                     GCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATAT                                                                                                         <  SRR3722129.1234416/151‑1 (MQ=60)
                                       gagacagTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAA                                                                                                       <  SRR3722129.760823/144‑1 (MQ=60)
                                         CCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATG                                                                                                     <  SRR3722129.1180323/151‑1 (MQ=60)
                                           ttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagCGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGA                                                                                                    <  SRR3722129.4529/82‑1 (MQ=60)
                                                     GTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAAC                                                                                         >  SRR3722129.447167/1‑151 (MQ=60)
                                                         ATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCctgtctcttatacacatctccgagc                                                                                     >  SRR3722129.1114912/1‑126 (MQ=60)
                                                         ATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACC                                                                                     <  SRR3722129.1019610/151‑1 (MQ=60)
                                                          tttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacaGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACC                                                                                     <  SRR3722129.903262/80‑1 (MQ=60)
                                                            TTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGT                                                                                  >  SRR3722129.570157/1‑151 (MQ=60)
                                                              CCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTAC                                                                                <  SRR3722129.1298168/151‑1 (MQ=60)
                                                                 CCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGC                                                                             <  SRR3722129.757854/151‑1 (MQ=60)
                                                                         ACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtctt                                                                     >  SRR3722129.203384/1‑92 (MQ=60)
                                                                         ACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAG                                                                     <  SRR3722129.843196/151‑1 (MQ=60)
                                                                         ACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGC                                                                       >  SRR3722129.123929/1‑149 (MQ=60)
                                                                           CTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCC                                                                   >  SRR3722129.1189719/1‑151 (MQ=60)
                                                                           CTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCC                                                                   >  SRR3722129.265630/1‑151 (MQ=60)
                                                                             ACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGT                                                                 >  SRR3722129.660532/1‑151 (MQ=60)
                                                                                  AGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTA                                                             <  SRR3722129.563738/150‑1 (MQ=60)
                                                                                    CCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGC                                                          <  SRR3722129.747742/151‑1 (MQ=60)
                                                                                    CCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGC                                                          <  SRR3722129.1064096/151‑1 (MQ=60)
                                                                                       GGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCA                                                         <  SRR3722129.591562/149‑1 (MQ=60)
                                                                                             CTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACtgtctcttatacacatctccgagcccacgagactaaggcg                                                 >  SRR3722129.629632/1‑111 (MQ=60)
                                                                                               GTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTG                                               >  SRR3722129.122401/1‑151 (MQ=60)
                                                                                                      GTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACC                                        >  SRR3722129.1259997/1‑151 (MQ=60)
                                                                                                      GTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACC                                        <  SRR3722129.811410/151‑1 (MQ=60)
                                                                                                             AACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGACGTATTATGCACCTGCAGCTGAACAACCTGCGCA                                  >  SRR3722129.991830/1‑150 (MQ=60)
                                                                                                                              TTGCTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaaaaaaa                 >  SRR3722129.768212/1‑62 (MQ=60)
                                                                                                                                 CTCCTGCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCtgtc             >  SRR3722129.644621/1‑147 (MQ=60)
                                                                                                                                      GCACCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCC          <  SRR3722129.147731/149‑1 (MQ=60)
                                                                                                                                         CCGGAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaacacaaagtcgagggtgacaggg     >  SRR3722129.239577/1‑51 (MQ=60)
                                                                                                                                            GAAGGCTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAAC  >  SRR3722129.230681/1‑151 (MQ=60)
                                                                                                                                                 |                                                                                                                                                 
AAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTGCTCCTGCACCGGAAGGTTACCCACAGCAGTCACAATATGCGCAGCCTGCAGTGCAATATAATGAGCCGCTGCAACAACCAGTACAGCCGCAGCAGCCGTATTATGCACCTGCAGCTGAACAACCTGCGCAACAGCCGTATTATGCCCCTGCGCCAGAACAAC  >  NZ_CP009273/929671‑929961

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