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 146 27.9 2327082 89.7% 2087392 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 1,841,200 +G intergenic (‑2736/‑2680) h16_A1684 ← / → ycaC conserved hypothetical protein/YcaC related amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604791,841,1971.G100.0% 110.1 / NA 29intergenic (‑2733/‑2683)h16_A1684/ycaCconserved hypothetical protein/YcaC related amidohydrolase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (13/16);  total (13/16)

GATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGGT  >  AM260479/1841067‑1841323
                                                                                                                                   |                                                                                                                              
gATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAggg                                                                                                                             >  1:1042724/1‑134 (MQ=255)
                            gcCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGaa                                                                                                  <  1:919188/134‑1 (MQ=255)
                              ccAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACt                                                                                                <  2:1042724/134‑1 (MQ=255)
                                      cGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGc                                                                                        <  2:998737/134‑1 (MQ=255)
                                                      gccggTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCg                                                                                              >  2:415973/1‑112 (MQ=255)
                                                      gccggTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCg                                                                                              <  1:415973/112‑1 (MQ=255)
                                                       ccggTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggt                                                                           <  1:994778/130‑1 (MQ=255)
                                                       ccggTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggt                                                                           >  2:994778/1‑130 (MQ=255)
                                                                         tCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggt                                                                           <  2:926852/112‑1 (MQ=255)
                                                                         tCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggt                                                                           >  1:926852/1‑112 (MQ=255)
                                                                         tCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTtctc                                                     >  1:133819/1‑134 (MQ=255)
                                                                             ggCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggt                                                                           <  2:913673/108‑1 (MQ=255)
                                                                             ggCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCtggt                                                                           >  1:913673/1‑108 (MQ=255)
                                                                              gCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATc                                                            <  2:910408/122‑1 (MQ=255)
                                                                              gCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATc                                                            >  1:910408/1‑122 (MQ=255)
                                                                                    gTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGcc                                          >  2:107968/1‑134 (MQ=255)
                                                                                     tAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGccc                                         >  1:109090/1‑134 (MQ=255)
                                                                                             cgccTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGa                                 <  1:302043/134‑1 (MQ=255)
                                                                                              gccTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAg                                <  1:713576/134‑1 (MQ=255)
                                                                                                cTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCa                                                                                       <  2:587687/77‑1 (MQ=255)
                                                                                                cTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCa                                                                                       >  1:587687/1‑77 (MQ=255)
                                                                                                 tGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTc                                                                       <  2:308879/92‑1 (MQ=255)
                                                                                                 tGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTc                                                                       >  1:308879/1‑92 (MQ=255)
                                                                                                   cAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAg                            <  1:951558/133‑1 (MQ=255)
                                                                                                   cAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGc                           >  1:46094/1‑134 (MQ=255)
                                                                                                                tAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGc              <  2:133819/134‑1 (MQ=255)
                                                                                                                  gTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCt                                                   <  1:796916/95‑1 (MQ=255)
                                                                                                                  gTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCt                                                   >  2:796916/1‑95 (MQ=255)
                                                                                                                          ggAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCgtg    >  1:188981/1‑134 (MQ=255)
                                                                                                                            aaGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCgtggt  <  2:109090/134‑1 (MQ=255)
                                                                                                                                   |                                                                                                                              
GATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGGT  >  AM260479/1841067‑1841323

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTTCTGGGCCTTGGCGATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGGTGGCGCGCCTGCGCTCGAGTAGTGGC  >  AM260479/1841052‑1841348
                                                                                                                                                     |                                                                                                                                                    
CTTCTGGGCCTTGGCGATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGG                                                                                                                                                     >  GWNJ‑0478:712:GW2002102894th:2:2210:9073:8707/1‑150 (MQ=60)
  TCTGGGCCTTGGCGATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGG                                                                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:2205:20152:18533/148‑1 (MQ=60)
                                          CGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGC                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:2204:15252:10318/150‑1 (MQ=60)
                                              agagacagGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGC                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1208:12530:68530/142‑1 (MQ=60)
                                                                  CCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATC                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:2207:5432:48732/149‑1 (MQ=60)
                                                                         GTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCtgtctct                                                                             >  GWNJ‑0478:712:GW2002102894th:2:2204:14951:56209/1‑143 (MQ=60)
                                                                         GTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCC                                                                             >  GWNJ‑0478:712:GW2002102894th:2:1109:21301:4462/1‑150 (MQ=60)
                                                                             GGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCtgtctcttat                                                                          >  GWNJ‑0478:712:GW2002102894th:2:2202:17238:76386/1‑139 (MQ=60)
                                                                              GCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCG                                                                        >  GWNJ‑0478:712:GW2002102894th:2:2202:8871:56497/1‑150 (MQ=60)
                                                                                     GCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCG                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1106:8133:64746/1‑150 (MQ=60)
                                                                                         gctcggagatgtgtataagagacaGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCA                                                             <  GWNJ‑0478:712:GW2002102894th:2:2112:19938:36862/126‑1 (MQ=60)
                                                                                                GCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCtgtctcttatacacatctccgagcccacgagacactcgcta                                                       >  GWNJ‑0478:712:GW2002102894th:2:2101:18794:31217/1‑108 (MQ=60)
                                                                                                 CTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTgtctcttatacacatctccgagccc                                                      >  GWNJ‑0478:712:GW2002102894th:2:1203:18768:40016/1‑124 (MQ=60)
                                                                                                   GTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCC                                                   >  GWNJ‑0478:712:GW2002102894th:2:1105:15347:85528/1‑150 (MQ=60)
                                                                                                           CCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGAT                                           <  GWNJ‑0478:712:GW2002102894th:2:1201:4863:96850/150‑1 (MQ=60)
                                                                                                            CGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATG                                          <  GWNJ‑0478:712:GW2002102894th:2:2108:4821:94486/150‑1 (MQ=60)
                                                                                                             gagtgtctcgtgggctcggagatgtgtataagagaCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGG                                        <  GWNJ‑0478:712:GW2002102894th:2:2102:1455:20822/115‑1 (MQ=60)
                                                                                                                 GCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCC                                      <  GWNJ‑0478:712:GW2002102894th:2:2209:16923:3615/149‑1 (MQ=60)
                                                                                                                          CATGATAGTCGCGCAGGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGG                            >  GWNJ‑0478:712:GW2002102894th:2:1113:16773:2101/1‑150 (MQ=60)
                                                                                                                                         GGAAGTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGGTGGCGCGCCTGCGCT             >  GWNJ‑0478:712:GW2002102894th:2:1210:2921:67109/1‑150 (MQ=60)
                                                                                                                                             GTTCAGGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCtgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtcttctgc         >  GWNJ‑0478:712:GW2002102894th:2:1206:13140:30159/1‑88 (MQ=60)
                                                                                                                                                    GGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCtgtctcttatacacatctccgagcccacgagacactc  >  GWNJ‑0478:712:GW2002102894th:2:1211:2368:48389/1‑113 (MQ=60)
                                                                                                                                                     |                                                                                                                                                    
CTTCTGGGCCTTGGCGATCTCGAAGAAGCGTCGGCGGCAATGCGCCCAGCATCCGCCCTCGATGATCCGGCCGGTACGGCAAAGCGCATCGTGGCCGCTGTAGGCAGCCGCCTGCAGGTAGCCATGATAGTCGCGCAGGAAGTTCAGGGATGGGCGCCGGCGCGGGATTCCGCGAACTCGAACAGCACCGCTGGCTGGTGGTCGACCCATAATCCGTTCTCCTGCCGCAGCCCGGCCACGAGGTAGCCCCATAGATGGGCCTGCCGCGTGGTGGCGCGCCTGCGCTCGAGTAGTGGC  >  AM260479/1841052‑1841348

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