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
A5 F25 I1 R1 61 74.4 3161000 98.6% 3116746 146.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 6,016,700 T→C intergenic (+90/+20) xpt → / ← PP_5266 xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC6,016,7000TC96.7% 99.0 / NA 30intergenic (+90/+20)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (20/9);  minor base A (1/0);  total (21/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAA  >  NC_002947_CJ‑RC/6016584‑6016845
                                                                                                                    |                                                                                                                                                 
ggggttgtggcCCGCTTCACCCAATAACCCCCACCCGGGCTGCGGGGGCCCGCCGGCGCCCCCTGGTCCCGCACTGGGCTCCACGGCTGGGCGGTACGTGGGGTCGGGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTgcc                                                                                                                   <  6:94176/144‑1 (MQ=255)
                          aaCGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACCGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTg                                                                                         <  2:350478/149‑1 (MQ=255)
                                           ggggccACCCCCGCCCCCACCGGCCCCGCAGTGCTCTCCACCGCTGGGCCTAACGTGGGGGCGGGGGCGTGCCCGTGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTaa                                                                        <  1:8481/144‑1 (MQ=255)
                                                 aCGCCCGCTCCCACAGGTGCAGCACAGGTCTCGCCCGCGGTGGCGTACCTGTGGGAGCGGGCGTGCCCGCGTATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACCACCAGCGCTGGTAAAGTTGTTCTTTAAGCCCCt                                                                  >  7:135457/1‑149 (MQ=38)
                                                       gCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCt                                                            >  6:297110/1‑149 (MQ=255)
                                                            cacaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                                        >  2:237369/1‑148 (MQ=255)
                                                            cacaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                                        >  2:77127/1‑148 (MQ=255)
                                                            cacaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                                        >  1:83490/1‑148 (MQ=255)
                                                            cacaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                                        >  5:5091/1‑148 (MQ=255)
                                                            cacaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGccc                                                        >  6:92495/1‑148 (MQ=21)
                                                             acaGGTGCAGCACAGGTCTCCAACGCTCTGCCGTACCTGTAGGAGCGGGCGTGCCCGCGAATGGCCCCAACCCTTTTATCTCCTTGCCGCCAGCCCCGCCAACAACAAGCGCTCGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                                                       >  4:263137/1‑148 (MQ=255)
                                                             acaGGTACGGCACAGCGTTCGAGACCTGTGCTGCACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                                                       >  6:370356/1‑148 (MQ=21)
                                                               aGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAAc                                                    >  5:8225/1‑149 (MQ=255)
                                                                       cacaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCg                                            >  6:170044/1‑149 (MQ=255)
                                                                        acaGGTATCGAACGCTGTGCCGTACCTGTGGCAGCGGGCGTGCCAGCGAATGGACGCAAAACTTTTATCTGTTTGACGCCAGCCACGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGCCTGCCCAAACTGCATCCGc                                           >  5:187894/1‑149 (MQ=255)
                                                                               tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATc                                              >  7:359944/1‑139 (MQ=255)
                                                                               tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATc                                              <  8:359944/139‑1 (MQ=255)
                                                                                 gAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGt                                    <  5:170044/147‑1 (MQ=255)
                                                                                 gAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTgg                                  >  6:163770/1‑149 (MQ=255)
                                                                                           ccGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATAc                         <  5:297110/148‑1 (MQ=255)
                                                                                            cGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCTACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCTAACTGCATCCGCGCCAGGTGGCTTTCATAccc                       <  6:187894/149‑1 (MQ=255)
                                                                                                 cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttc                                                                                <  7:48550/87‑1 (MQ=255)
                                                                                                 cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGttc                                                                                >  8:48550/1‑87 (MQ=255)
                                                                                                 cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCg                  >  6:397234/1‑149 (MQ=255)
                                                                                                 cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCg                  >  7:160987/1‑149 (MQ=255)
                                                                                                 cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCg                  >  8:151276/1‑149 (MQ=255)
                                                                                                         gCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGt          <  8:234967/149‑1 (MQ=255)
                                                                                                              cGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATc     >  3:33587/1‑149 (MQ=255)
                                                                                                                tGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATcc    >  5:188147/1‑148 (MQ=255)
                                                                                                                 gCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCaa  >  5:238462/1‑149 (MQ=255)
                                                                                                                    |                                                                                                                                                 
GGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAA  >  NC_002947_CJ‑RC/6016584‑6016845

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGC  >  NC_002947_CJ‑RC/6016609‑6016848
                                                                                           |                                                                                                                                                    
acccccccaccgggctccggggccccccccgccccaactgggcccccacggttcgccacgggtggccggtacggggggggggGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGG                                                                                           <  NB501016:282:H2V7GAFX2:3:11609:18879:8384/69‑1 (MQ=60)
     cccccggggcggggggcccccccccccccccgcggggccccccgggtgccccccggtggcgcctgtgggggggggggggggGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAG                                                                                      <  NB501016:282:H2V7GAFX2:3:21404:13116:4816/70‑1 (MQ=60)
             ggcgcggccccgccccccccccctcgggccccccccgggccccccgggcgcgccccccgggggggggggggggGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTT                                                                              <  NB501016:282:H2V7GAFX2:4:11405:23486:9558/78‑1 (MQ=60)
                 cggggccacccccgcccccaccggccccgcagtgctctccaccgctgggcctaacgtgggggcggGGGCGTGCCCGTGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAG                                                                          <  NB501016:282:H2V7GAFX2:1:11102:13993:9511/86‑1 (MQ=60)
                       CACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGCCCGCGGTGGCGTACCTGTGGGAGCGGGCGTGCCCGCGTATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACCACCAGCGCTGGTAAAGTTGTTCTTTAAGCCCCTG                                                                    >  NB501016:282:H2V7GAFX2:4:11512:21723:11248/1‑151 (MQ=60)
                         cgcccgctcccacaggtacggcccagcgttcgagacctgtggtgccCCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCCACCCGTTTATCTGGTTTCCGCCCGCCCCGCCAACCACAAGCGCTGGGAAAGTTCTTCTTTAAGCCCCTGCG                                                                  >  NB501016:282:H2V7GAFX2:3:11410:6686:10956/47‑151 (MQ=60)
                                  CCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCC                                                          >  NB501016:282:H2V7GAFX2:1:11203:18006:6186/1‑150 (MQ=60)
                                  CCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCC                                                          >  NB501016:282:H2V7GAFX2:3:11401:15683:19041/1‑150 (MQ=60)
                                    cgggcccccccatgtcacccccggggggcccctacgggggggggggggcgGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAAC                                                       <  NB501016:282:H2V7GAFX2:4:21603:21995:5658/101‑1 (MQ=60)
                                     CAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACT                                                      >  NB501016:282:H2V7GAFX2:3:11402:22521:10864/1‑151 (MQ=60)
                                              AACAGGTATCGAACGCTGTGCCGTACCTGTGGCAGCGGGCGTGCCAGCGAATGGACGCAAAACTTTTATCTGTTTGACGCCAGCCACGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGCCTGCCCAAACTGCATCCGCC                                             >  NB501016:282:H2V7GAFX2:3:11610:1296:7383/1‑151 (MQ=60)
                                                     CTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCC                                                >  NB501016:282:H2V7GAFX2:4:21604:16112:1711/1‑141 (MQ=60)
                                                       cgaccggggtgcggaccgggggggcggggggGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTG                                      <  NB501016:282:H2V7GAFX2:4:11502:1536:18506/118‑1 (MQ=60)
                                                       CGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTG                                      <  NB501016:282:H2V7GAFX2:3:11607:5820:6642/149‑1 (MQ=60)
                                                                 GCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACC                           <  NB501016:282:H2V7GAFX2:3:21506:13722:15384/150‑1 (MQ=60)
                                                                       CCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGG                    >  NB501016:282:H2V7GAFX2:4:11605:20676:5648/1‑151 (MQ=60)
                                                                       CCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCT                                                                                  <  NB501016:282:H2V7GAFX2:4:11409:14468:5751/89‑1 (MQ=60)
                                                                                    GCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCC       >  NB501016:282:H2V7GAFX2:2:11107:5251:7429/1‑151 (MQ=60)
                                                                                      GTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCA      >  NB501016:282:H2V7GAFX2:3:11610:15375:8261/1‑150 (MQ=60)
                                                                                       TGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAAC    >  NB501016:282:H2V7GAFX2:3:21407:20081:16382/1‑151 (MQ=60)
                                                                                        GCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACG   <  NB501016:282:H2V7GAFX2:3:11505:20842:15965/151‑1 (MQ=60)
                                                                                         CCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGC  >  NB501016:282:H2V7GAFX2:4:11406:19803:16493/1‑151 (MQ=60)
                                                                                         CCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGC  >  NB501016:282:H2V7GAFX2:3:11412:4034:3661/1‑151 (MQ=60)
                                                                                           |                                                                                                                                                    
AAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGC  >  NC_002947_CJ‑RC/6016609‑6016848

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: