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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,849,816 T→C E117G (GAA→GGA)  ygbA ← nitrous oxide‑stimulated promoter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,849,8160TC96.8% 194.8 / ‑1.1 63E117G (GAA→GGA) ygbAnitrous oxide‑stimulated promoter family protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (30/31);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTT  >  NZ_CP009273/2849730‑2849893
                                                                                      |                                                                             
gTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCAt                                                                            <  2:497437/90‑1 (MQ=255)
  cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                          >  2:794018/1‑90 (MQ=255)
  cACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGa                                                                          >  1:1210493/1‑90 (MQ=255)
       aTGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGttt                                                                     >  1:291698/1‑90 (MQ=255)
          aaGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGtttttt                                                                  >  2:702748/1‑90 (MQ=255)
          aaGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGtttttt                                                                  >  1:778434/1‑90 (MQ=255)
          aaGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGtttttt                                                                  >  2:1475104/1‑90 (MQ=255)
          aaGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGtttttt                                                                  >  1:1422583/1‑90 (MQ=255)
           aGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTc                                                                 <  2:205369/90‑1 (MQ=255)
              ggCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGa                                                              >  2:542323/1‑90 (MQ=255)
              ggCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGa                                                              >  2:324815/1‑90 (MQ=255)
              ggCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGa                                                              <  1:917292/90‑1 (MQ=255)
                   ggAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGata                                                         <  1:1221436/90‑1 (MQ=255)
                         gATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTc                                                   >  2:714473/1‑90 (MQ=255)
                         gATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTc                                                   >  1:1168854/1‑90 (MQ=255)
                              aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                              <  1:1390681/90‑1 (MQ=255)
                              aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                              <  1:1364311/90‑1 (MQ=255)
                              aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                              >  2:829156/1‑90 (MQ=255)
                              aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                              <  1:1012615/90‑1 (MQ=255)
                              aTCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCa                                              <  2:970163/90‑1 (MQ=255)
                                      tGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGc                                      <  2:1104519/90‑1 (MQ=255)
                                        aTTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCAc                                    <  2:657999/90‑1 (MQ=255)
                                           tAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCgg                                 >  1:847592/1‑90 (MQ=255)
                                            aaCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                >  1:421588/1‑90 (MQ=255)
                                            aaCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                >  2:492500/1‑90 (MQ=255)
                                            aaCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                >  1:54148/1‑90 (MQ=255)
                                            aaCTGCAAATTGCCAGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGa                                <  1:1512811/90‑1 (MQ=255)
                                             aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGAc                               <  1:829927/90‑1 (MQ=255)
                                             aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGAc                               <  2:778434/90‑1 (MQ=255)
                                             aCTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGAc                               <  2:181819/90‑1 (MQ=255)
                                                gCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCt                            >  1:819992/1‑90 (MQ=255)
                                                gCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCAATGAGGTTTTTTCGGACGATATGTTTCCGGCAGTTCTGGCACCGGACGCt                            >  2:1190817/1‑90 (MQ=255)
                                                 cAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCtt                           <  1:790282/90‑1 (MQ=255)
                                                 cAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCtt                           <  2:1210493/90‑1 (MQ=255)
                                                       gCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATc                     >  1:1375415/1‑90 (MQ=255)
                                                            aCAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatga                <  1:542323/90‑1 (MQ=255)
                                                            aCAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatga                <  1:324815/90‑1 (MQ=255)
                                                               gATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagat             >  2:648483/1‑90 (MQ=255)
                                                               gATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatgagat             >  1:677922/1‑90 (MQ=255)
                                                                    gCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc        <  1:702748/90‑1 (MQ=255)
                                                                    gCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc        >  1:745703/1‑90 (MQ=255)
                                                                    gCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGc        <  1:601318/90‑1 (MQ=255)
                                                                     ccTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa       <  2:1449673/90‑1 (MQ=255)
                                                                     ccTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa       <  1:273211/90‑1 (MQ=255)
                                                                     ccTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCa       <  2:452736/90‑1 (MQ=255)
                                                                      cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      <  2:568372/90‑1 (MQ=255)
                                                                      cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      <  1:1119859/90‑1 (MQ=255)
                                                                      cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      <  2:847592/90‑1 (MQ=255)
                                                                      cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      <  1:894295/90‑1 (MQ=255)
                                                                      cTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      <  1:714473/90‑1 (MQ=255)
                                                                        gTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACgg    >  1:574768/1‑90 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatg                 <  2:1163428/76‑1 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCgatg                 >  1:1163428/1‑76 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt   >  2:209778/1‑90 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt   >  2:1472290/1‑90 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt   >  2:1281995/1‑90 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt   >  2:1276391/1‑90 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt   >  1:143117/1‑90 (MQ=255)
                                                                         tCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGt   >  1:1189130/1‑90 (MQ=255)
                                                                          ccGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGtt  <  2:1442250/90‑1 (MQ=255)
                                                                          ccGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGtt  <  1:794018/90‑1 (MQ=255)
                                                                                  cTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      <  2:466306/78‑1 (MQ=255)
                                                                                  cTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAc      >  1:466306/1‑78 (MQ=255)
                                                                                      |                                                                             
GTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTT  >  NZ_CP009273/2849730‑2849893

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ 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
GCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCA  >  NZ_CP009273/2849717‑2849915
                                                                                                   |                                                                                                   
GCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGTTTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATC                                                                                                     <  SRR3722078.837090/100‑1 (MQ=60)
 CACGACCCTGCTGTCACGGATGAAGGGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCC                                                                                                    <  SRR3722078.322938/100‑1 (MQ=60)
    GACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATG                                                                                                 <  SRR3722078.353244/100‑1 (MQ=60)
    GACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATG                                                                                                 <  SRR3722078.348493/100‑1 (MQ=60)
    GACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATG                                                                                                 <  SRR3722078.1339193/100‑1 (MQ=60)
     ACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGA                                                                                                >  SRR3722078.1062806/1‑100 (MQ=60)
     ACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGA                                                                                                >  SRR3722078.1229516/1‑100 (MQ=60)
      CCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAG                                                                                               <  SRR3722078.800335/100‑1 (MQ=60)
          GCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTT                                                                                           >  SRR3722078.295259/1‑100 (MQ=60)
             GTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTT                                                                                        >  SRR3722078.788532/1‑100 (MQ=60)
             GTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTT                                                                                        >  SRR3722078.1445690/1‑100 (MQ=60)
                           GGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTT                                                                          <  SRR3722078.929773/100‑1 (MQ=60)
                            GCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTC                                                                         >  SRR3722078.1187020/1‑100 (MQ=60)
                                GGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGC                                                                     <  SRR3722078.1240656/100‑1 (MQ=60)
                                           cacgaccgccatcgtcggcagcgtcagatgtgtataagagacagCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCAC                                                          <  SRR3722078.116358/56‑1 (MQ=60)
                                           ATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCAC                                                          <  SRR3722078.1413148/100‑1 (MQ=60)
                                           ATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCAC                                                          <  SRR3722078.1386260/100‑1 (MQ=60)
                                           ATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCAC                                                          <  SRR3722078.1027522/100‑1 (MQ=60)
                                              GGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGG                                                       >  SRR3722078.858910/1‑100 (MQ=60)
                                               GCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGA                                                      >  SRR3722078.426673/1‑100 (MQ=60)
                                               GCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGA                                                      >  SRR3722078.54792/1‑100 (MQ=60)
                                                   TGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCT                                                  >  SRR3722078.830841/1‑100 (MQ=60)
                                                         AACTGCAAATTGCCAGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCAT                                            <  SRR3722078.1537484/100‑1 (MQ=60)
                                                          ACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATC                                           <  SRR3722078.840947/100‑1 (MQ=60)
                                                          ACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATC                                           >  SRR3722078.1397586/1‑100 (MQ=60)
                                                              CAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATG                                       <  SRR3722078.800621/100‑1 (MQ=60)
                                                                  TTGCCGGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGAT                                   >  SRR3722078.685764/1‑100 (MQ=60)
                                                                       GGACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGC                              >  SRR3722078.755172/1‑100 (MQ=60)
                                                                         ACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAC                            <  SRR3722078.328797/100‑1 (MQ=60)
                                                                         ACAGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCAC                            <  SRR3722078.548726/100‑1 (MQ=60)
                                                                           GGATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGG                          >  SRR3722078.581477/1‑100 (MQ=60)
                                                                            GATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACTGT                         >  SRR3722078.1181475/1‑100 (MQ=60)
                                                                            GATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGT                         >  SRR3722078.1207716/1‑100 (MQ=60)
                                                                            GATCTGCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGT                         >  SRR3722078.144779/1‑100 (MQ=60)
                                                                                 GCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGA                    <  SRR3722078.608270/100‑1 (MQ=60)
                                                                                 GCCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGA                    <  SRR3722078.711225/100‑1 (MQ=60)
                                                                                  CCTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGAT                   <  SRR3722078.276539/100‑1 (MQ=60)
                                                                                   CTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATC                  <  SRR3722078.906314/100‑1 (MQ=60)
                                                                                   CTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATC                  <  SRR3722078.723251/100‑1 (MQ=60)
                                                                                   CTGTCCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATC                  <  SRR3722078.1136956/100‑1 (MQ=60)
                                                                                     GTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCct                >  SRR3722078.471867/1‑98 (MQ=60)
                                                                                       CCGGCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGAT              <  SRR3722078.804425/100‑1 (MQ=60)
                                                                                          GCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGAC           >  SRR3722078.172379/1‑100 (MQ=60)
                                                                                          GCATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGAC           >  SRR3722078.1151325/1‑100 (MQ=60)
                                                                                            ATACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGC         >  SRR3722078.722563/1‑100 (MQ=60)
                                                                                              ACTATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCA       >  SRR3722078.398009/1‑100 (MQ=60)
                                                                                                 ATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAG    >  SRR3722078.762722/1‑100 (MQ=60)
                                                                                                 ATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAG    >  SRR3722078.1478589/1‑100 (MQ=60)
                                                                                                 ATCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAG    >  SRR3722078.916221/1‑100 (MQ=60)
                                                                                                  TCCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGC   >  SRR3722078.1395989/1‑100 (MQ=60)
                                                                                                   CCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCctgtctcttatacaca  >  SRR3722078.550771/1‑84 (MQ=60)
                                                                                                   |                                                                                                   
GCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGCATTGATTTAACTGCAAATTGCCGGACAGATCTGCCTGTCCGGCATACTATTCATGAGGTTTTTTCGGACGATATTTTTCCGGCAGTTCTGGCACCGGACGCTTGTCATCGATGAGATGACGCACGGTTAAGATCGGATGACGCCACAGCA  >  NZ_CP009273/2849717‑2849915

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