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
A4 F18 I1 R1 54 68.8 3145800 98.6% 3101758 145.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 1,419,549 +C intergenic (‑38/+98) PP_1242 ← / ← PP_1244 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC1,419,5481.C100.0% 91.5 / NA 27intergenic (‑37/+99)PP_1242/PP_1244hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (7/20);  total (7/20)

CTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGG  >  NC_002947_CJ‑RC/1419401‑1419620
                                                                                                                                                    |                                                                        
cTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctggc                                                                         <  1:248259/149‑1 (MQ=255)
      ggggATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctggccact                                                                     >  7:279775/1‑143 (MQ=255)
          aTGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGCGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCgg                                                                <  2:289917/149‑1 (MQ=255)
          aTGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGCGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCgg                                                                <  2:295240/149‑1 (MQ=255)
                  ttCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGccc                                                        <  6:300365/149‑1 (MQ=255)
                      aTGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTc                                                    >  1:125121/1‑149 (MQ=255)
                                 cGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACt                                         <  1:318645/149‑1 (MQ=255)
                                        ctgttgtTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAg                                   <  3:243310/147‑1 (MQ=255)
                                         tgttgtTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACt                                         >  1:115923/1‑141 (MQ=255)
                                         tgttgtTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACt                                         <  2:115923/141‑1 (MQ=255)
                                         tgttgtTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTg                                 >  2:291697/1‑149 (MQ=255)
                                           ttgttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAg                                   >  4:243310/1‑145 (MQ=255)
                                              ttttCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCa                                    >  8:72333/1‑141 (MQ=255)
                                              ttttCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCa                                    <  7:72333/141‑1 (MQ=255)
                                                ttCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACg                          >  4:115763/1‑149 (MQ=255)
                                                  cGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACgg                         <  2:159338/148‑1 (MQ=255)
                                                     cTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgt                     <  3:115763/149‑1 (MQ=255)
                                                        ggTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgtg                    <  4:38684/147‑1 (MQ=255)
                                                         gTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtgtg                                      <  2:387897/128‑1 (MQ=255)
                                                         gTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtgtg                                      >  1:387897/1‑128 (MQ=255)
                                                          tAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGa                <  4:308558/149‑1 (MQ=255)
                                                            gAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACAGTGCAGTGCGTGACGGCTGTGGTGAc               <  2:91731/148‑1 (MQ=255)
                                                                 aGGTCTTGCCGACTGTGGGGGGTTCTGTTGTGTTGTCGGTGGTCTTTGTTGTGTTTCCCTGGATTCTGGGTGGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtgtg         <  2:125121/149‑1 (MQ=38)
                                                                  ggTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTgg        <  3:135609/149‑1 (MQ=255)
                                                                  ggTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTgg        <  3:141307/149‑1 (MQ=255)
                                                                  ggTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTgg        <  2:159528/149‑1 (MQ=255)
                                                                   gTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTggg       <  7:153964/149‑1 (MQ=255)
                                                                    tCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGa      <  1:291697/149‑1 (MQ=255)
                                                                         ccGACTGTGGTGGGATCTGTTGTGTTGTCGGTGCGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCGCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCgg  <  3:288779/148‑1 (MQ=38)
                                                                                                                                                    |                                                                        
CTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGG  >  NC_002947_CJ‑RC/1419401‑1419620

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAG  >  NC_002947_CJ‑RC/1419400‑1419663
                                                                                                                                                      |                                                                                                                  
CCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCC                                                                                                                   <  NB501016:282:H2V7GAFX2:1:21109:25564:14130/151‑1 (MQ=60)
      AGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCcactt                                                                                                               >  NB501016:282:H2V7GAFX2:4:21501:16597:18022/1‑144 (MQ=60)
                      GATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCC                                                                                              >  NB501016:282:H2V7GAFX2:1:11211:26768:1423/1‑151 (MQ=60)
                                 GCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTG                                                                                   <  NB501016:282:H2V7GAFX2:1:21210:13035:15679/151‑1 (MQ=60)
                                        tcTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGT                                                                             <  NB501016:282:H2V7GAFX2:2:21109:25982:7009/148‑1 (MQ=60)
                                         GTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTG                                                                                   >  NB501016:282:H2V7GAFX2:1:11209:9568:7858/1‑143 (MQ=60)
                                         GTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCctcttatggggggggcccccaaccaccacagatgtcttgt                                                                            >  NB501016:282:H2V7GAFX2:4:11409:6944:6103/1‑109 (MQ=60)
                                             TGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTACCACAGtacctgtgggaggggg                                                                         >  NB501016:282:H2V7GAFX2:4:11608:21436:3297/1‑133 (MQ=60)
                                              GTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAG                                                                              <  NB501016:282:H2V7GAFX2:4:11501:10090:6059/143‑1 (MQ=60)
                                                     CCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTG                                                               <  NB501016:282:H2V7GAFX2:2:11210:9292:4174/151‑1 (MQ=60)
                                                         GGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGC                                                                                >  NB501016:282:H2V7GAFX2:1:21312:22701:17378/1‑130 (MQ=60)
                                                               AGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCcgcctggcggggggggcccccccccccccccgggggggtgggggggggggggggggggggtg                                                      >  NB501016:282:H2V7GAFX2:3:21512:1517:16041/1‑87 (MQ=60)
                                                                  AGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGG                                                  <  NB501016:282:H2V7GAFX2:2:11301:1327:14409/151‑1 (MQ=60)
                                                                  AGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGG                                                  <  NB501016:282:H2V7GAFX2:2:11302:26770:14397/151‑1 (MQ=60)
                                                                   GGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGA                                                 <  NB501016:282:H2V7GAFX2:4:11603:18138:16442/151‑1 (MQ=60)
                                                                    GTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAG                                                <  NB501016:282:H2V7GAFX2:1:21205:14711:15176/151‑1 (MQ=60)
                                                                         GCCGACTGTGGTGGGATCTGTTGTGTTGTCGGTGCGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCGCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGG                                            <  NB501016:282:H2V7GAFX2:2:21207:9916:14566/150‑1 (MQ=60)
                                                                                                           GGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCcgcggggagggggggggcccccccccccccccccgaggtgtgtggggggggggggggggggggggggggggggggggggggggcgggccagggggcccgggcgg            >  NB501016:282:H2V7GAFX2:1:21212:15566:7854/1‑43 (MQ=60)
                                                                                                                 GTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCcgcgatgccgggggggcccccccccccccccccggggtggggggggggggggggggggggggcggggggggggggggggggggcccgacgcgggcgcggcggatagaaggggtt  >  NB501016:282:H2V7GAFX2:1:21202:21003:16928/1‑37 (MQ=60)
                                                                                                                 GTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCccgattgagggggggccccccccccccccccaggtgtggtggggggggggggggggggggggcgggggggggggggggggggcggcgcgcgggggcgggcgggggagtgattgt  >  NB501016:282:H2V7GAFX2:1:11201:8853:13829/1‑37 (MQ=60)
                                                                                                                                                      |                                                                                                                  
CCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAG  >  NC_002947_CJ‑RC/1419400‑1419663

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: