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
A3 F13 I0 R1 211 68.8 2930004 98.2% 2877263 146.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947_CJ‑RC 6,016,700 T→C 100% 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,7000TC100.0% 105.5 / ‑6.8 31intergenic (+90/+20)xpt/PP_5266xanthine phosphoribosyltransferase/acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (24/6);  minor base A (1/0);  total (25/6)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAA  >  NC_002947_CJ‑RC/6016608‑6016845
                                                                                            |                                                                                                                                                 
tAAACGCGACACCTGCTGCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCTAACGCTGTGCCGTACCTGGGGGAGGGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAgcgc                                                                                           <  6:345098/149‑1 (MQ=38)
                           gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCCACCCTTTTATCTGGTTTCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTg                                                                 >  7:352078/1‑148 (MQ=255)
                           gCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCGTTTATCTGGTTGCCGCCCGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGc                                                                >  5:202779/1‑149 (MQ=255)
                            cccGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCg                                                               >  1:33246/1‑149 (MQ=255)
                               gCTCCCACAGGAGCAGCACATGTCTCGAACGCTGTGCCGTACCTGTGGCAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCACGCCACGCCAACAACAAGCGCTGGTAAAGATGTTCTTTAAGCACCTGCGGca                                                            >  4:240354/1‑148 (MQ=255)
                                  cccACTGGTGAAGCAAAGGTTTCTAACGGTGTGCGGAACGTGTGGTTGGGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                                                          <  8:352078/148‑1 (MQ=255)
                                     acaGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGcccc                                                       >  2:214867/1‑148 (MQ=255)
                                               cacaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCg                                            >  1:332741/1‑149 (MQ=255)
                                               cacaGGCCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATcc                                             >  7:88549/1‑148 (MQ=255)
                                                acaGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCAGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGc                                           >  4:123632/1‑149 (MQ=255)
                                                aacggTCCCGAACGCTGTGCCGTCCGGGTGGGCGGGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATc                                              <  8:343797/143‑1 (MQ=255)
                                                       tcGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGt                                    >  5:133196/1‑149 (MQ=255)
                                                        cGAACGCTGTGACGTACCTGTGGGAGCGGGCGTGCCCACGAATGGCCGCAACCCTTTTTTCTGGTTGCCGCCAGCCCCGCCAACACAAAGCGCTGGTAAAGTTGTTCTTTAAGCCCCTGCGGCTGCCCCAAATGCATCCGCCCCAGGTg                                   >  2:148676/1‑149 (MQ=255)
                                                               tgtgCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTc                             >  7:307309/1‑148 (MQ=255)
                                                               tgtgCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTc                             >  3:344870/1‑148 (MQ=255)
                                                               tgtgCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCa                            >  5:51584/1‑149 (MQ=255)
                                                               ggtgCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCa                            <  6:327717/148‑1 (MQ=255)
                                                                  gCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATAc                         <  2:81860/149‑1 (MQ=255)
                                                                   ccGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATAc                         >  8:86434/1‑148 (MQ=255)
                                                                   ccGTACCTGTGGGAGCGGGAGTGCCCGAGAATGGCCGAAACCCTTTTATCTGGTTGCCGACAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCAATGCGGCTGCCACAACTGCATCCGCCCCAGGTGGATTTCATAc                         >  7:209471/1‑148 (MQ=255)
                                                                         cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTcc                   >  3:24474/1‑148 (MQ=255)
                                                                         cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTcc                   >  8:160398/1‑148 (MQ=255)
                                                                         cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTcc                   >  8:180934/1‑148 (MQ=255)
                                                                         cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTcc                   <  4:24474/148‑1 (MQ=255)
                                                                         cTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCg                  >  5:170833/1‑149 (MQ=255)
                                                                          tgtgGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCgc                 >  2:164976/1‑148 (MQ=255)
                                                                            tgGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGcc               >  2:130543/1‑149 (MQ=255)
                                                                             gggAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGccc              >  2:63750/1‑149 (MQ=255)
                                                                               gAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGcccc             >  7:361142/1‑148 (MQ=255)
                                                                                    ggCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCa       >  2:55721/1‑149 (MQ=255)
                                                                                         gCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCaa  >  4:31301/1‑149 (MQ=255)
                                                                                            |                                                                                                                                                 
TAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAA  >  NC_002947_CJ‑RC/6016608‑6016845

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ 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
CGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGG  >  NC_002947_CJ‑RC/6016614‑6016836
                                                                                      |                                                                                                                                        
ccaccggtgctgcggggccccccccgccccaacgggggacccccgggccccaccggtgtgcgtccaggtggggggggggggGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAG                                                                          <  NB501016:282:H2V7GAFX2:1:11108:3860:13509/70‑1 (MQ=60)
                    CGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCCACCCTTTTATCTGGTTTCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCG                                                      >  NB501016:282:H2V7GAFX2:4:21608:10690:11929/1‑151 (MQ=60)
                    CGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCGTTTATCTGGTTGCCGCCCGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTTAGCCCCTGCG                                                      >  NB501016:282:H2V7GAFX2:3:21404:23321:6699/1‑151 (MQ=60)
                     GCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGG                                                     >  NB501016:282:H2V7GAFX2:1:11106:10658:19238/1‑151 (MQ=60)
                              CACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAA                                                                                                                                     >  NB501016:282:H2V7GAFX2:2:21212:7705:2450/1‑62 (MQ=60)
                                     cccccaccggtgcccacccctggcccgaccgggggggcggGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATC                                     <  NB501016:282:H2V7GAFX2:1:11209:24396:17597/111‑1 (MQ=60)
                                        GCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGC                                  >  NB501016:282:H2V7GAFX2:1:21306:19126:4936/1‑151 (MQ=60)
                                        GCACAGGCCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGC                                  >  NB501016:282:H2V7GAFX2:4:11505:14743:13882/1‑151 (MQ=60)
                                             ggcccctatcggttgccgtaccgtggggggggGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCC                               <  NB501016:282:H2V7GAFX2:4:21612:8953:17119/117‑1 (MQ=60)
                                                CTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTG                          >  NB501016:282:H2V7GAFX2:3:11603:12864:1362/1‑151 (MQ=60)
                                                   caccccggggcgcagggtggggggggggggGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCT                       <  NB501016:282:H2V7GAFX2:3:11511:10547:19633/121‑1 (MQ=60)
                                                        CTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCAT                  >  NB501016:282:H2V7GAFX2:3:11411:15108:14487/1‑151 (MQ=60)
                                                        CTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCAT                  >  NB501016:282:H2V7GAFX2:4:21511:13833:15103/1‑151 (MQ=60)
                                                        CTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCA                   >  NB501016:282:H2V7GAFX2:2:21310:17264:8731/1‑150 (MQ=60)
                                                            GCCGTACCTGTGGGAGCGGGAGTGCCCGAGAATGGCCGAAACCCTTTTATCTGGTTGCCGACAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCAATGCGGCTGCCACAACTGCATCCGCCCCAGGTGGATTTCATACCC              >  NB501016:282:H2V7GAFX2:4:21405:23635:3619/1‑151 (MQ=60)
                                                                  CCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGG        >  NB501016:282:H2V7GAFX2:3:11609:22455:15598/1‑151 (MQ=60)
                                                                  CCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCG         >  NB501016:282:H2V7GAFX2:2:11105:11335:18018/1‑150 (MQ=60)
                                                                        GGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGG  >  NB501016:282:H2V7GAFX2:4:21610:15283:8002/1‑151 (MQ=60)
                                                                                      |                                                                                                                                        
CGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGG  >  NC_002947_CJ‑RC/6016614‑6016836

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: