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 149 26.7 2337422 87.5% 2045244 105.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,302,920 +C coding (1199/1251 nt) dfp → phosphopantothenoylcysteine synthetase/decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,302,9181.C100.0% 131.8 / NA 34R399R (CGG→CGCdfpphosphopantothenoylcysteine synthetase/decarboxylase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (16/18);  total (16/18)

CTGGACGACAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATGC  >  AM260479/3302793‑3303041
                                                                                                                              |                                                                                                                           
cTGGACGACAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATg                                                                                                                      >  2:914588/1‑134 (MQ=255)
  ggACGACAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATggc                                                                                                                    >  1:852188/1‑134 (MQ=255)
   gacgacAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATggcg                                                                                                                   <  1:285469/134‑1 (MQ=255)
                  gTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTa                                                                                                    >  1:96251/1‑134 (MQ=255)
                       gTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATg                                                                                               >  2:501168/1‑134 (MQ=255)
                             cgccgcCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGc                                                                                         <  1:380252/134‑1 (MQ=255)
                                                            gACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTccccc                                                          <  2:96251/134‑1 (MQ=255)
                                                            gACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTccccc                                                          <  1:671778/134‑1 (MQ=255)
                                                             aCAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCa                                                         <  2:215318/134‑1 (MQ=255)
                                                                        cTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATg                                              <  1:914588/134‑1 (MQ=255)
                                                                        cTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATg                                              <  1:501168/134‑1 (MQ=255)
                                                                           gcgcgcCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGAcc                                           >  2:845498/1‑134 (MQ=255)
                                                                               gcCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGAt                                                                 >  2:1140056/1‑108 (MQ=255)
                                                                               gcCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGAt                                                                 <  1:1140056/108‑1 (MQ=255)
                                                                               gcCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGc                                       >  2:711771/1‑134 (MQ=255)
                                                                                ccAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGgcgc                                                                                                              <  2:1092546/62‑1 (MQ=255)
                                                                                ccAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGgcgc                                                                                                              >  1:1092546/1‑62 (MQ=255)
                                                                                                                 cGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGa     <  1:711771/134‑1 (MQ=255)
                                                                                                                   ggCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAAcc                              <  1:548736/107‑1 (MQ=255)
                                                                                                                   ggCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAAcc                              >  2:548736/1‑107 (MQ=255)
                                                                                                                   ggCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAAc                               <  2:735566/106‑1 (MQ=255)
                                                                                                                   ggCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAAc                               >  1:735566/1‑106 (MQ=255)
                                                                                                                   ggCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCaa                           >  2:581017/1‑110 (MQ=255)
                                                                                                                   ggCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCaa                           <  1:581017/110‑1 (MQ=255)
                                                                                                                   ggCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTc       >  2:975626/1‑130 (MQ=255)
                                                                                                                   ggCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTc       <  1:975626/130‑1 (MQ=255)
                                                                                                                   ggCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATg   >  1:631415/1‑134 (MQ=255)
                                                                                                                    gcgcACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATgc  <  1:845498/134‑1 (MQ=255)
                                                                                                                       cacaCGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCAcc                                                                                  <  2:382623/51‑1 (MQ=39)
                                                                                                                       cacaCGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCAcc                                                                                  <  2:574234/51‑1 (MQ=39)
                                                                                                                       cacaCGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCAcc                                                                                  >  1:574234/1‑51 (MQ=39)
                                                                                                                       cacaCGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCAcc                                                                                  <  2:1060639/51‑1 (MQ=39)
                                                                                                                       cacaCGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCAcc                                                                                  >  1:382623/1‑51 (MQ=39)
                                                                                                                       cacaCGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCAcc                                                                                  >  1:1060639/1‑51 (MQ=39)
                                                                                                                              |                                                                                                                           
CTGGACGACAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATGC  >  AM260479/3302793‑3303041

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CCATACCTTCGGCCTGGACGACAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATGCGCGCCTGAATGAATGGGG  >  AM260479/3302780‑3303059
                                                                                                                                             |                                                                                                                                           
CCATACCTTCGGCCTGGACGACAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCG                                                                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:2116:8488:19990/1‑150 (MQ=60)
               GGACGACAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTA                                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1201:10434:9580/150‑1 (MQ=60)
                GACGACAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTAT                                                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:1106:13595:2361/1‑150 (MQ=60)
                                         ACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGT                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:1207:9781:100586/150‑1 (MQ=60)
                                          ggttgtgggtggattctgctttttttttttcaagcagaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacaggcCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGT                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:2211:3194:18296/51‑1 (MQ=60)
                                                                  ataagagacagAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCC                                                                    <  GWNJ‑0478:712:GW2002102894th:2:2215:4284:76269/138‑1 (MQ=60)
                                                                        CGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACC                                                             <  GWNJ‑0478:712:GW2002102894th:2:2107:14840:92720/150‑1 (MQ=60)
                                                                              GCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCtgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtct                                                       >  GWNJ‑0478:712:GW2002102894th:2:2213:5792:59504/1‑93 (MQ=60)
                                                                                    ACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAAC                                                 <  GWNJ‑0478:712:GW2002102894th:2:1212:14130:54681/150‑1 (MQ=60)
                                                                                    ACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAAC                                                 <  GWNJ‑0478:712:GW2002102894th:2:2202:20460:82517/150‑1 (MQ=60)
                                                                                                   cgagtgtctcgtgggctcggagatgtgtataagagacaggcCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATC                                  <  GWNJ‑0478:712:GW2002102894th:2:1105:19199:99571/109‑1 (MQ=60)
                                                                                                    gtataagagacaGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCA                                 <  GWNJ‑0478:712:GW2002102894th:2:1214:2487:78571/138‑1 (MQ=60)
                                                                                                        aagagacaGGCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGG                              <  GWNJ‑0478:712:GW2002102894th:2:2101:3693:42400/141‑1 (MQ=60)
                                                                                                                 GCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCctgtctcttatac                     >  GWNJ‑0478:712:GW2002102894th:2:2111:10531:48069/1‑136 (MQ=60)
                                                                                                                 GCCAGCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATGC                    >  GWNJ‑0478:712:GW2002102894th:2:2105:11136:57514/1‑150 (MQ=60)
                                                                                                                     GCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaa                 >  GWNJ‑0478:712:GW2002102894th:2:1216:15118:98806/1‑81 (MQ=60)
                                                                                                                     GCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaa                 >  GWNJ‑0478:712:GW2002102894th:2:2212:7948:2680/1‑81 (MQ=60)
                                                                                                                     GCGCCTGCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaa                 >  GWNJ‑0478:712:GW2002102894th:2:1206:9563:15006/1‑81 (MQ=60)
                                                                                                                           GCCCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATGCGCGCCTGAAT          <  GWNJ‑0478:712:GW2002102894th:2:2206:10911:37494/150‑1 (MQ=60)
                                                                                                                             CCGGGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATGCGCGCCTGAATGA        <  GWNJ‑0478:712:GW2002102894th:2:2109:6466:22165/150‑1 (MQ=60)
                                                                                                                                GGCGCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATGCGCGCCTGAATGAATG     <  GWNJ‑0478:712:GW2002102894th:2:2114:5695:89019/150‑1 (MQ=60)
                                                                                                                                   GCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaataaaaacatattgcccatgaacactct  >  GWNJ‑0478:712:GW2002102894th:2:2214:9979:67420/1‑40 (MQ=60)
                                                                                                                                   GCACACGGCCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATGCGCGCCTGAATGAATGGGG  >  GWNJ‑0478:712:GW2002102894th:2:2114:18784:12486/1‑150 (MQ=60)
                                                                                                                                             |                                                                                                                                           
CCATACCTTCGGCCTGGACGACAACGAAATCGTGCTGTTCGACGCCGCCGGCATGACGCGGCTGCCGCGCGCCGACAAGCTGTCACTGGCGCGCCAGCTGGTGGCGGCGATCGGCCAGCGCCTGCCCGGGCGCACACGGCCTGATGGCGGGCGCGCGTCCGGTATCATGCGCGGCTGCCACCTATTCCGTTGAACCGATTTCCCCCAGGATTCCCATGACCCAGCCACTGAACCCAACCGTCGAAATCAAGGTGCTCGATGCGCGCCTGAATGAATGGGG  >  AM260479/3302780‑3303059

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