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,882,413 A→C V349G (GTT→GGT)  cysI ← assimilatory sulfite reductase (NADPH) hemoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,882,4130AC98.1% 170.9 / ‑4.2 52V349G (GTT→GGT) cysIassimilatory sulfite reductase (NADPH) hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (26/25);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGC  >  NZ_CP009273/2882326‑2882498
                                                                                       |                                                                                     
ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                     >  1:1106299/1‑90 (MQ=255)
ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                     <  2:1392305/90‑1 (MQ=255)
ccGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                     >  1:311248/1‑90 (MQ=255)
  ggTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCa                                                                                   <  1:446811/90‑1 (MQ=255)
       tCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCaa                                                                              <  1:1326323/90‑1 (MQ=255)
        cAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAAt                                                                             >  1:916720/1‑90 (MQ=255)
          gCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGGGCCAGTTATCATCAATGCCCTTACCCCAGCCAATAc                                                                           >  1:744189/1‑90 (MQ=255)
            ggACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGa                                                                         >  2:1473969/1‑90 (MQ=255)
            ggACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGa                                                                         >  2:197782/1‑90 (MQ=255)
                 cgcCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGcc                                                                    <  2:540901/90‑1 (MQ=255)
                 cgcCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGcc                                                                    <  2:1443466/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  2:1245599/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  1:89263/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  1:468380/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  2:188291/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  2:267520/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  2:509360/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  2:618901/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  2:916720/90‑1 (MQ=255)
                              aGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  1:1048928/90‑1 (MQ=255)
                               ggATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAGTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       <  2:173541/89‑1 (MQ=255)
                               ggATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAGTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAAc                                                       >  1:173541/1‑89 (MQ=255)
                               ggATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACt                                                      >  1:1158862/1‑90 (MQ=255)
                                 aTGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCa                                                    >  1:922861/1‑90 (MQ=255)
                                 aTGCGACCATTTTCGATAAACAGCGTCAGGGGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCa                                                    >  1:1445727/1‑90 (MQ=255)
                                   gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                     >  1:1062458/1‑55 (MQ=255)
                                   gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAccc                                                                                     <  2:1062458/55‑1 (MQ=255)
                                   gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                                  <  2:1333493/90‑1 (MQ=255)
                                   gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                                  >  1:906744/1‑90 (MQ=255)
                                   gCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCata                                                  >  1:634879/1‑90 (MQ=255)
                                    cGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCatat                                                 >  1:284171/1‑90 (MQ=255)
                                     gACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATg                                                >  2:162190/1‑90 (MQ=255)
                                             tCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTAACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc                                        <  2:1410061/90‑1 (MQ=255)
                                                  aaaCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTc                                   <  2:532449/90‑1 (MQ=255)
                                                  aaaCAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTc                                   <  2:284171/90‑1 (MQ=255)
                                                       gCGTCAGGTGCCAGTTATCATCAATGCCCTTAcccc                                                                                    <  2:210131/36‑1 (MQ=37)
                                                       gCGTCAGGTGCCAGTTATCATCAATGCCCTTAcccc                                                                                    >  1:210131/1‑36 (MQ=37)
                                                       gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAg                                                                                  >  2:1420539/1‑38 (MQ=38)
                                                       gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAg                                                                                  <  1:1420539/38‑1 (MQ=38)
                                                       gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAAtt                              >  2:871746/1‑90 (MQ=255)
                                                       gCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAAtt                              >  2:262277/1‑90 (MQ=255)
                                                        cGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAAttt                             >  2:672618/1‑90 (MQ=255)
                                                           cAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGAt                          >  2:231294/1‑90 (MQ=255)
                                                            aGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATc                         >  2:946484/1‑90 (MQ=255)
                                                             ggTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATcc                        >  2:911396/1‑90 (MQ=255)
                                                                  cAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc                                        >  2:475029/1‑69 (MQ=255)
                                                                  cAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATc                                        <  1:475029/69‑1 (MQ=255)
                                                                   aGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgcgc                  <  2:341467/90‑1 (MQ=255)
                                                                   aGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgcgc                  <  1:911396/90‑1 (MQ=255)
                                                                   aGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCcgcgc                  <  1:112213/90‑1 (MQ=255)
                                                                         catcaATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCt            <  1:271295/90‑1 (MQ=255)
                                                                         catcaATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCt            <  1:1210931/90‑1 (MQ=255)
                                                                              aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt       >  2:241418/1‑90 (MQ=255)
                                                                              aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt       >  2:416263/1‑90 (MQ=255)
                                                                              aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt       >  2:503816/1‑90 (MQ=255)
                                                                              aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt       >  1:372300/1‑90 (MQ=255)
                                                                              aTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACt       >  2:1070749/1‑90 (MQ=255)
                                                                                 cccTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTcc    <  2:1445727/90‑1 (MQ=255)
                                                                                   cTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGc  >  1:1443269/1‑90 (MQ=255)
                                                                                       |                                                                                     
CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGC  >  NZ_CP009273/2882326‑2882498

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ 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
ATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCA  >  NZ_CP009273/2882314‑2882511
                                                                                                   |                                                                                                  
ATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAC                                                                                                    >  SRR3722078.190522/1‑100 (MQ=60)
  CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGCCCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                                  <  SRR3722078.962294/100‑1 (MQ=60)
  CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                                  >  SRR3722078.315062/1‑100 (MQ=60)
  CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                                  >  SRR3722078.1123180/1‑100 (MQ=60)
  CTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCC                                                                                                  <  SRR3722078.1206032/100‑1 (MQ=60)
   TCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCC                                                                                                 <  SRR3722078.1388874/100‑1 (MQ=60)
          GGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAAT                                                                                          >  SRR3722078.929184/1‑100 (MQ=60)
            CCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGGGCCAGTTATCATCAATGCCCTTACCCCAGCCAATAC                                                                                        >  SRR3722078.753617/1‑100 (MQ=60)
              GGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGA                                                                                      <  SRR3722078.452154/100‑1 (MQ=60)
               cgaccgccatcgtcggcagcgtcagatgtgtataagagacagTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGAT                                                                                     <  SRR3722078.1443604/58‑1 (MQ=60)
                   TCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCC                                                                                 <  SRR3722078.1347529/100‑1 (MQ=60)
                                 GGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAGTGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACT                                                                   >  SRR3722078.175576/1‑100 (MQ=60)
                                 GGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACT                                                                   >  SRR3722078.1176803/1‑100 (MQ=60)
                                   ATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCA                                                                 >  SRR3722078.935496/1‑100 (MQ=60)
                                   ATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGGGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCA                                                                 >  SRR3722078.1469291/1‑100 (MQ=60)
                                     AATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCctgtctcttatacacatctgacgctgccga                                                               >  SRR3722078.675982/1‑70 (MQ=60)
                                     AATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACctgtctcttatacacatctgacgct                                                               >  SRR3722078.1078453/1‑75 (MQ=60)
                                     AATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATA                                                               >  SRR3722078.918943/1‑100 (MQ=60)
                                     AATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATA                                                               >  SRR3722078.642174/1‑100 (MQ=60)
                                      ATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATAT                                                              >  SRR3722078.287638/1‑100 (MQ=60)
                                          AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                          <  SRR3722078.90315/100‑1 (MQ=60)
                                          AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                          <  SRR3722078.473955/100‑1 (MQ=60)
                                          AGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGAC                                                          <  SRR3722078.1064647/100‑1 (MQ=60)
                                                         ataagagACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATT                                           <  SRR3722078.480663/93‑1 (MQ=60)
                                                         TCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTACCCCAGCCAATACGctgtctcttatacacatctgacgctgccgacgatggcggtcgtg                                           >  SRR3722078.212670/1‑56 (MQ=60)
                                                                               AGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCAC                     <  SRR3722078.923710/100‑1 (MQ=60)
                                                                               AGTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCAC                     <  SRR3722078.113538/100‑1 (MQ=60)
                                                                                GTTATCATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACT                    >  SRR3722078.376856/1‑100 (MQ=60)
                                                                                     CATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGC               <  SRR3722078.274606/100‑1 (MQ=60)
                                                                                     CATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGC               >  SRR3722078.1466783/1‑100 (MQ=60)
                                                                                     CATCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGC               <  SRR3722078.1229965/100‑1 (MQ=60)
                                                                                       TCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTT             >  SRR3722078.626136/1‑100 (MQ=60)
                                                                                       TCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTT             >  SRR3722078.150976/1‑100 (MQ=60)
                                                                                       TCAATGCCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTT             >  SRR3722078.1036305/1‑100 (MQ=60)
                                                                                             CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG       >  SRR3722078.551659/1‑100 (MQ=60)
                                                                                             CCCTTACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACG       >  SRR3722078.1513799/1‑100 (MQ=60)
                                                                                                  ACCCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCA  >  SRR3722078.814891/1‑100 (MQ=60)
                                                                                                   |                                                                                                  
ATCTCCAGCAGGCCGGTTTTCAGCGGACGCGCCGGATAATCAAGGATGCGACCATTTTCGATAAACAGCGTCAGGTGCCAGTTATCATCAATGCCCTTAACCCAGCCAATACGATCGCCTCGTCCGGTGAACTCATATGGACGGATCGGTTCAAATTTGATCCCCGCGCGACGCTCCACTTCCGCTTTAAACGTCTCA  >  NZ_CP009273/2882314‑2882511

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