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,678,723 A→C G48G (GGT→GGG glyA ← serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,678,7230AC96.5% 181.6 / ‑1.8 57G48G (GGT→GGGglyAserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (33/22);  total (34/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

GTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTG  >  NZ_CP009273/2678638‑2678806
                                                                                     |                                                                                   
tttgttAAACAATACCAACATCCCCCCACCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCt                                                                                 <  1:426065/89‑1 (MQ=255)
   gttCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTgcg                                                                              >  1:463938/1‑90 (MQ=255)
       aaCGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTg                                                                          >  1:215427/1‑90 (MQ=255)
          gATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCAt                                                                       >  1:337034/1‑90 (MQ=255)
             atCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAc                                                                    >  1:1251382/1‑90 (MQ=255)
             atCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAc                                                                    >  2:479858/1‑90 (MQ=255)
                    tACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCggg                                                             >  1:910468/1‑90 (MQ=255)
                     aCTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGc                                                            >  1:824566/1‑90 (MQ=255)
                      cTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCt                                                           >  2:226802/1‑90 (MQ=255)
                       tCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTg                                                          >  2:248783/1‑90 (MQ=255)
                        cGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTgg                                                         >  2:1445614/1‑90 (MQ=255)
                           aaCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGGCCCCTGCGCCTGCATTACGCGCGGGCTGgtgt                                                      >  2:921248/1‑90 (MQ=255)
                           aaCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGc                                                            <  2:451366/84‑1 (MQ=255)
                           aaCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGc                                                            >  1:451366/1‑84 (MQ=255)
                           aaCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCt                                                           <  1:1308512/85‑1 (MQ=255)
                           aaCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCt                                                           >  2:1308512/1‑85 (MQ=255)
                           aaCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGgtgt                                                      >  2:747515/1‑90 (MQ=255)
                                    agtagCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCgg                                             <  2:1230391/90‑1 (MQ=255)
                                    agtagCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCgg                                             <  1:409862/90‑1 (MQ=255)
                                    agtagCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCgg                                             <  1:409863/90‑1 (MQ=255)
                                       agCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAgg                                          <  1:842943/90‑1 (MQ=255)
                                          gTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGa                                       <  1:643231/90‑1 (MQ=255)
                                          gTTTGCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGAt                                      >  2:1337167/1‑90 (MQ=255)
                                           tttGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGAt                                      <  2:1251382/90‑1 (MQ=255)
                                           tttGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGAt                                      <  2:1052527/90‑1 (MQ=255)
                                            ttGCGCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATc                                     >  2:537777/1‑90 (MQ=255)
                                            ttGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATc                                     >  1:629008/1‑90 (MQ=255)
                                                ccGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGtt                                 <  2:101559/90‑1 (MQ=255)
                                                   gATAACCTTCAGCATATTTGTTGGTCAGCTGAGGCCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGa                              <  1:921248/90‑1 (MQ=255)
                                                   gATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGa                              <  2:461953/90‑1 (MQ=255)
                                                            cAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCtt                     >  1:1400748/1‑90 (MQ=255)
                                                              gCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTcc                   <  2:239389/90‑1 (MQ=255)
                                                                 tatTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATc                                     <  2:1091488/69‑1 (MQ=255)
                                                                 tatTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGa                >  2:1295178/1‑90 (MQ=255)
                                                                 tatTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGGGTAGTTTTCGGAGGCGATc                                     >  1:1091488/1‑69 (MQ=255)
                                                                 tatTTGTTGGTCAGCTGACACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGa                >  1:182581/1‑90 (MQ=255)
                                                                  atTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAg                                                    >  2:785224/1‑53 (MQ=255)
                                                                  atTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAg                                                    <  1:785224/53‑1 (MQ=255)
                                                                      gttgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTAc           >  2:93341/1‑90 (MQ=255)
                                                                        tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCt                                                           <  2:1371716/40‑1 (MQ=255)
                                                                        tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCt                                                           >  1:1371716/1‑40 (MQ=255)
                                                                        tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAgg                                          <  2:1330339/57‑1 (MQ=255)
                                                                        tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAgg                                          >  1:1330339/1‑57 (MQ=255)
                                                                        tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGa                                       >  2:1269487/1‑60 (MQ=255)
                                                                        tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGa                                       <  1:1269487/60‑1 (MQ=255)
                                                                        tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACtt         >  1:287986/1‑90 (MQ=255)
                                                                        tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACtt         >  1:632535/1‑90 (MQ=255)
                                                                        tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACtt         >  2:707511/1‑90 (MQ=255)
                                                                        tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACtt         >  1:122900/1‑90 (MQ=255)
                                                                        tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACtt         >  1:973979/1‑90 (MQ=255)
                                                                        tgGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGGTTTCGGAGGCGATCAGGTCGATGTGCTCTTCCTGACGTACtt         >  1:298932/1‑90 (MQ=255)
                                                                            cAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTc     >  1:1006190/1‑90 (MQ=255)
                                                                             aGCTGAGACCCCTTCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTcc    <  1:1372305/90‑1 (MQ=255)
                                                                             aGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTcc    <  2:910468/90‑1 (MQ=255)
                                                                              gCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCt   >  2:298740/1‑90 (MQ=255)
                                                                              gCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCt   <  1:67736/90‑1 (MQ=255)
                                                                               cTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTg  <  2:1006190/90‑1 (MQ=255)
                                                                                     |                                                                                   
GTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTG  >  NZ_CP009273/2678638‑2678806

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ 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
CGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCAC  >  NZ_CP009273/2678625‑2678822
                                                                                                  |                                                                                                   
ggctggctctttttttttgcaactatataaacatcaccgcatCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAAC                                                                                                    <  SRR3722078.785686/58‑1 (MQ=60)
ggcggcttaagttgtttaaaagaaacacCAACATCCTCCCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAAC                                                                                                    <  SRR3722078.518924/72‑1 (MQ=60)
cgatcgctccccttttttcaaaaGATATAAACATCCTCCCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACC                                                                                                    <  SRR3722078.428620/77‑1 (MQ=60)
CGATCGATCGCCTGTTTTCAAACAATATCAACATCCTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACC                                                                                                    <  SRR3722078.471119/100‑1 (MQ=60)
      ATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCG                                                                                              >  SRR3722078.469477/1‑100 (MQ=60)
          CCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTG                                                                                          >  SRR3722078.218034/1‑100 (MQ=60)
            tgttttacaacaataccaacatccccccacCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCA                                                                                        <  SRR3722078.1060614/70‑1 (MQ=60)
             tttgttaaacaataccaacatccccccacCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCAT                                                                                       <  SRR3722078.431198/71‑1 (MQ=60)
             GTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCAT                                                                                       >  SRR3722078.341154/1‑100 (MQ=60)
                GTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAC                                                                                    >  SRR3722078.1271191/1‑100 (MQ=60)
                       GATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGG                                                                             >  SRR3722078.922754/1‑100 (MQ=60)
                        ATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGC                                                                            >  SRR3722078.835492/1‑100 (MQ=60)
                              ACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGT                                                                      >  SRR3722078.456760/1‑100 (MQ=60)
                                atccccccacCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAG                                                                    <  SRR3722078.956387/90‑1 (MQ=60)
                                acccccccacCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAG                                                                    <  SRR3722078.449142/90‑1 (MQ=60)
                                   CTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTT                                                                 <  SRR3722078.1329398/100‑1 (MQ=60)
                                          gcagcgtcagatgtgtataagagacagCCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGG                                                          <  SRR3722078.795459/73‑1 (MQ=60)
                                               GTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATC                                                     >  SRR3722078.636244/1‑100 (MQ=60)
                                                 AGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAG                                                   <  SRR3722078.414818/100‑1 (MQ=60)
                                                 AGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAG                                                   <  SRR3722078.414817/100‑1 (MQ=60)
                                                    AGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTC                                                <  SRR3722078.854187/100‑1 (MQ=60)
                                                       cagatgtgtataagagacagGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGAT                                             <  SRR3722078.1289630/80‑1 (MQ=60)
                                                       GTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGAT                                             <  SRR3722078.650622/100‑1 (MQ=60)
                                                               GGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTT                                     >  SRR3722078.1423436/1‑100 (MQ=60)
                                                                GATAACCTTCAGCATATTTGTTGGTCAGCTGAGGCCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTC                                    <  SRR3722078.933840/100‑1 (MQ=60)
                                                                    ACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGGGTAGTTTTCGGAGGCGATCAGTTCGATGTctgtctcttat                                >  SRR3722078.1108055/1‑89 (MQ=60)
                                                                    ACCTTCAGCATATTTGTTGGTCAGCTGACACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGA                                >  SRR3722078.184744/1‑100 (MQ=60)
                                                                           GCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTctgtctcttatacacatctgacgctgacgacgatggcggg                         >  SRR3722078.1393808/1‑60 (MQ=60)
                                                                           GCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCctgtctcttatacacatctgacg                         >  SRR3722078.1351630/1‑77 (MQ=60)
                                                                           GCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTT                         >  SRR3722078.639807/1‑100 (MQ=60)
                                                                           GCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTT                         >  SRR3722078.987953/1‑100 (MQ=60)
                                                                           GCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTT                         >  SRR3722078.124331/1‑100 (MQ=60)
                                                                           GCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTT                         >  SRR3722078.291499/1‑100 (MQ=60)
                                                                           GCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGGTTTCGGAGGCGATCAGGTCGATGTGCTCTTCCTGACGTACTT                         >  SRR3722078.302589/1‑100 (MQ=60)
                                                                               ATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTC                     >  SRR3722078.1020947/1‑100 (MQ=60)
                                                                                          AGCTGAGACCCCTTCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAG          <  SRR3722078.1394416/100‑1 (MQ=60)
                                                                                           GCTGAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGC         <  SRR3722078.68529/100‑1 (MQ=60)
                                                                                           GCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGC         >  SRR3722078.1392826/1‑100 (MQ=60)
                                                                                              GAGACCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTG      >  SRR3722078.570339/1‑100 (MQ=60)
                                                                                                  CCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCAC  >  SRR3722078.200099/1‑100 (MQ=60)
                                                                                                  CCCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCAC  <  SRR3722078.337805/100‑1 (MQ=60)
                                                                                                  |                                                                                                   
CGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGCTGGTGTAGTTTTCGGAGGCGATCAGTTCGATGTGCTCTTCCTGACGTACTTTTTCCTGCTCCATAGCCTGCCAC  >  NZ_CP009273/2678625‑2678822

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