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 599,729 C→T A131T (GCA→ACA)  ybdG ← mechanosensitive ion channel YbdG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273599,7290CT89.3% 176.1 / 9.1 56A131T (GCA→ACA) ybdGmechanosensitive ion channel YbdG
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base T (20/30);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

ACCAAGACCGCTGATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGCCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATAC  >  NZ_CP009273/599601‑599864
                                                                                                                                |                                                                                                                                       
aCCAAGACCGCTGATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTc                                                                                                                               >  1:1089285/1‑139 (MQ=255)
   aaGACCGCTGATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTgc                                                                                                                            >  2:303146/1‑139 (MQ=255)
          cTGATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGAt                                                                                                                     >  1:832401/1‑139 (MQ=255)
          cTGATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGAt                                                                                                                     >  1:266212/1‑139 (MQ=255)
                    tCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAAc                                                                                                           <  1:541543/139‑1 (MQ=255)
                       cTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATc                                                                                                        >  1:1195175/1‑139 (MQ=255)
                       cTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGCCGGGAATTTCTGCGCCAGATTCAAAATAACATc                                                                                                        >  1:397171/1‑139 (MQ=255)
                           cGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGc                                                                                                    <  1:507770/139‑1 (MQ=255)
                                  ccAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACgaga                                                                                             <  1:597235/139‑1 (MQ=255)
                                  ccAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGCCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACgaga                                                                                             <  2:397171/139‑1 (MQ=255)
                                    aaTCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGagaag                                                                                           <  1:275956/139‑1 (MQ=255)
                                           agcGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGa                                                                                      >  2:810623/1‑137 (MQ=255)
                                            gcGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAag                                                                                   <  2:941925/139‑1 (MQ=255)
                                            gcGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAag                                                                                   <  2:459649/139‑1 (MQ=255)
                                            gcGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAag                                                                                   <  2:249388/139‑1 (MQ=255)
                                             cGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAag                                                                                   >  1:408605/1‑138 (MQ=255)
                                               agaTCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAagcag                                                                                <  1:339809/139‑1 (MQ=255)
                                                gaTCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAagcagc                                                                               <  1:360092/139‑1 (MQ=255)
                                                    atcaAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCa                                                                            <  2:1009905/138‑1 (MQ=255)
                                                    atcaAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCa                                                                            <  2:408605/138‑1 (MQ=255)
                                                    atcaAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCa                                                                            <  2:266212/138‑1 (MQ=255)
                                                    atcaAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGCCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCa                                                                            <  2:294256/138‑1 (MQ=255)
                                                          aTGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATAcatca                                                                     <  1:665187/139‑1 (MQ=255)
                                                          aTGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATAcatca                                                                     <  1:703605/139‑1 (MQ=255)
                                                          aTGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATAcatca                                                                     <  1:980905/139‑1 (MQ=255)
                                                          aTGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTGTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATAcatca                                                                     <  1:547327/139‑1 (MQ=255)
                                                             ccGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATg                                                                   <  2:832401/138‑1 (MQ=255)
                                                                         gcgcCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTgcg                                                      <  2:31164/139‑1 (MQ=255)
                                                                          cgcCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTgcgc                                                     >  1:1000715/1‑139 (MQ=255)
                                                                          cgcCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTgcgc                                                     >  1:44285/1‑139 (MQ=255)
                                                                              gATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAg                                                 <  1:159151/139‑1 (MQ=255)
                                                                               aTCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAgg                                                <  2:475881/139‑1 (MQ=255)
                                                                                         tCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGaa                                       >  2:20354/1‑138 (MQ=255)
                                                                                              tGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGcc                                                                                                                         <  1:1163325/51‑1 (MQ=255)
                                                                                              tGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGcc                                                                                                                         >  2:1163325/1‑51 (MQ=255)
                                                                                              tGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACgaga                                                                                             <  1:835702/79‑1 (MQ=255)
                                                                                              tGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACgaga                                                                                             >  2:835702/1‑79 (MQ=255)
                                                                                              tGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCa                                                                            >  1:1004318/1‑96 (MQ=255)
                                                                                              tGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCa                                                                            <  2:1004318/96‑1 (MQ=255)
                                                                                              tGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTg                                 >  2:952610/1‑139 (MQ=255)
                                                                                              tGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTg                                 >  1:214578/1‑139 (MQ=255)
                                                                                              tGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTg                                 >  2:443081/1‑139 (MQ=255)
                                                                                                      cccTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCg                         >  2:811904/1‑139 (MQ=255)
                                                                                                      cccTTTCAGCGGTAACTGAGATGCTGCCGGGAATTTCTGCGCCAGATTCaaaa                                                                                                               <  1:60307/53‑1 (MQ=255)
                                                                                                      cccTTTCAGCGGTAACTGAGATGCTGCCGGGAATTTCTGCGCCAGATTCaaaa                                                                                                               >  2:60307/1‑53 (MQ=255)
                                                                                                          ttCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATAcatca                                                                     >  1:684320/1‑91 (MQ=255)
                                                                                                          ttCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATAcatca                                                                     <  2:684320/91‑1 (MQ=255)
                                                                                                          ttCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGc                     >  1:231116/1‑139 (MQ=255)
                                                                                                           tCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGcc                                                                                                                         <  2:159186/38‑1 (MQ=255)
                                                                                                           tCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGcc                                                                                                                         >  1:159186/1‑38 (MQ=255)
                                                                                                               cGGTAACTGAGATGCTGCCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCtttt                >  2:669418/1‑139 (MQ=255)
                                                                                                                  tAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGc             <  2:708754/139‑1 (MQ=255)
                                                                                                                  tAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGc             <  1:811904/139‑1 (MQ=255)
                                                                                                                  tAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGc             <  1:230394/139‑1 (MQ=255)
                                                                                                                  tAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGc             <  2:1000715/139‑1 (MQ=255)
                                                                                                                     cTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGc          <  1:439419/139‑1 (MQ=255)
                                                                                                                           tgctgtCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATAcatca                                                                     >  2:669964/1‑74 (MQ=255)
                                                                                                                           tgctgtCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATAcatca                                                                     <  1:669964/74‑1 (MQ=255)
                                                                                                                             ctgtCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATAc  <  1:810623/139‑1 (MQ=255)
                                                                                                                             ctgcCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATAc  <  1:669418/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                                       
ACCAAGACCGCTGATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGCCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATAC  >  NZ_CP009273/599601‑599864

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ 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
GCAGCCATTGCACCAAGACCGCTGATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGCCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATACCGCCTGAATATTGA  >  NZ_CP009273/599590‑599878
                                                                                                                                           |                                                                                                                                                     
GCAGCCATTGCACCAAGACCGCTGATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCT                                                                                                                                            >  SRR3722129.1104527/1‑151 (MQ=60)
          CACCAAGACCGCTGATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATT                                                                                                                                  >  SRR3722129.844361/1‑151 (MQ=60)
          CACCAAGACCGCTGATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATT                                                                                                                                  >  SRR3722129.269963/1‑151 (MQ=60)
                 atattttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAAT                                                                                                                            <  SRR3722129.1179549/73‑1 (MQ=60)
                       GATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCC                                                                                                                     >  SRR3722129.1211803/1‑151 (MQ=60)
                       GATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGCCGGGAATTTCTGCGCCAGATTCAAAATAACATCC                                                                                                                     >  SRR3722129.402952/1‑151 (MQ=60)
                           tttttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGCCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGC                                                                                                                  <  SRR3722129.61110/75‑1 (MQ=60)
                              ATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACG                                                                                                              <  SRR3722129.549470/151‑1 (MQ=60)
                                     GCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGAC                                                                                                       <  SRR3722129.515208/151‑1 (MQ=60)
                                            ACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGC                                                                                                <  SRR3722129.606151/151‑1 (MQ=60)
                                             acgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGC                                                                                                <  SRR3722129.847701/101‑1 (MQ=60)
                                             CCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGC                                                                                                >  SRR3722129.414529/1‑150 (MQ=60)
                                              CAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAG                                                                                              <  SRR3722129.279921/151‑1 (MQ=60)
                                                       tttttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACA                                                                                      <  SRR3722129.1156589/75‑1 (MQ=60)
                                                         GAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCA                                                                                   <  SRR3722129.344715/151‑1 (MQ=60)
                                                          AGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCAT                                                                                  <  SRR3722129.365286/151‑1 (MQ=60)
                                                                    cggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAAC                                                                        <  SRR3722129.679870/96‑1 (MQ=60)
                                                                    AATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAAC                                                                        <  SRR3722129.675039/151‑1 (MQ=60)
                                                                    AATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAAC                                                                        <  SRR3722129.713965/151‑1 (MQ=60)
                                                                    AATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAAC                                                                        <  SRR3722129.994779/151‑1 (MQ=60)
                                                                    AATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTGTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAAC                                                                        <  SRR3722129.555378/151‑1 (MQ=60)
                                                                          GACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCG                                                                  >  SRR3722129.1014821/1‑151 (MQ=60)
                                                                          GACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCG                                                                  >  SRR3722129.44839/1‑151 (MQ=60)
                                                                                        CGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAAT                                                    <  SRR3722129.161364/151‑1 (MQ=60)
                                                                                              GTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCtgtctcttatacacatctccgagcccacg                                                 >  SRR3722129.1018457/1‑119 (MQ=60)
                                                                                              GTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGC                                              >  SRR3722129.217571/1‑151 (MQ=60)
                                                                                                          GAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACctgtctcttatacacatctccgagcccacgagacta                                    >  SRR3722129.694416/1‑113 (MQ=60)
                                                                                                          GAAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCC                                  >  SRR3722129.234405/1‑151 (MQ=60)
                                                                                                           AAATATCCCTTTCAGCGGTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaaaaagact                                  >  SRR3722129.161399/1‑60 (MQ=60)
                                                                                                                            GTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATAC                <  SRR3722129.233668/151‑1 (MQ=60)
                                                                                                                            GTAACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATAC                <  SRR3722129.823579/151‑1 (MQ=60)
                                                                                                                               ACTGAGATGCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATACCGC             <  SRR3722129.445794/151‑1 (MQ=60)
                                                                                                                                       GCTGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATACCGCCTGAATAT     <  SRR3722129.822287/151‑1 (MQ=60)
                                                                                                                                       GCTGCCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATACCGCCTGAATAT     <  SRR3722129.679322/151‑1 (MQ=60)
                                                                                                                                         TGTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATACCGCCTGAATATTG   >  SRR3722129.409076/1‑151 (MQ=60)
                                                                                                                                          GTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTctgtctcttatacacatctccgagcccacgagactaaggcgaatct  >  SRR3722129.220319/1‑105 (MQ=60)
                                                                                                                                          GTCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATACCGCCTGAATATTG   >  SRR3722129.759953/1‑150 (MQ=60)
                                                                                                                                           |                                                                                                                                                     
GCAGCCATTGCACCAAGACCGCTGATCAGAATCGCTGGCGACTGACCAATCAGCAGCGAGATCATCAAAATGCCGACCAGAATCGCGCCGATCAGTTTAATCCCCTGAAATATCCCTTTCAGCGGTAACTGAGATGCTGCCGGGAATTTCTGCGCCAGATTCAAAATAACATCCAGCAACGAGAAGACTGAAAGCAGCGCATACATCATGATCCACAACTGCGCGCAGGTAGTCAGAATATCTGCCGCTTCGGTGCCTTTTTGCAGCCAGAATACCGCCTGAATATTGA  >  NZ_CP009273/599590‑599878

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