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
A0 F0 I1 R1 52 69.5 2964196 98.4% 2916768 145.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 1,126,647 +C intergenic (‑101/+48) tdcG‑II ← / ← gcvP‑I L‑serine dehydratase/glycine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC1,126,6451.C100.0% 89.6 / NA 26intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (20/6);  total (20/6)

TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  NC_002947_CJ‑RC/1126578‑1126788
                                                                    |                                                                                                                                               
ttGAGCCCTGTAGGCGCGGGTTTATCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACggg                                                                 >  6:352603/1‑149 (MQ=255)
                  ggTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTg                                               >  7:353841/1‑149 (MQ=255)
                         ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGATGGGCAGGCGCACACCAGGTTGCGGTCGc                                        >  5:126532/1‑149 (MQ=255)
                         ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                        >  5:284242/1‑149 (MQ=255)
                         ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                        >  3:252086/1‑149 (MQ=255)
                           gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                       >  8:251000/1‑148 (MQ=255)
                           gcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGcc                                       >  8:142022/1‑148 (MQ=255)
                                   ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                              >  1:330642/1‑149 (MQ=255)
                                    gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACg                                                                   >  2:77695/1‑111 (MQ=255)
                                    gAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACg                                                                   <  1:77695/111‑1 (MQ=255)
                                     ttttataCCCCCCATCCCATTGGGGGTTAAACCCGCTCCTAAAAATACCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGt                            <  7:108501/146‑1 (MQ=255)
                                      ggTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTc                           >  3:96007/1‑149 (MQ=255)
                                         ataGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCCAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCAAGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAc                        >  1:354920/1‑149 (MQ=39)
                                             ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                    >  3:324958/1‑149 (MQ=255)
                                              caccaTCGAATTGGGGGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc                   <  2:264602/149‑1 (MQ=255)
                                               accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc                  >  2:40723/1‑149 (MQ=255)
                                                 caTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGa                >  2:153605/1‑149 (MQ=255)
                                                      cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggc           >  7:297166/1‑149 (MQ=255)
                                                       aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg          >  8:163113/1‑149 (MQ=255)
                                                         tCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc        >  2:80351/1‑149 (MQ=255)
                                                          cgcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc        <  5:352603/148‑1 (MQ=255)
                                                           gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCggcggc        >  6:333664/1‑147 (MQ=255)
                                                           gcgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCa      <  7:163113/149‑1 (MQ=255)
                                                              ggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTa   <  1:153605/149‑1 (MQ=255)
                                                              ggTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTa   >  5:360590/1‑149 (MQ=255)
                                                               gTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAc  >  2:299508/1‑149 (MQ=255)
                                                                    |                                                                                                                                               
TTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  NC_002947_CJ‑RC/1126578‑1126788

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGC  >  NC_002947_CJ‑RC/1126541‑1126792
                                                                                                           |                                                                                                                                                 
cccgccggcccggccgcgggggggtgtgggtgccccgcccccggtcggggggggggtccccccgccctcggggttcccccccccccccccccgggggctccaacccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGG                                                                                                            <  NB501016:282:H2V7GAFX2:2:11303:26107:14389/41‑1 (MQ=60)
   tggtggccccccgtcgttgtgggtggctgaaacagggcgcggtgcgccggggttttcggtctccttggcttgctacccccccctggcgctggggtataaacccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGG                                                                                                        <  NB501016:282:H2V7GAFX2:2:11204:6176:4906/45‑1 (MQ=60)
           cctctgcccggtttttggggcttatttcccccctgtgggggggggtacccccccctatggggtttaatcccccccctcccctggcgcttaaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCG                                                                                             <  NB501016:282:H2V7GAFX2:4:11401:26402:5517/56‑1 (MQ=60)
              ccccccgggtggggggggccactcccccccgtggggggggggttaccccccctagggcggcttaacccccccccccccccgggggtaaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCC                                                                                          <  NB501016:282:H2V7GAFX2:2:11304:15742:4700/59‑1 (MQ=60)
                gccctggttgtgggtggcacttgcccccgggtgggggggttttccccggctaggggggtgtaaacccccccaccccttgggggttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTG                                                                                        <  NB501016:282:H2V7GAFX2:4:11509:12781:11164/61‑1 (MQ=60)
                             cctccccccacccccgtctcgggcggctaccccccgctatgcggttttaacccccccaaccccttgccctttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGAT                                                                            <  NB501016:282:H2V7GAFX2:3:21606:13043:18384/73‑1 (MQ=60)
                                 cacttggcgccgtgcggggggggttaccccggctaatgcgttgtaaccccccccagcccagccggtttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACG                                                                       <  NB501016:282:H2V7GAFX2:3:21606:6058:19064/78‑1 (MQ=60)
                                                gggggggtttaccccggcttcggggctctaacccccacaccccattggggttaaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCA                                                        <  NB501016:282:H2V7GAFX2:3:21508:16606:8181/93‑1 (MQ=60)
                                                      GGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGC                                                  >  NB501016:282:H2V7GAFX2:4:21606:20662:18881/1‑151 (MQ=60)
                                                         ttcccccggttgggggtaggaatccccccctggctttgcggggaaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGG                                                <  NB501016:282:H2V7GAFX2:1:11105:6255:12667/101‑1 (MQ=60)
                                                         ttaccccgctatggggttttaaaccccaccctgcattgggggtttaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACTCCAGGTTGCGG                                                <  NB501016:282:H2V7GAFX2:4:11610:20971:5927/101‑1 (MQ=60)
                                                         atccccggtcttgggggtttaaaccccccccccctttgggggttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGG                                                <  NB501016:282:H2V7GAFX2:4:11604:8045:15902/101‑1 (MQ=60)
                                                         acaccgcgatatggcgtattaatcccccaaatccattgcggttaaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGTCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGG                                                <  NB501016:282:H2V7GAFX2:1:21206:2129:6099/101‑1 (MQ=60)
                                                             CCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGATGGGCAGGCGCACACCAGGTTGCGGTCGCC                                           >  NB501016:282:H2V7GAFX2:3:11601:25368:8768/1‑151 (MQ=60)
                                                             CCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCC                                           >  NB501016:282:H2V7GAFX2:2:21203:22292:9949/1‑151 (MQ=60)
                                                             CCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCC                                           >  NB501016:282:H2V7GAFX2:3:21508:8052:6837/1‑151 (MQ=60)
                                                                 ataaggggttgtaaacccacccaccccatgggggttaaaacCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAAC                                       <  NB501016:282:H2V7GAFX2:3:21503:17144:9316/110‑1 (MQ=60)
                                                                       gggtgtaaaccccccccccccatgcggggtaaaacCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTG                                 <  NB501016:282:H2V7GAFX2:4:11607:12710:18066/116‑1 (MQ=60)
                                                                       CGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTG                                 >  NB501016:282:H2V7GAFX2:1:21305:13379:10698/1‑151 (MQ=60)
                                                                        GGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGG                                                                      <  NB501016:282:H2V7GAFX2:1:11203:22257:6068/113‑1 (MQ=60)
                                                                          tttttcccccccccgcgccggggcgttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTC                               <  NB501016:282:H2V7GAFX2:2:11202:17088:8576/118‑1 (MQ=60)
                                                                          tgtaaccccccccaccgatggggtgttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTC                               <  NB501016:282:H2V7GAFX2:3:11402:1724:9729/118‑1 (MQ=60)
                                                                          AGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCG                              >  NB501016:282:H2V7GAFX2:2:11207:24044:16225/1‑151 (MQ=60)
                                                                             TATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCCAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCAAGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACC                           >  NB501016:282:H2V7GAFX2:1:21310:3500:10569/1‑151 (MQ=60)
                                                                                 GCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGC                       >  NB501016:282:H2V7GAFX2:2:21306:1876:14039/1‑151 (MQ=60)
                                                                                   CAACATCGCATTCGAGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGATCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCG                     >  NB501016:282:H2V7GAFX2:3:11606:26574:14369/1‑151 (MQ=60)
                                                                                          GCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGCCCGGCG              >  NB501016:282:H2V7GAFX2:4:21508:25225:2438/1‑151 (MQ=60)
                                                                                              TCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCC           <  NB501016:282:H2V7GAFX2:3:21609:7447:5520/150‑1 (MQ=60)
                                                                                               CGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAG         <  NB501016:282:H2V7GAFX2:4:11607:25319:5928/151‑1 (MQ=60)
                                                                                                  GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCTGGGGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC      <  NB501016:282:H2V7GAFX2:1:11306:20369:2685/151‑1 (MQ=60)
                                                                                                  GGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC      >  NB501016:282:H2V7GAFX2:3:21610:6759:18539/1‑151 (MQ=60)
                                                                                                     aaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTG   <  NB501016:282:H2V7GAFX2:4:11408:24195:18378/146‑1 (MQ=60)
                                                                                                      aaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGC  <  NB501016:282:H2V7GAFX2:2:11107:25714:9298/147‑1 (MQ=60)
                                                                                                           |                                                                                                                                                 
TGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGC  >  NC_002947_CJ‑RC/1126541‑1126792

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