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
A1 F2 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,073,163 A→T D26V (GAT→GTT)  BW25113_RS20185 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,073,1630AT100.0% 332.5 / NA 84D26V (GAT→GTT) BW25113_RS20185hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (49/35);  total (49/35)

GAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCA  >  NZ_CP009273/4073078‑4073248
                                                                                     |                                                                                     
gagaTATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTtctc                                                                                   >  2:465953/1‑90 (MQ=255)
 agaTATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCg                                                                                  >  1:207039/1‑90 (MQ=255)
  gaTATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGa                                                                                 <  1:86007/90‑1 (MQ=255)
  gaTATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGa                                                                                 >  1:448270/1‑90 (MQ=255)
  gaTATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGa                                                                                 >  2:1185549/1‑90 (MQ=255)
  gaTATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGa                                                                                 >  2:443128/1‑90 (MQ=255)
  gaTATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGa                                                                                 <  1:346043/90‑1 (MQ=255)
       ttATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                            <  1:1049606/90‑1 (MQ=255)
       ttATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                            <  2:764946/90‑1 (MQ=255)
       ttATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                            <  1:947034/90‑1 (MQ=255)
             aCTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTa                                                                      >  2:749155/1‑90 (MQ=255)
             aCTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTa                                                                      >  1:201300/1‑90 (MQ=255)
                ctctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATc                                                                   >  1:758297/1‑90 (MQ=255)
                ctctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATc                                                                   >  2:754843/1‑90 (MQ=255)
                ctctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATc                                                                   >  2:1176734/1‑90 (MQ=255)
                  ctGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCtt                                                                 >  2:38848/1‑90 (MQ=255)
                        ggTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGt                                                           >  2:1078398/1‑90 (MQ=255)
                        ggTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGt                                                           >  2:67147/1‑90 (MQ=255)
                          tataGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGc                                                         <  1:38848/90‑1 (MQ=255)
                          tataGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGc                                                         >  2:458639/1‑90 (MQ=255)
                          tataGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGc                                                         >  2:763415/1‑90 (MQ=255)
                          tataGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGc                                                         <  1:754843/90‑1 (MQ=255)
                            taGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCgg                                                       <  1:201173/90‑1 (MQ=255)
                            taGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCgg                                                       <  1:1211037/90‑1 (MQ=255)
                            taGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCgg                                                       <  2:1031244/90‑1 (MQ=255)
                              gATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAc                                                     <  2:448270/90‑1 (MQ=255)
                              gATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAc                                                     <  1:489604/90‑1 (MQ=255)
                               atatGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAcc                                                    <  1:295722/90‑1 (MQ=255)
                               atatGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAcc                                                    <  2:260147/90‑1 (MQ=255)
                               atatGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAcc                                                    <  2:467560/90‑1 (MQ=255)
                                  tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTa                                                                      <  2:944168/69‑1 (MQ=255)
                                  tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTa                                                                      >  1:944168/1‑69 (MQ=255)
                                  tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                                 >  2:1036419/1‑90 (MQ=255)
                                  tGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAt                                                 >  2:577031/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATTGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                            >  2:855514/1‑62 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATTGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                            <  1:855514/62‑1 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGa                                                                                 <  2:755145/57‑1 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGa                                                                                 >  1:755145/1‑57 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                            >  2:391387/1‑62 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGc                                                                            <  1:391387/62‑1 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCGGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  2:1189012/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAcc                                                    <  2:562345/86‑1 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGAcc                                                    >  1:562345/1‑86 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  2:501000/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  2:460440/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  1:192807/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  1:1131512/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  1:211270/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  2:345615/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  2:80245/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  2:148272/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  1:1104216/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  2:1117202/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  1:774945/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  1:1002467/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  1:845775/1‑90 (MQ=255)
                                   gggCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTAtt                                                >  1:566823/1‑90 (MQ=255)
                                        gAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAggg                                           >  1:1092336/1‑90 (MQ=255)
                                        gAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAggg                                           >  1:793098/1‑90 (MQ=255)
                                        gAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAggg                                           >  2:766072/1‑90 (MQ=255)
                                         aaTTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGa                                          <  1:401371/90‑1 (MQ=255)
                                         aaTTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGa                                          <  1:167700/90‑1 (MQ=255)
                                          aTTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGaa                                         <  1:67147/90‑1 (MQ=255)
                                             ttACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCg                                      >  1:1204613/1‑90 (MQ=255)
                                               aCTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGt                                    <  2:1213263/90‑1 (MQ=255)
                                               aCTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGt                                    <  1:577031/90‑1 (MQ=255)
                                               aCTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGt                                    <  1:63894/90‑1 (MQ=255)
                                                cTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTa                                   <  2:567107/90‑1 (MQ=255)
                                                    aTTTTTCGGATGAGGGGATTAAGCAACTTGATGTTTTCGAGGGGTAGTGTAATTTTCCTTGGGAGGAACTATTTAGGGAAGCGGTTGAta                               >  1:979927/1‑89 (MQ=255)
                                                    aTTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATc                               >  1:145079/1‑90 (MQ=255)
                                                    aTTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATc                               <  1:416722/90‑1 (MQ=255)
                                                       tATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAAt                            <  2:1046691/90‑1 (MQ=255)
                                                              tGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGAt                     >  1:526834/1‑90 (MQ=255)
                                                                 ggTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATaaa                  >  1:558603/1‑90 (MQ=255)
                                                                     aTTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAaatcaa              >  1:1230020/1‑90 (MQ=255)
                                                                     aTTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTCCCTCGTGCGGAACTATTAAGGGAAGCGGAAGATCAATACCTGATAaatcaa              >  1:859376/1‑90 (MQ=255)
                                                                         aGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGc          >  1:701649/1‑90 (MQ=255)
                                                                           cAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCaa        <  2:701649/90‑1 (MQ=255)
                                                                           cAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCaa        <  1:859064/90‑1 (MQ=255)
                                                                            aaCTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCaaa       <  1:11234/90‑1 (MQ=255)
                                                                            aaCTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCaaa       <  1:327587/90‑1 (MQ=255)
                                                                            aaCTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCaaa       <  1:794445/90‑1 (MQ=255)
                                                                              cTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACa     >  2:219534/1‑90 (MQ=255)
                                                                               ttGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAg    <  1:1049813/90‑1 (MQ=255)
                                                                                 gatgttCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCa  <  2:566823/90‑1 (MQ=255)
                                                                                 gatgttCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCa  <  2:793098/90‑1 (MQ=255)
                                                                                     |                                                                                     
GAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCA  >  NZ_CP009273/4073078‑4073248

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GATATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTCC  >  NZ_CP009273/4073069‑4073262
                                                                                              |                                                                                                   
GATATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCG                                                                                                >  SRR3722089.209643/1‑100 (MQ=60)
 ATATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGA                                                                                               >  SRR3722089.453651/1‑100 (MQ=60)
           GATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGT                                                                                     <  SRR3722089.87068/100‑1 (MQ=60)
           GATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGT                                                                                     <  SRR3722089.350344/100‑1 (MQ=60)
            ATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTA                                                                                    >  SRR3722089.203838/1‑100 (MQ=60)
               TTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATC                                                                                 >  SRR3722089.768416/1‑100 (MQ=60)
                gatgtgtataagagacgGGTATAGATATGGGCAGAATTTTACTCGATTTATTGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCT                                                                                <  SRR3722089.867720/83‑1 (MQ=60)
                gatgtgtataagagaCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCT                                                                                <  SRR3722089.396157/85‑1 (MQ=60)
                TTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCT                                                                                <  SRR3722089.1065340/100‑1 (MQ=60)
                TTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCT                                                                                <  SRR3722089.960931/100‑1 (MQ=60)
                                 GGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGctgtctcttat                                                               >  SRR3722089.958006/1‑89 (MQ=60)
                                  GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTctgtctcttatacacatctgacg                                                              >  SRR3722089.765197/1‑77 (MQ=60)
                                  GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722089.785453/1‑100 (MQ=60)
                                  GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722089.1120912/1‑100 (MQ=60)
                                  GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722089.857798/1‑100 (MQ=60)
                                  GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722089.568926/1‑100 (MQ=60)
                                  GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722089.573497/1‑100 (MQ=60)
                                  GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722089.1148718/1‑100 (MQ=60)
                                  GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722089.1017371/1‑100 (MQ=60)
                                  GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722089.213937/1‑100 (MQ=60)
                                  GTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATT                                                              >  SRR3722089.195230/1‑100 (MQ=60)
                                   TATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTA                                                             <  SRR3722089.39331/100‑1 (MQ=60)
                                   TATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTA                                                             <  SRR3722089.764892/100‑1 (MQ=60)
                                     TAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAG                                                           <  SRR3722089.1229452/100‑1 (MQ=60)
                                     TAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAG                                                           <  SRR3722089.203709/100‑1 (MQ=60)
                                       GATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGG                                                         >  SRR3722089.1108828/1‑100 (MQ=60)
                                       GATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGG                                                         >  SRR3722089.804024/1‑100 (MQ=60)
                                       GATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGG                                                         <  SRR3722089.495408/100‑1 (MQ=60)
                                        ATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGA                                                        <  SRR3722089.299497/100‑1 (MQ=60)
                                            GGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCG                                                    >  SRR3722089.1222925/1‑100 (MQ=60)
                                                  AATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGAT                                              <  SRR3722089.406245/100‑1 (MQ=60)
                                                  AATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGAT                                              <  SRR3722089.169793/100‑1 (MQ=60)
                                                   ATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATC                                             <  SRR3722089.67961/100‑1 (MQ=60)
                                                   ATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATC                                             >  SRR3722089.146882/1‑100 (MQ=60)
                                                        ACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATAC                                        <  SRR3722089.583887/100‑1 (MQ=60)
                                                        ACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATAC                                        <  SRR3722089.64674/100‑1 (MQ=60)
                                                             ATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGAT                                   >  SRR3722089.532986/1‑100 (MQ=60)
                                                             ATTTATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGAT                                   <  SRR3722089.421762/100‑1 (MQ=60)
                                                                TATCGAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAA                                >  SRR3722089.565101/1‑100 (MQ=60)
                                                                    GAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAA                            >  SRR3722089.1248716/1‑100 (MQ=60)
                                                                    GAATGAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTCCCTCGTGCGGAACTATTAAGGGAAGCGGAAGATCAATACCTGATAAATCAA                            >  SRR3722089.871655/1‑100 (MQ=60)
                                                                        GAGGTGATTAAGCAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGC                        >  SRR3722089.710794/1‑100 (MQ=60)
                                                                                    CAACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAA            <  SRR3722089.871338/100‑1 (MQ=60)
                                                                                     AACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAAC           <  SRR3722089.805391/100‑1 (MQ=60)
                                                                                     AACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAAC           <  SRR3722089.11364/100‑1 (MQ=60)
                                                                                     AACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAAC           <  SRR3722089.331702/100‑1 (MQ=60)
                                                                                      ACTTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACC          >  SRR3722089.972455/1‑100 (MQ=60)
                                                                                        TTGATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAG        <  SRR3722089.1065550/100‑1 (MQ=60)
                                                                                           ATGTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGT     >  SRR3722089.403639/1‑100 (MQ=60)
                                                                                             GTTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTC   <  SRR3722089.535841/100‑1 (MQ=60)
                                                                                              TTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTCC  >  SRR3722089.484410/1‑100 (MQ=60)
                                                                                              TTCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTCC  >  SRR3722089.1178864/1‑100 (MQ=60)
                                                                                              |                                                                                                   
GATATTGGCGAGATATTTATGAACTCTCTGGCAGGTATAGATATGGGCAGAATTTTACTCGATTTATCGAATGAGGTGATTAAGCAACTTGATGATCTCGAGGTGCAGCGTAATCTTCCTCGTGCGGACCTATTAAGGGAAGCGGTAGATCAATACCTGATAAATCAATCGCAAACAGCAAGAACCAGTGTTCC  >  NZ_CP009273/4073069‑4073262

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