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
A1 F1 I1 R1 154 26.6 2069536 90.2% 1866721 108.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260480 1,166,676 +G intergenic (+150/‑109) h16_B1032 → / → h16_B1033 conserved hypothetical protein/probable extra‑cytoplasmic solute receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604801,166,6741.G100.0% 115.9 / NA 30intergenic (+148/‑111)h16_B1032/h16_B1033conserved hypothetical protein/probable extra‑cytoplasmic solute receptor
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (16/14);  total (16/14)

TGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGCGA  >  AM260480/1166571‑1166796
                                                                                                        |                                                                                                                          
tgtgATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGac                                                                                                >  1:852821/1‑133 (MQ=255)
          aaGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGc                                                                                     >  2:380962/1‑134 (MQ=255)
            ggCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGt                                                                                   >  2:508754/1‑134 (MQ=255)
                  aCCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCgg                                                                             <  1:200425/134‑1 (MQ=255)
                     cgcgACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGacaa                                                                                              >  2:920095/1‑114 (MQ=255)
                     cgcgACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGacaa                                                                                              <  1:920095/114‑1 (MQ=255)
                            gtgtGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGt                                                                   <  1:896925/134‑1 (MQ=255)
                                        ctctTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCg                                                       <  2:999045/134‑1 (MQ=255)
                                           tttCGCCTTGCTGCCGTGGACCCCCCAGAACTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCctc                                                                                                                 <  1:26980/73‑1 (MQ=255)
                                           tttCGCCTTGCTGCCGTGGACCCCCCAGAACTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCctc                                                                                                                 >  2:26980/1‑73 (MQ=255)
                                               gCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCga                                                 <  2:852821/133‑1 (MQ=255)
                                                       gcCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTgag                                         >  1:1004540/1‑133 (MQ=255)
                                                       gcCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTgag                                         >  2:778856/1‑133 (MQ=255)
                                                                 ccccAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCa                              <  2:800659/134‑1 (MQ=255)
                                                                  cccAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCaa                             <  1:994505/134‑1 (MQ=255)
                                                                   ccAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAAt                            <  1:508754/134‑1 (MQ=255)
                                                                      gACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATc                           >  1:445273/1‑132 (MQ=255)
                                                                      gACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATc                           >  1:3009/1‑132 (MQ=255)
                                                                      gACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATc                           >  2:590214/1‑132 (MQ=255)
                                                                         cTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAagag                      >  2:826312/1‑134 (MQ=255)
                                                                           cAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAagag                      >  2:801230/1‑132 (MQ=255)
                                                                           cAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAagag                      <  1:801230/132‑1 (MQ=255)
                                                                                gCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGa               >  2:553931/1‑134 (MQ=255)
                                                                                       ggCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGTGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAg         >  1:677666/1‑133 (MQ=255)
                                                                                       ggCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAg         >  1:911364/1‑133 (MQ=255)
                                                                                         ccGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGaa      <  2:3009/134‑1 (MQ=255)
                                                                                          cGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGTGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAg     <  2:677666/134‑1 (MQ=255)
                                                                                             tACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGCGa  <  1:590214/134‑1 (MQ=255)
                                                                                                  cgAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGt                    <  1:945332/111‑1 (MQ=255)
                                                                                                  cgAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGt                    >  2:945332/1‑111 (MQ=255)
                                                                                                        |                                                                                                                          
TGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGCGA  >  AM260480/1166571‑1166796

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 37 ≤ 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
GGAGGAACTGACGGCCCGGCATCCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGCGAACCAATAACCTATGGTGGCGCGGATACGC  >  AM260480/1166536‑1166825
                                                                                                                                             |                                                                                                                                                     
GGAGGAACTGACGGCCCGGCATCCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCC                                                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1104:16843:90059/149‑1 (MQ=60)
 CAGGAACTGACGGCCCGGCATCCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTC                                                                                                                                              <  GWNJ‑0478:712:GW2002102894th:2:1110:18583:76413/150‑1 (MQ=60)
                gagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGAACTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCG                                                                                                                               <  GWNJ‑0478:712:GW2002102894th:2:1103:3509:90502/104‑1 (MQ=60)
                    ATCCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACA                                                                                                                            >  GWNJ‑0478:712:GW2002102894th:2:2209:9713:38432/1‑149 (MQ=60)
                                   agacagGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATC                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:2211:19827:57453/144‑1 (MQ=60)
                                                    GACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACG                                                                                           <  GWNJ‑0478:712:GW2002102894th:2:1112:20603:80318/150‑1 (MQ=60)
                                                              GGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATC                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2208:15945:14343/150‑1 (MQ=60)
                                                                           CTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGA                                                                     >  GWNJ‑0478:712:GW2002102894th:2:2215:10369:34793/1‑149 (MQ=60)
                                                                                       catacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaggtGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATC                                                        <  GWNJ‑0478:712:GW2002102894th:2:2107:16286:82662/96‑1 (MQ=60)
                                                                                          GCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCA                                                       >  GWNJ‑0478:712:GW2002102894th:2:1213:9705:97263/1‑148 (MQ=60)
                                                                                          GCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCA                                                       >  GWNJ‑0478:712:GW2002102894th:2:1101:6822:24473/1‑148 (MQ=60)
                                                                                                    CCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGAC                                           <  GWNJ‑0478:712:GW2002102894th:2:2214:16915:81168/150‑1 (MQ=60)
                                                                                                     ggctcggagatgtgtataagagacaGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACA                                          <  GWNJ‑0478:712:GW2002102894th:2:2208:15742:82214/125‑1 (MQ=60)
                                                                                                     CCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACA                                          <  GWNJ‑0478:712:GW2002102894th:2:1216:19109:41088/150‑1 (MQ=60)
                                                                                                           CCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGTGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGG                                     >  GWNJ‑0478:712:GW2002102894th:2:2110:15361:50837/1‑149 (MQ=60)
                                                                                                           CCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGA                                    <  GWNJ‑0478:712:GW2002102894th:2:2203:20397:94937/150‑1 (MQ=60)
                                                                                                           CCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGG                                     >  GWNJ‑0478:712:GW2002102894th:2:2211:17421:4868/1‑149 (MQ=60)
                                                                                                             aagagacagGTGTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAG                                  <  GWNJ‑0478:712:GW2002102894th:2:2212:6314:6428/141‑1 (MQ=60)
                                                                                                                        GTGGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACActgtctcttatacacatct                       >  GWNJ‑0478:712:GW2002102894th:2:1104:13718:83067/1‑131 (MQ=60)
                                                                                                                          GGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACActgtctcttatacacatctcc                     >  GWNJ‑0478:712:GW2002102894th:2:1113:16937:64798/1‑129 (MQ=60)
                                                                                                                          GGCCGATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGCGAACCAATAACC                     >  GWNJ‑0478:712:GW2002102894th:2:2110:7261:55710/1‑150 (MQ=60)
                                                                                                                               ATACCGCGAGCAGGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGCGAACCAATAACCTATGG                <  GWNJ‑0478:712:GW2002102894th:2:2107:3394:5038/150‑1 (MQ=60)
                                                                                                                                             GTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCAGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGctgtctcttatacacatctccgagcccacgagacg  >  GWNJ‑0478:712:GW2002102894th:2:2207:2259:24712/1‑115 (MQ=60)
                                                                                                                                             |                                                                                                                                                     
GGAGGAACTGACGGCCCGGCATCCCGGACCCAACATGTGATGGCGAAGGCGGGACCCGCGACGGTGTGGGTCCCGCTCTTTCGCCTTGCTGCCGTGGACCCCCCAGACCTCAATTGCCGTGTGGCCGATACCGCGAGCAGGTTTTCCCTCTGTTTGCGCGAATCGACAACAACTTGGCGTGATCGGTCAAGCCGGTTGACGCGTGGCGATCCGAGAATTGAGAAAAGACGGCAATCAAGAGGTACAGACATGAGGAAGCGAACCAATAACCTATGGTGGCGCGGATACGC  >  AM260480/1166536‑1166825

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