Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,811,134 A→G H122R (CAC→CGC)  mprA → DNA‑binding transcriptional repressor MprA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,811,1340AG100.0% 220.1 / NA 61H122R (CAC→CGC) mprADNA‑binding transcriptional repressor MprA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/34);  total (27/34)

AGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACA  >  NC_000913/2811001‑2811269
                                                                                                                                     |                                                                                                                                       
aGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGtt                                                                                                                                   >  1:1249545/1‑140 (MQ=255)
      gCTCTTGGATCATCCCGTACAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTgcgc                                                                                                                            <  2:819859/140‑1 (MQ=255)
       ctctTGGATCATCCCGTACCAACGCAACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTgcgc                                                                                                                            >  1:1177002/1‑140 (MQ=255)
            ggATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGt                                                                                                                       >  1:854737/1‑140 (MQ=255)
                 atcCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTac                                                                                                                  <  1:988932/140‑1 (MQ=255)
                  tcCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTgcgc                                                                                                                            >  2:613754/1‑129 (MQ=255)
                  tcCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTgcgc                                                                                                                            <  1:613754/129‑1 (MQ=255)
                  tcCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTacc                                                                                                                 >  1:994853/1‑140 (MQ=255)
                    ccGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTaccac                                                                                                               >  2:421075/1‑140 (MQ=255)
                    ccGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTaccac                                                                                                               <  2:1275949/140‑1 (MQ=255)
                    ccGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTaccac                                                                                                               <  2:910306/140‑1 (MQ=255)
                      gTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCg                                                                                                             <  1:579587/140‑1 (MQ=255)
                      gTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCg                                                                                                             <  2:9503/140‑1 (MQ=255)
                      gTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTGTTTGCGCGAGGTTTTACCACCg                                                                                                             <  2:672717/140‑1 (MQ=255)
                       tACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTACCACCgca                                                                                                           <  1:37962/140‑1 (MQ=255)
                           aaCGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCgcagca                                                                                                        <  2:1318791/140‑1 (MQ=255)
                             cgcgtcgcgtATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATa                                                                                                      >  1:1158108/1‑140 (MQ=255)
                             cgcgacgcgTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATa                                                                                                      >  2:488499/1‑140 (MQ=255)
                             cgcgacgcgTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATa                                                                                                      >  2:101556/1‑140 (MQ=255)
                             cgcgacgcgTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATa                                                                                                      >  2:1029813/1‑140 (MQ=255)
                             cgcgacgcgTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATa                                                                                                      >  2:718309/1‑140 (MQ=255)
                                  cgcgTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAActgc                                                                                                 >  1:157139/1‑140 (MQ=255)
                                   gcgTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAActgcc                                                                                                <  1:718309/140‑1 (MQ=255)
                                       aTTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCa                                                                                            >  2:281037/1‑140 (MQ=255)
                                        ttGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCAt                                                                                           <  1:421075/140‑1 (MQ=255)
                                        ttGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCAt                                                                                           <  1:985614/140‑1 (MQ=255)
                                        ttGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCAt                                                                                           <  1:488499/140‑1 (MQ=255)
                                        ttGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCAt                                                                                           <  2:821791/140‑1 (MQ=255)
                                            cGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCAt                                                                                           <  1:835301/136‑1 (MQ=255)
                                            cGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCAt                                                                                           >  2:835301/1‑136 (MQ=255)
                                              aTGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAActc                                                                                     >  2:911922/1‑140 (MQ=255)
                                                   cTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTc                                                                                >  2:482390/1‑140 (MQ=255)
                                                     ggAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCcg                                                                              >  2:420995/1‑140 (MQ=255)
                                                       aaaaaCGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCcgcg                                                                            <  1:1216647/140‑1 (MQ=255)
                                                         aaaCGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCt                                                                          <  2:715801/140‑1 (MQ=255)
                                                            cgcgGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAg                                                                       <  2:1158108/140‑1 (MQ=255)
                                                             gcgGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGc                                                                      >  1:209293/1‑140 (MQ=255)
                                                                           cgtcgtGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGa                                                        >  2:532196/1‑140 (MQ=255)
                                                                               gtGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAg                                                    >  2:1076309/1‑140 (MQ=255)
                                                                                tGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGc                                                   <  1:1029813/140‑1 (MQ=255)
                                                                                tGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGc                                                   >  2:1098036/1‑140 (MQ=255)
                                                                                 gAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCt                                                  <  2:591133/140‑1 (MQ=255)
                                                                                  aaaGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACGACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTc                                                 >  2:1113760/1‑140 (MQ=255)
                                                                                  aaaGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACTGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTc                                                 <  2:1249545/140‑1 (MQ=255)
                                                                                  aaaGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTc                                                 >  2:873684/1‑140 (MQ=255)
                                                                                                    gcTGCCTGCATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCaaa                               <  2:157139/140‑1 (MQ=255)
                                                                                                       gcctgcATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTg                            <  2:1319/140‑1 (MQ=255)
                                                                                                       gcctgcATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTg                            <  1:482390/140‑1 (MQ=255)
                                                                                                       gcctgcATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTg                            <  1:420995/140‑1 (MQ=255)
                                                                                                       gcctgcATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTg                            >  2:1210707/1‑140 (MQ=255)
                                                                                                           gcATCTGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGctct                        <  2:693190/140‑1 (MQ=255)
                                                                                                                tGCAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGt                   <  1:873684/140‑1 (MQ=255)
                                                                                                                  cAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCtctc                       >  1:736540/1‑134 (MQ=255)
                                                                                                                  cAATTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCtctc                       <  2:736540/134‑1 (MQ=255)
                                                                                                                    aTTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGtct                 <  1:534546/138‑1 (MQ=255)
                                                                                                                    aTTAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGtct                 >  2:534546/1‑138 (MQ=255)
                                                                                                                      tAACGGAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGAc             <  1:532196/140‑1 (MQ=255)
                                                                                                                          ggAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGa         <  2:854737/140‑1 (MQ=255)
                                                                                                                          ggAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGa         <  2:209293/140‑1 (MQ=255)
                                                                                                                          ggAAAAAGGTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGa         >  2:325668/1‑140 (MQ=255)
                                                                                                                                 ggTCGCGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACa  <  1:911922/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                       
AGTTGTGCTCTTGGATCATCCCGTACCAACGCGACGCGTATTGCCGATGAACTGGAAAAACGCGGTTGGATCGAACGTCGTGAAAGCGATAACGATCGCCGCTGCCTGCATCTGCAATTAACGGAAAAAGGTCACGAGTTTTTGCGCGAGGTTTTACCACCGCAGCATAACTGCCTGCATCAACTCTGGTCCGCGCTCAGCACAACAGAAAAAGATCAGCTCGAGCAAATCACCCGCAAATTGCTCTCCCGTCTCGACCAGATGGAACA  >  NC_000913/2811001‑2811269

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