Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,853,037 A→G 7.0% S688S (AGT→AGC dnaX ← DNA polymerase III subunit gamma

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,853,0370AG 7.0% 161.0 / 2.9 57S688S (AGT→AGCdnaXDNA polymerase III subunit gamma
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base G (2/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

CTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTGCCGCCGCCTGAGCCG  >  NC_002947/4852890‑4853185
                                                                                                                                                    |                                                                                                                                                    
cTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGAc                                                                                                                                                     >  8:9831/1‑149 (MQ=255)
     tcttctGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGc                                                                                                                                                 >  6:398681/1‑148 (MQ=255)
       ttctGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCTGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAg                                                                                                                                              >  1:69890/1‑149 (MQ=255)
          tGCATTTTTTCCTGCATCTGCTGGGACTGCTTCATCAGGCCGGCAATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGc                                                                                                                                           >  7:129786/1‑149 (MQ=255)
            cATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAggg                                                                                                                                            <  6:74486/146‑1 (MQ=255)
            cATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGc                                                                                                                                           <  1:237529/147‑1 (MQ=255)
            cATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTAAACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCtt                                                                                                                                         <  6:37321/149‑1 (MQ=255)
            cATTTTTTCCTGCATCTGCTGGGCCTGCTACATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAggg                                                                                                                                            <  6:132808/146‑1 (MQ=255)
             aTTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTc                                                                                                                                        >  3:37661/1‑149 (MQ=255)
              ttttttCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTc                                                                                                                                        >  5:163023/1‑148 (MQ=255)
                  ttCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGGCTAGTCAGGGccagtaat                                                                                                                                    <  7:50016/148‑8 (MQ=255)
                   tCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAg                                                                                                                                   <  5:339/148‑1 (MQ=255)
                      tGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGGCTAGTCAGGGccag                                                                                                                                        >  6:303453/1‑137 (MQ=255)
                      tGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGGCTAGTCAGGGccag                                                                                                                                        <  5:303453/140‑4 (MQ=255)
                      tGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCt                                                                                                                                >  7:150368/1‑148 (MQ=255)
                       gCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTc                                                                                                                               >  8:50120/1‑148 (MQ=255)
                        cATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCaa                                                                                                                             <  1:175690/149‑1 (MQ=255)
                            tgctgGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTGATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGt                                                                                                                         <  2:37717/149‑1 (MQ=255)
                              ctgGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATa                                                                                                                           <  6:366807/145‑1 (MQ=255)
                                       cTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCa                                                                                                              >  6:10850/1‑149 (MQ=255)
                                              aGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCTGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGgc                                                                                                       <  2:69890/149‑1 (MQ=255)
                                               ggcgggccATTCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGgc                                                                                                       <  8:129786/148‑1 (MQ=255)
                                                     ccATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGaac                                                                                                <  6:163023/149‑1 (MQ=255)
                                                          ccACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACt                                                                                            <  5:210117/148‑1 (MQ=255)
                                                            aCCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGa                                                                                         <  5:133087/149‑1 (MQ=255)
                                                             ccTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGat                                                                                        <  3:305829/149‑1 (MQ=255)
                                                              cTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGat                                                                                        >  4:360887/1‑148 (MQ=255)
                                                                   atcatGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATct                                                                                   <  8:150368/148‑1 (MQ=255)
                                                                      atGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATctgct                                                                                <  7:50120/148‑1 (MQ=255)
                                                                              taCCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAAc                                                                        <  7:310110/148‑1 (MQ=255)
                                                                                 cTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATctgc                                                                                 <  3:285447/136‑1 (MQ=255)
                                                                                 cTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATctgc                                                                                 >  4:285447/1‑136 (MQ=255)
                                                                                 cTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGAt                                                                    >  5:362933/1‑149 (MQ=255)
                                                                                 cTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGAt                                                                    >  1:68798/1‑149 (MQ=255)
                                                                                 cTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGAt                                                                    >  5:75760/1‑149 (MQ=255)
                                                                                    gATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCt                                                                 >  6:134378/1‑149 (MQ=255)
                                                                                         ggTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCtc                                                             >  1:405694/1‑148 (MQ=255)
                                                                                         ggTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACATTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCtc                                                             >  8:161426/1‑148 (MQ=255)
                                                                                           tCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCtcgat                                                          <  5:119716/149‑1 (MQ=255)
                                                                                           tCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCtcgat                                                          >  6:119716/1‑149 (MQ=255)
                                                                                            cATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCtcgatc                                                         >  7:224805/1‑149 (MQ=255)
                                                                                              tgtgATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCtcgatcga                                                       >  2:264491/1‑149 (MQ=255)
                                                                                                  aTGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATcgacgc                                                    <  6:75760/148‑1 (MQ=255)
                                                                                                  aTGCGGCCGGCCCGGTGCCGGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATcgacgc                                                    >  2:309964/1‑148 (MQ=255)
                                                                                                       gCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGc                                              >  5:345726/1‑149 (MQ=255)
                                                                                                          cGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGt                                            >  4:145624/1‑148 (MQ=255)
                                                                                                          cGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTc                                           <  1:264491/149‑1 (MQ=255)
                                                                                                          cGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTc                                           <  7:9831/149‑1 (MQ=255)
                                                                                                                       ggccgCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCgctc                               >  4:240447/1‑149 (MQ=255)
                                                                                                                        gccgCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCgctcg                              <  1:355719/149‑1 (MQ=255)
                                                                                                                        gccgCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCgctcg                              <  6:345726/149‑1 (MQ=255)
                                                                                                                            cATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCgctcggctc                          <  3:240447/149‑1 (MQ=255)
                                                                                                                                      cATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGATCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTgc                >  3:199018/1‑149 (MQ=255)
                                                                                                                                       aTTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTgcc               >  6:109194/1‑149 (MQ=255)
                                                                                                                                        gtagaagCCCTGGCTAGTCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTgccg              >  2:192133/7‑149 (MQ=255)
                                                                                                                                           cTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTgccgcc            <  2:405694/148‑1 (MQ=255)
                                                                                                                                                ccTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTGCCGCCGCCt        >  1:113119/1‑147 (MQ=255)
                                                                                                                                                    aCTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATAGAAGCGACGGAGTCCTGCTGGCGCTCGGCTCGCTTGCGTGCCGCCGCCTGAGCCg  >  6:51290/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGCATCTTCTGCATTTTTTCCTGCATCTGCTGGGCCTGCTTCATCAGGCCGGCCATGCCACCTTTCATCATGGGGTATACCTCGATTGGGTCATGTGATGCGGCCCGGCACCGGGCCGGCCGCATGGTTATCCATTACTGGCCCTGACTAGCCAGGGCTTCTACAGGCTCAATAGTATCCTGCCGCACCTTGGCGCCGAACAACTTGATCATCTGCTGGATCAACGGATCCTGCTCGATCGACGCGACGGCGTCCTGCTGGCGCTCGGCTCGCTTGCGTGCCGCCGCCTGAGCCG  >  NC_002947/4852890‑4853185

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