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 F3 I152 R1 8 44.1 2758600 95.6% 2637221 84.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,754,231 G→T R127S (CGC→AGC)  selA ← L‑seryl‑tRNA(Sec) selenium transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,754,2310GT100.0% 300.9 / NA 75R127S (CGC→AGC) selAL‑seryl‑tRNA(Sec) selenium transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (28/47);  total (28/47)

CGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  NZ_CP009273/3754145‑3754316
                                                                                      |                                                                                     
cgccgcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTcc                                                                                    <  2:894476/90‑1 (MQ=255)
cgccgcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTcc                                                                                    <  1:736726/90‑1 (MQ=255)
    gcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcg                                                                                <  2:1377422/90‑1 (MQ=255)
     cgttGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:850371/88‑1 (MQ=255)
     cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATCCGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:847997/90‑1 (MQ=255)
     cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  2:372733/90‑1 (MQ=255)
     cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:994088/90‑1 (MQ=255)
          ttATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:472172/90‑1 (MQ=255)
          ttATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:611851/90‑1 (MQ=255)
          ttATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  1:915140/90‑1 (MQ=255)
            aTTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAg                                                                        <  2:325115/90‑1 (MQ=255)
            aTTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAg                                                                        <  1:1357885/90‑1 (MQ=255)
              tGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGAt                                                                      >  2:382287/1‑90 (MQ=255)
                 cGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCAt                                                                   <  1:971313/90‑1 (MQ=255)
                 cGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCAt                                                                   <  2:17223/90‑1 (MQ=255)
                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  1:477840/81‑1 (MQ=255)
                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          >  2:313270/1‑81 (MQ=255)
                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          >  1:509296/1‑81 (MQ=255)
                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          >  2:477840/1‑81 (MQ=255)
                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:509296/81‑1 (MQ=255)
                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  1:313270/81‑1 (MQ=255)
                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  2:594993/1‑90 (MQ=255)
                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  2:1295920/1‑90 (MQ=255)
                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  2:1186700/1‑90 (MQ=255)
                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:722908/1‑90 (MQ=255)
                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:64942/1‑90 (MQ=255)
                   aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:119414/1‑90 (MQ=255)
                             cTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCa                                                       >  1:84349/1‑90 (MQ=255)
                             cTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCa                                                       >  1:702085/1‑90 (MQ=255)
                                ccGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTg                                                    <  1:672279/90‑1 (MQ=255)
                                  gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCa                                                                    <  1:944280/72‑1 (MQ=255)
                                  gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCa                                                                    >  2:944280/1‑72 (MQ=255)
                                  gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGc                                                  <  2:1277048/90‑1 (MQ=255)
                                  gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGc                                                  >  2:1234850/1‑90 (MQ=255)
                                   cccccGAAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCg                                                 <  1:231597/90‑1 (MQ=255)
                                    ccccGTACTACGGCACAGCAACTGCGCCAGCGCCCGTTAGCTATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:381079/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:1265082/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:1254319/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:78030/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:458630/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:1001050/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:559213/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:1234850/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:135701/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:624284/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:64942/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:667339/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:705576/90‑1 (MQ=255)
                                    ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:720751/90‑1 (MQ=255)
                                        gTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGc                                            <  2:449148/90‑1 (MQ=255)
                                              cGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGcc                                      >  2:1174904/1‑90 (MQ=255)
                                                gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATa                                                                  >  2:752147/1‑60 (MQ=255)
                                                gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATa                                                                  <  1:752147/60‑1 (MQ=255)
                                                gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTg                                    >  1:860808/1‑90 (MQ=255)
                                                   cagcagCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               >  1:10088/1‑44 (MQ=255)
                                                   cagcagCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  2:10088/44‑1 (MQ=255)
                                                     gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:912800/73‑1 (MQ=255)
                                                     gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                >  2:912800/1‑73 (MQ=255)
                                                     gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                               <  1:1266956/90‑1 (MQ=255)
                                                        gcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               >  1:968904/1‑39 (MQ=255)
                                                        gcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  2:968904/39‑1 (MQ=38)
                                                        gcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:234336/44‑1 (MQ=255)
                                                        gcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          >  1:234336/1‑44 (MQ=255)
                                                          tGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCt                          >  1:56555/1‑90 (MQ=255)
                                                          tGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCt                          >  2:814993/1‑90 (MQ=255)
                                                             gcCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAg                                                                        <  1:609094/41‑1 (MQ=255)
                                                             gcCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAg                                                                        >  2:609094/1‑41 (MQ=255)
                                                                 gcgcCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGAt                                                                      >  1:601810/1‑39 (MQ=39)
                                                                 gcgcCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGAt                                                                      <  2:601810/39‑1 (MQ=39)
                                                                       gATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCtt             <  2:1357520/90‑1 (MQ=255)
                                                                       gATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCtt             <  1:500480/90‑1 (MQ=255)
                                                                           gcgATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGc         <  2:629537/90‑1 (MQ=255)
                                                                                  ccGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTaaa  >  1:340756/1‑90 (MQ=255)
                                                                                  ccGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTaaa  >  1:1225962/1‑90 (MQ=255)
                                                                                  ccGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTaaa  >  1:922055/1‑90 (MQ=255)
                                                                                      |                                                                                     
CGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  NZ_CP009273/3754145‑3754316

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  NZ_CP009273/3754133‑3754316
                                                                                                  |                                                                                     
CAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTT                                                                                      >  SRR3722022.1073124/1‑100 (MQ=60)
 AATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTC                                                                                     <  SRR3722022.508581/100‑1 (MQ=60)
 AATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCACCAGCGCCCGATCGCGATGTCCGCTTC                                                                                     <  SRR3722022.1021999/100‑1 (MQ=60)
    AACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGG                                                                                  <  SRR3722022.1156998/100‑1 (MQ=60)
            CGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCC                                                                          <  SRR3722022.749496/100‑1 (MQ=60)
                gcactgctcTTGACGATGCAGGCCCCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGAT                                                                      <  SRR3722022.719905/91‑1 (MQ=60)
                 CGTTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATC                                                                     <  SRR3722022.865752/100‑1 (MQ=60)
                 CATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATCCGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATC                                                                     <  SRR3722022.863321/100‑1 (MQ=60)
                 CATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATC                                                                     <  SRR3722022.1012924/100‑1 (MQ=60)
                     GTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAC                                                                 >  SRR3722022.517462/1‑100 (MQ=60)
                     GTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAC                                                                 >  SRR3722022.65988/1‑100 (MQ=60)
                     GTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAC                                                                 >  SRR3722022.735386/1‑100 (MQ=60)
                     GTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAC                                                                 >  SRR3722022.121306/1‑100 (MQ=60)
                      TTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACT                                                                <  SRR3722022.318492/100‑1 (MQ=60)
                      TTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACT                                                                <  SRR3722022.485601/100‑1 (MQ=60)
                      TTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACN                                                                <  SRR3722022.932198/100‑1 (MQ=60)
                        aactcctcgtcggcagcgtcagatgtgtataagagacagCAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCG                                                              <  SRR3722022.618998/61‑1 (MQ=60)
                        ATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCG                                                              <  SRR3722022.1384559/100‑1 (MQ=60)
                            agagacagGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGG                                                          <  SRR3722022.962023/92‑1 (MQ=60)
                             CGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGT                                                         <  SRR3722022.989661/100‑1 (MQ=60)
                              cagatgtgtataagagacagCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTC                                                        <  SRR3722022.765194/80‑1 (MQ=60)
                               ATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCA                                                       >  SRR3722022.714137/1‑100 (MQ=60)
                               ATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCA                                                       >  SRR3722022.85694/1‑100 (MQ=60)
                                            CCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCAT                                          <  SRR3722022.683716/100‑1 (MQ=60)
                                               CCCCCGAAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGC                                       <  SRR3722022.235484/100‑1 (MQ=60)
                                                gagacagATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCC                                      <  SRR3722022.929795/93‑1 (MQ=60)
                                                CCCCGTACTACGGCACAGCAACTGCGCCAGCGCCCGTTAGCTATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATATCC                                      <  SRR3722022.387339/100‑1 (MQ=60)
                                                CCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCC                                      <  SRR3722022.1259120/100‑1 (MQ=60)
                                                CCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCC                                      <  SRR3722022.137869/100‑1 (MQ=60)
                                                  CCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTG                                    >  SRR3722022.876473/1‑100 (MQ=60)
                                                     TAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCctgtctcttatacacatctgacgctgccgacgagga                                 >  SRR3722022.10261/1‑64 (MQ=60)
                                                          CGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCctgtctcttatacacatctgacgctgccgacgaggagttcc                            >  SRR3722022.987211/1‑59 (MQ=60)
                                                          CGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCtgtctcttatacacatctgacgctgccgacgagga                            >  SRR3722022.238274/1‑65 (MQ=60)
                                                            GCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCT                          >  SRR3722022.57480/1‑100 (MQ=60)
                                                                 GCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCAC                     <  SRR3722022.1291898/100‑1 (MQ=60)
                                                                   AGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGctgtctcttatacacatctgacgctgccgacgaggagttcc                   >  SRR3722022.611517/1‑59 (MQ=60)
                                                                                   GATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAA   <  SRR3722022.508530/100‑1 (MQ=60)
                                                                                    ATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  SRR3722022.346412/1‑100 (MQ=60)
                                                                                    ATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  SRR3722022.939270/1‑100 (MQ=60)
                                                                                    ATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  SRR3722022.1250024/1‑100 (MQ=60)
                                                                                                  |                                                                                     
CAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  NZ_CP009273/3754133‑3754316

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