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 721,666 A→G V276A (GTC→GCC)  kdpB ← potassium‑transporting ATPase subunit KdpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273721,6660AG96.4% 183.8 / ‑0.9 56V276A (GTC→GCC) kdpBpotassium‑transporting ATPase subunit KdpB
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (23/31);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CGACGCCTGACGGTTACCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCG  >  NZ_CP009273/721533‑721780
                                                                                                                                     |                                                                                                                  
cGACGCCTGACGGTTACCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGcc                                                                                                                >  1:748204/1‑138 (MQ=255)
      cTGACGGTTACCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCt                                                                                                         <  1:96412/139‑1 (MQ=255)
       tGACGGTTACCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCt                                                                                                         >  1:663733/1‑138 (MQ=255)
               aCCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAgg                                                                                                >  1:584988/1‑139 (MQ=255)
                     tgtgATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGTCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAg                                                                                                 <  2:436904/132‑1 (MQ=255)
                     tgtgATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGTCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAg                                                                                                 >  1:436904/1‑132 (MQ=255)
                         aTGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATa                                                                                       <  1:548274/138‑1 (MQ=255)
                         aTGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAg                                                                                      <  1:132278/139‑1 (MQ=255)
                             tGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGt                                                                                  <  1:724898/139‑1 (MQ=255)
                             tGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGt                                                                                  >  1:710094/1‑139 (MQ=255)
                             tGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGt                                                                                  >  2:727200/1‑139 (MQ=255)
                             tGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGt                                                                                  >  2:254477/1‑139 (MQ=255)
                             tGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGt                                                                                  <  1:1156321/139‑1 (MQ=255)
                               ccGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAg                                                                                                 >  1:382860/1‑122 (MQ=255)
                               ccGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAg                                                                                                 <  2:382860/122‑1 (MQ=255)
                               ccGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGgcc                                                                                              <  2:1280031/125‑1 (MQ=255)
                               ccGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGgcc                                                                                              >  1:1280031/1‑125 (MQ=255)
                               ccGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGt                                                                                  >  1:72958/1‑137 (MQ=255)
                               ccGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGt                                                                                  <  2:72958/137‑1 (MQ=255)
                               ccGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTg                                                                                >  2:1134192/1‑139 (MQ=255)
                               ccGGATCTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATAACATTCCCGCCTATCCTCCGCCTCATCCCGGCCGCACCGATCGCAGCCAACCGTCCGCCATGAGTGGGTg                                                                                >  2:993106/1‑139 (MQ=255)
                                 ggTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTggg                                                                              <  2:1057503/139‑1 (MQ=255)
                                 ggTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTgg                                                                               >  1:689204/1‑138 (MQ=255)
                                        tCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATcaga                                                                        <  1:727200/138‑1 (MQ=255)
                                        tCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATcaga                                                                        <  1:955143/138‑1 (MQ=255)
                                        tCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATcaga                                                                        <  2:584988/138‑1 (MQ=255)
                                             tAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAg                                                                                                 >  1:1247631/1‑108 (MQ=255)
                                             tAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAg                                                                                                 <  2:1247631/108‑1 (MQ=255)
                                               gCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACca                                                                 <  1:290267/138‑1 (MQ=255)
                                                  gAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATcaga                                                                        >  2:359423/1‑128 (MQ=255)
                                                  gAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATcaga                                                                        <  1:359423/128‑1 (MQ=255)
                                                   aacgtcaacgtcACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCagcagc                                                            <  2:873953/139‑1 (MQ=255)
                                                   aacgtcaacgtcACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCagcagc                                                            <  2:352717/139‑1 (MQ=255)
                                                           cgtcaCCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAg                                                    <  2:583777/139‑1 (MQ=255)
                                                                       cTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCt                                        <  2:710094/139‑1 (MQ=255)
                                                                       cTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCt                                        <  2:748204/139‑1 (MQ=255)
                                                                       cTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCt                                        <  2:948869/139‑1 (MQ=255)
                                                                       cTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCt                                        <  1:1044488/139‑1 (MQ=255)
                                                                       cTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCt                                        <  1:861301/139‑1 (MQ=255)
                                                                          cAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGAc                                     >  1:693677/1‑139 (MQ=255)
                                                                            aCTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTg                                   <  1:758522/139‑1 (MQ=255)
                                                                                cACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCAt                                >  2:313226/1‑138 (MQ=255)
                                                                                     ccGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGc                          >  1:1275407/1‑139 (MQ=255)
                                                                                       gCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGccc                        <  1:79016/139‑1 (MQ=255)
                                                                                        cTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGccc                        >  2:77633/1‑138 (MQ=255)
                                                                                           gtggCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCAcg                    <  2:1238209/139‑1 (MQ=255)
                                                                                                 aTCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGaa               >  2:885196/1‑138 (MQ=255)
                                                                                                      aTTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc         <  1:1162294/139‑1 (MQ=255)
                                                                                                         cgcgCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCca       <  1:313226/138‑1 (MQ=255)
                                                                                                           cgcCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCca       >  2:590981/1‑136 (MQ=255)
                                                                                                           cgcCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCca       <  1:590981/136‑1 (MQ=255)
                                                                                                             ccTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCg  <  2:1275407/139‑1 (MQ=255)
                                                                                                              cTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCca       >  1:885027/1‑133 (MQ=255)
                                                                                                              cTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCca       <  2:885027/133‑1 (MQ=255)
                                                                                                                        gCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCt                                        >  1:61075/1‑90 (MQ=255)
                                                                                                                        gCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCt                                        <  2:61075/90‑1 (MQ=255)
                                                                                                                                     |                                                                                                                  
CGACGCCTGACGGTTACCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCG  >  NZ_CP009273/721533‑721780

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
GGGATAAACTCCGACGCCTGACGGTTACCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGC  >  NZ_CP009273/721522‑721814
                                                                                                                                                |                                                                                                                                                    
GGGATAAACTCCGACGCCTGACGGTTACCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCG                                                                                                                                                 >  SRR3722129.759129/1‑150 (MQ=60)
       ACTCCGACGCCTGACGGTTACCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTG                                                                                                                                          >  SRR3722129.673574/1‑150 (MQ=60)
               GCCTGACGGTTACCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGC                                                                                                                                 >  SRR3722129.593632/1‑151 (MQ=60)
                CCTGACGGTTACCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCC                                                                                                                                <  SRR3722129.97718/151‑1 (MQ=60)
                     CGGTTACCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGTCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCA                                                                                                                            >  SRR3722129.443243/1‑150 (MQ=60)
                             GAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTT                                                                                                                   >  SRR3722129.720530/1‑151 (MQ=60)
                               GTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGctgtctc                                                                                                                 >  SRR3722129.388452/1‑144 (MQ=60)
                               GTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGctgt                                                                                                                 >  SRR3722129.1297795/1‑147 (MQ=60)
                               GTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGG                                                                                                                 >  SRR3722129.73908/1‑151 (MQ=60)
                                 GTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGG                                                                                                                >  SRR3722129.699374/1‑150 (MQ=60)
                                   GATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATC                                                                                                             <  SRR3722129.134064/151‑1 (MQ=60)
                                   GATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGAT                                                                                                              <  SRR3722129.556343/150‑1 (MQ=60)
                                       GTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGAC                                                                                                         <  SRR3722129.735519/151‑1 (MQ=60)
                                       GTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGAC                                                                                                         <  SRR3722129.1172455/151‑1 (MQ=60)
                                             GTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGctgtctcttatacacatctcc                                                                                                   >  SRR3722129.1265015/1‑130 (MQ=60)
                                                 GATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAG                                                                                               <  SRR3722129.364604/151‑1 (MQ=60)
                                                  ATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAG                                                                                               <  SRR3722129.968638/150‑1 (MQ=60)
                                                  ATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAG                                                                                               <  SRR3722129.737847/150‑1 (MQ=60)
                                                         AGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACC                                                                                        <  SRR3722129.294432/150‑1 (MQ=60)
                                                                          ACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACT                                                                      >  SRR3722129.703921/1‑151 (MQ=60)
                                                                                 GCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTAC                                                               <  SRR3722129.1059158/151‑1 (MQ=60)
                                                                                 GCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTAC                                                               <  SRR3722129.873575/151‑1 (MQ=60)
                                                                                     CAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCC                                                           >  SRR3722129.1293118/1‑151 (MQ=60)
                                                                                      AACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCC                                                          <  SRR3722129.769570/151‑1 (MQ=60)
                                                                                                 CGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAA                                               <  SRR3722129.80047/151‑1 (MQ=60)
                                                                                                              CACATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTG                                  >  SRR3722129.897632/1‑151 (MQ=60)
                                                                                                                CATTCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGC                                <  SRR3722129.1178505/151‑1 (MQ=60)
                                                                                                                   TCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG                              <  SRR3722129.599718/150‑1 (MQ=60)
                                                                                                                   TCGCGCCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGG                              <  SRR3722129.317700/150‑1 (MQ=60)
                                                                                                                        CCTAGCATCCGGCTCATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCtgtctcttatacacatctccgagcccacgagactaagg                        >  SRR3722129.61884/1‑113 (MQ=60)
                                                                                                                                      CATCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGAT           >  SRR3722129.980575/1‑150 (MQ=60)
                                                                                                                                        TCCCGGCGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAACGGCCACAGCGTGGCGGTTGGCAGGAAAAAGAagaaaggc       >  SRR3722129.285534/1‑143 (MQ=60)
                                                                                                                                              CGGCGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGC  <  SRR3722129.552670/151‑1 (MQ=60)
                                                                                                                                                |                                                                                                                                                    
GGGATAAACTCCGACGCCTGACGGTTACCGAGTGTGATGGTGCCGGTTTTATCCAGTAGCAGAACGTCAACGTCACCTGCCGCTTCAACTGCACGTCCGCTGGTGGCAATCACATTCGCGCCTAGCATCCGGCTCATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGCGGTTGCCAGTAAAAAGACGATAGTCAGGGC  >  NZ_CP009273/721522‑721814

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