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 F5 I97 R1 275 44.0 1733530 96.5% 1672856 122.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 923,801 A→G V284A (GTC→GCC)  cydC ← cysteine/glutathione ABC transporter ATP‑binding protein/permease CydC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273923,8010AG100.0% 165.9 / NA 44V284A (GTC→GCC) cydCcysteine/glutathione ABC transporter ATP‑binding protein/permease CydC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/25);  total (19/25)

TATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGA  >  NZ_CP009273/923667‑923935
                                                                                                                                      |                                                                                                                                      
tATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCaaa                                                                                                                                    <  1:277172/139‑1 (MQ=255)
   cAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAg                                                                                                                                 >  1:378990/1‑139 (MQ=255)
              aCCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCaaa                                                                                                                                    >  1:100123/1‑125 (MQ=255)
              aCCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCaaa                                                                                                                                    <  2:100123/125‑1 (MQ=255)
                cTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAAcgcg                                                                                                                    <  1:127115/139‑1 (MQ=255)
                cTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAAcgcg                                                                                                                    <  1:841917/139‑1 (MQ=255)
                cTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAAcgcg                                                                                                                    <  2:357043/139‑1 (MQ=255)
                cTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAAcgcg                                                                                                                    <  2:380739/139‑1 (MQ=255)
                        tttGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGtt                                                                                                             <  2:290065/138‑1 (MQ=255)
                                         gagaTACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAAcgc                                                                                           >  1:160233/1‑139 (MQ=255)
                                                gtacgGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGa                                                                                    >  1:653700/1‑139 (MQ=255)
                                                    ggCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGcc                                                                                <  2:493106/139‑1 (MQ=255)
                                                    ggCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGc                                                                                 >  1:59196/1‑138 (MQ=255)
                                                         agGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAAcgcg                                                                                                                    >  1:618464/1‑98 (MQ=255)
                                                         agGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAAcgcg                                                                                                                    <  2:618464/98‑1 (MQ=255)
                                                         agGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCca                                                                           >  1:626151/1‑139 (MQ=255)
                                                         agGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCca                                                                           <  2:609104/139‑1 (MQ=255)
                                                                     gCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCa                                                                               >  1:292821/1‑123 (MQ=255)
                                                                     gCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCa                                                                               <  2:292821/123‑1 (MQ=255)
                                                                     gCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAgg                                                                 <  1:344398/137‑1 (MQ=255)
                                                                     gCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGAt                                                               >  2:42489/1‑139 (MQ=255)
                                                                      ccccAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGAt                                                               <  2:331147/138‑1 (MQ=255)
                                                                                tGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAAcg                                                    >  2:231792/1‑139 (MQ=255)
                                                                                  aaaTGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACg                                                                                  <  1:606455/107‑1 (MQ=255)
                                                                                  aaaTGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACg                                                                                  >  2:606455/1‑107 (MQ=255)
                                                                                  aaaTGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAAcgcg                                                  <  2:59196/139‑1 (MQ=255)
                                                                                  aaaTGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAAcg                                                    <  2:153401/137‑1 (MQ=255)
                                                                                         cccGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGc                                                                       <  1:637663/111‑1 (MQ=255)
                                                                                         cccGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGc                                                                       >  2:637663/1‑111 (MQ=255)
                                                                                          ccGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCATACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGag                                          <  2:240488/139‑1 (MQ=255)
                                                                                          ccGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGag                                          <  2:653700/139‑1 (MQ=255)
                                                                                          ccGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGag                                          <  2:378990/139‑1 (MQ=255)
                                                                                          ccGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGag                                          <  2:589178/139‑1 (MQ=255)
                                                                                            gTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGAt                                                               >  1:824416/1‑116 (MQ=255)
                                                                                            gTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGAt                                                               <  2:824416/116‑1 (MQ=255)
                                                                                                 tGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCAt                                    >  2:827529/1‑138 (MQ=255)
                                                                                                   gTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTa                                  <  2:626151/138‑1 (MQ=255)
                                                                                                    tGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATc                                >  2:29058/1‑139 (MQ=255)
                                                                                                          tGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGc                          >  1:684368/1‑139 (MQ=255)
                                                                                                              tCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGAc                       >  1:819625/1‑138 (MQ=255)
                                                                                                                      ggCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCa                                                                                                        >  2:503525/1‑49 (MQ=255)
                                                                                                                      ggCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCa                                                                                                        <  1:503525/49‑1 (MQ=255)
                                                                                                                           aCGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGtt          <  1:295664/138‑1 (MQ=255)
                                                                                                                                cAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGAtt    >  2:551156/1‑139 (MQ=255)
                                                                                                                                  gaaggCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGa  <  1:231792/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
TATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGA  >  NZ_CP009273/923667‑923935

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGC  >  NZ_CP009273/923659‑923946
                                                                                                                                              |                                                                                                                                                 
AGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGG                                                                                                                                           >  SRR3722120.384575/1‑151 (MQ=60)
       GTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTA                                                                                                                                    <  SRR3722120.281167/151‑1 (MQ=60)
           CAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTActgt                                                                                                                                >  SRR3722120.101462/1‑147 (MQ=60)
                       CCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGC                                                                                                                    <  SRR3722120.128834/151‑1 (MQ=60)
                       CCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGC                                                                                                                    <  SRR3722120.853242/151‑1 (MQ=60)
                                    cgatgtttttttacaagcagaaggcggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAAC                                                                                                        <  SRR3722120.510727/71‑1 (MQ=60)
                                      CCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCC                                                                                                     >  SRR3722120.162425/1‑151 (MQ=60)
                                             GTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGAC                                                                                              >  SRR3722120.662696/1‑151 (MQ=60)
                                                 GAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCC                                                                                           >  SRR3722120.59922/1‑150 (MQ=60)
                                                      ACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCctgtctcttatacacatctccgagcccacg                                                                                      >  SRR3722120.627067/1‑120 (MQ=60)
                                                      ACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCAC                                                                                     >  SRR3722120.634836/1‑151 (MQ=60)
                                                         cggagatgtgtataagagacagCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGC                                                                                  <  SRR3722120.614874/129‑1 (MQ=60)
                                                                  GGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCctgtct                                                                         >  SRR3722120.297014/1‑145 (MQ=60)
                                                                    gatgtgtataagagacagGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCAC                                                                       <  SRR3722120.646459/133‑1 (MQ=60)
                                                                            TGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAA                                                                 <  SRR3722120.349380/149‑1 (MQ=60)
                                                                                         GAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCtgtctcttatac                                                  >  SRR3722120.835499/1‑139 (MQ=60)
                                                                                                       ACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCG                                    >  SRR3722120.693744/1‑151 (MQ=60)
                                                                                                           GTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACA                                 >  SRR3722120.830652/1‑150 (MQ=60)
                                                                                                                                  AACGCGCAGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGG          <  SRR3722120.299888/150‑1 (MQ=60)
                                                                                                                                         AGAAGGCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGC  <  SRR3722120.235063/151‑1 (MQ=60)
                                                                                                                                              |                                                                                                                                                 
AGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGC  >  NZ_CP009273/923659‑923946

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