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 F2 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,211,170 (T)8→7 coding (1246/2187 nt) espL4 ← type III secretion system effector EspL4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,211,1630T.100.0% 204.0 / NA 47coding (1253/2187 nt)espL4type III secretion system effector EspL4
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (18/29);  total (18/29)

TACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTAT  >  NZ_CP009273/4211078‑4211248
                                                                                     |                                                                                     
tACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAttttt                                                                                   >  2:275782/1‑90 (MQ=255)
  caacaTTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAttttttt                                                                                 >  2:1185021/1‑90 (MQ=255)
  caacaTTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAttttttt                                                                                 >  2:1104186/1‑90 (MQ=255)
    acaTTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAtttttttcg                                                                               <  1:1046208/90‑3 (MQ=255)
    acaTTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAtttttttcg                                                                               <  1:1016311/90‑3 (MQ=255)
      aTTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGAtttttttcgag                                                                             >  1:733807/1‑86 (MQ=255)
        ttttATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGagag                                                                          >  2:736362/1‑90 (MQ=255)
           tATCCTTATGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTaa                                                                       <  2:1168089/90‑1 (MQ=255)
           tATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTaa                                                                       <  1:736362/90‑1 (MQ=255)
            aTCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTaa                                                                       <  2:533775/89‑1 (MQ=255)
            aTCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTaa                                                                       >  1:533775/1‑89 (MQ=255)
            aTCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAAt                                                                      >  1:64120/1‑90 (MQ=255)
             tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATc                                                                     <  2:412096/90‑1 (MQ=255)
             tCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATc                                                                     <  2:257377/90‑1 (MQ=255)
                ttACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCtt                                                                  >  1:152474/1‑90 (MQ=255)
                 tACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTc                                                                 >  2:324085/1‑90 (MQ=255)
                 tACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTc                                                                 >  2:718219/1‑90 (MQ=255)
                  aCGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCa                                                                <  1:214946/90‑1 (MQ=255)
                       atatCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATcc                                                           <  1:460023/90‑1 (MQ=255)
                       atatCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATcc                                                           <  1:1088923/90‑1 (MQ=255)
                       atatCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATcc                                                           <  2:12728/90‑1 (MQ=255)
                          tCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGt                                                        >  1:555383/1‑90 (MQ=255)
                          tCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGt                                                        >  1:1227369/1‑90 (MQ=255)
                           cGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGtt                                                       <  1:244438/90‑1 (MQ=255)
                            gCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGttt                                                      >  1:803325/1‑90 (MQ=255)
                            gCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGttt                                                      >  2:536471/1‑90 (MQ=255)
                                   ataatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACaa                                               <  2:852505/90‑1 (MQ=255)
                                   ataatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACaa                                               <  2:608048/90‑1 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                             <  1:929874/90‑1 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                             >  1:183606/1‑90 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                             >  2:57336/1‑90 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                             <  2:589763/90‑1 (MQ=255)
                                     aatCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCACGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATg                                             <  1:598927/90‑1 (MQ=255)
                                         cAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAAtt                                         <  1:15626/90‑1 (MQ=255)
                                           gAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGa                                       >  2:613386/1‑90 (MQ=255)
                                              aaaCCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGaaaa                                    <  2:1003799/90‑1 (MQ=255)
                                               aaCCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAAt                                   <  1:973233/90‑1 (MQ=255)
                                                        tATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAAc                          <  1:863253/90‑1 (MQ=255)
                                                         aTCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACg                         >  1:49882/1‑90 (MQ=255)
                                                           cTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGat                       >  1:547404/1‑90 (MQ=255)
                                                            ttttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGata                      <  2:350409/90‑1 (MQ=255)
                                                            ttttGCCCCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGata                      >  1:811064/1‑90 (MQ=255)
                                                              ttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATc                    >  2:65415/1‑90 (MQ=255)
                                                              ttGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATc                    <  2:74234/90‑1 (MQ=255)
                                                               tGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATcc                   <  1:324085/90‑1 (MQ=255)
                                                                  cTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCa                <  2:806140/90‑1 (MQ=255)
                                                                   tCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCAt               <  1:277482/90‑1 (MQ=255)
                                                                    ccAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCAtt              <  2:365833/90‑1 (MQ=255)
                                                                      aGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTAc            <  2:938217/90‑1 (MQ=255)
                                                                          tatGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATa        <  1:1109363/90‑1 (MQ=255)
                                                                          tatGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATa        <  2:90510/90‑1 (MQ=255)
                                                                           atGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTATGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATaa       >  2:563954/1‑90 (MQ=255)
                                                                                tGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTtat  <  1:1104186/90‑1 (MQ=255)
                                                                                     |                                                                                     
TACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTAT  >  NZ_CP009273/4211078‑4211248

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGGTTACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACC  >  NZ_CP009273/4211074‑4211263
                                                                                                |                                                                                             
AGGTTACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAG                                                                                           >  SRR3722089.743374/1‑100 (MQ=60)
      CAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAAT                                                                                     >  SRR3722089.64902/1‑100 (MQ=60)
      CAACATTATTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAAT                                                                                     >  SRR3722089.539993/1‑100 (MQ=60)
        ACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCC                                                                                   <  SRR3722089.1061890/100‑1 (MQ=60)
        ACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCC                                                                                   <  SRR3722089.1031454/100‑1 (MQ=60)
          ATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTT                                                                                 >  SRR3722089.154371/1‑100 (MQ=60)
               TATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATAT                                                                            <  SRR3722089.745987/100‑1 (MQ=60)
                    TTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGT                                                                       >  SRR3722089.1246017/1‑100 (MQ=60)
                    TTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGT                                                                       >  SRR3722089.561831/1‑100 (MQ=60)
                      ACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTT                                                                     >  SRR3722089.814480/1‑100 (MQ=60)
                      ACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTT                                                                     <  SRR3722089.217665/100‑1 (MQ=60)
                           ATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGAC                                                                <  SRR3722089.465518/100‑1 (MQ=60)
                           ATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGAC                                                                <  SRR3722089.1105347/100‑1 (MQ=60)
                               CGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATG                                                            <  SRR3722089.247551/100‑1 (MQ=60)
                               CGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATG                                                            >  SRR3722089.185899/1‑100 (MQ=60)
                                         AATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAAT                                                  <  SRR3722089.943453/100‑1 (MQ=60)
                                         AATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCACGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAAT                                                  <  SRR3722089.606277/100‑1 (MQ=60)
                                             CAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTC                                              <  SRR3722089.15822/100‑1 (MQ=60)
                                                   AACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATGGTCTCAACG                                        <  SRR3722089.987609/100‑1 (MQ=60)
                                                   AACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACG                                        >  SRR3722089.50502/1‑100 (MQ=60)
                                                     CCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGAT                                      >  SRR3722089.553769/1‑100 (MQ=60)
                                                      CATTTTTATCTTTTGCCCCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATA                                     >  SRR3722089.822396/1‑100 (MQ=60)
                                                            TATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCA                               <  SRR3722089.875604/100‑1 (MQ=60)
                                                                   TGCCTCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCAT                        <  SRR3722089.328161/100‑1 (MQ=60)
                                                                       TCCAGAATATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATT                    <  SRR3722089.281027/100‑1 (MQ=60)
                                                                              TATGCATGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCC             <  SRR3722089.1126149/100‑1 (MQ=60)
                                                                                    TGAGATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATAT       <  SRR3722089.1120882/100‑1 (MQ=60)
                                                                                         ATTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAAC   <  SRR3722089.310729/100‑1 (MQ=60)
                                                                                          TTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACC  <  SRR3722089.694303/100‑1 (MQ=60)
                                                                                                |                                                                                             
AGGTTACAACATTTTTATCCTTACGTGATATCGCTAAAAATAATCCAGAAAAACCATTTTTATCTTTTGCCTCCAGAATATGCATGAGATTTTTTTTCGAGAGTAATCCTTCATATCCTGTTTCTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACC  >  NZ_CP009273/4211074‑4211263

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