Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 838,766 A→G 6.7% F312S (TTC→TCC)  prs ← ribose‑phosphate pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947838,7660AG 6.7% 177.4 / 2.6 61F312S (TTC→TCC) prsribose‑phosphate pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base G (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

AATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCGTCCAGTACCGAC  >  NC_002947/838630‑838915
                                                                                                                                        |                                                                                                                                                     
aaTCCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGAt                                                                                                                                          <  1:1128956/150‑1 (MQ=255)
         aTTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCCTCGCGGAACATCGCGCTCATCGActcttc                                                                                                                                 >  2:582861/1‑150 (MQ=255)
              ggttCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACATCCGTTGGGGCGGGGCGCTTCACATCAACACGTCTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGc                                                                                                                            <  2:583764/148‑1 (MQ=255)
                   tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGa                                                                                                                                           >  2:818716/1‑130 (MQ=255)
                   tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGa                                                                                                                                           <  1:818716/130‑1 (MQ=255)
                   tcCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATg                                                                                                                       >  2:1645405/1‑150 (MQ=255)
                              cccccGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGa                                                                                                                         >  1:221801/1‑137 (MQ=255)
                              cccccGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACAGCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGa                                                                                                                         <  2:221801/137‑1 (MQ=255)
                               ccccGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCt                                                                                                           >  1:1497510/1‑150 (MQ=255)
                                    aGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGc                                                                                                      <  2:547253/150‑1 (MQ=255)
                                       gCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGAc                                                                                                   <  1:1153095/150‑1 (MQ=255)
                                         gTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGCTGCAGCTGACCGCTTCAGCGACaa                                                                                                 >  1:1318067/1‑150 (MQ=255)
                                          tttGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCGGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCATCGACAAc                                                                                                >  2:739526/1‑150 (MQ=255)
                                                         aaCGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAg                                                                                 >  1:891527/1‑150 (MQ=255)
                                                               gATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGTCg                                                                           <  1:1645405/150‑1 (MQ=255)
                                                                aTGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCgg                                                                          >  1:240398/1‑150 (MQ=255)
                                                                 tGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGa                                                                                        <  1:1122538/135‑1 (MQ=255)
                                                                 aggCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGa                                                                                        >  2:1122538/2‑135 (MQ=255)
                                                                     aaaaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGt                                                                                                                               <  1:772956/92‑1 (MQ=255)
                                                                     aaaaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGt                                                                                                                               >  2:772956/1‑92 (MQ=255)
                                                                      aaaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGc                                                                                          <  2:736196/128‑1 (MQ=255)
                                                                      aaaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGc                                                                                          >  1:736196/1‑128 (MQ=255)
                                                                      aaaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCCTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACgg                                                                    <  1:582861/150‑1 (MQ=255)
                                                                       aaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGaa                                                                                                                                                      <  1:1083310/67‑1 (MQ=255)
                                                                       aaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGaa                                                                                                                                                      >  2:1083310/1‑67 (MQ=255)
                                                                       aaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGc                                                                                          <  1:619135/127‑1 (MQ=255)
                                                                       aaaaGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGc                                                                                          >  2:619135/1‑127 (MQ=255)
                                                                            ccccgccccAACCTGTGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGa                                                                         <  2:1445148/139‑1 (MQ=255)
                                                                            ccccgccccAACCTGTGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGa                                                                         >  1:1445148/1‑139 (MQ=255)
                                                                             cccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAgg                                                             >  2:1356989/1‑150 (MQ=255)
                                                                             cccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAgg                                                             >  2:1554884/1‑150 (MQ=255)
                                                                              ccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTAAGCGGGACac                                                                                                                                                   >  2:621981/1‑62 (MQ=255)
                                                                              ccgccccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTAAGCGGGACac                                                                                                                                                   <  1:621981/63‑2 (MQ=255)
                                                                                   ccAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGag                                                       <  2:1308166/150‑1 (MQ=255)
                                                                                    cAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCg                                                                           <  1:577578/129‑1 (MQ=255)
                                                                                    cAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCg                                                                           >  2:577578/1‑129 (MQ=255)
                                                                                    cAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGagc                                                      >  2:72796/1‑150 (MQ=255)
                                                                                        acaGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCg                                                  >  1:20265/1‑150 (MQ=255)
                                                                                           ggTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCa                                               <  2:240398/150‑1 (MQ=255)
                                                                                            gTTGGGGCGGGGCGCTTCACATCAACACGTGTTCCGCTAAGCGGGACac                                                                                                                                                   <  2:646528/49‑2 (MQ=21)
                                                                                            gTTGGGGCGGGGCGCTTCACATCAACACGTGTTCCGCTAAGCGGGACac                                                                                                                                                   >  1:646528/1‑48 (MQ=21)
                                                                                            gTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGt                                                                   >  2:1283776/1‑129 (MQ=255)
                                                                                            gTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGt                                                                   <  1:1283776/129‑1 (MQ=255)
                                                                                             ttGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGCACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGa                                             <  1:589501/150‑1 (MQ=255)
                                                                                                  gcggggcgCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCg                                        <  2:20265/150‑1 (MQ=255)
                                                                                                  gcggggcgCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCg                                        >  2:1224502/1‑150 (MQ=255)
                                                                                                           ttCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGt                                                                   >  2:595396/1‑114 (MQ=255)
                                                                                                           ttCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGt                                                                   <  1:595396/114‑1 (MQ=255)
                                                                                                           ttCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGgtgt                               <  1:1421402/150‑1 (MQ=255)
                                                                                                                aTCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCt                                                           <  1:1520887/117‑1 (MQ=255)
                                                                                                                aTCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCt                                                           >  2:1520887/1‑117 (MQ=255)
                                                                                                                     cacGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAAc                                    <  1:118966/135‑1 (MQ=255)
                                                                                                                     cacGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAAc                                    >  2:118966/1‑135 (MQ=255)
                                                                                                                            cccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCg              >  2:434673/1‑150 (MQ=255)
                                                                                                                            cccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCg              >  2:1246477/1‑150 (MQ=255)
                                                                                                                             ccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGc                                                                                >  1:59874/1‑83 (MQ=255)
                                                                                                                             ccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGc                                                                                <  2:59874/83‑1 (MQ=255)
                                                                                                                             ccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGc                                         >  2:280277/1‑122 (MQ=255)
                                                                                                                             ccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGc                                         <  1:280277/122‑1 (MQ=255)
                                                                                                                             ccGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCGt             <  1:1246477/150‑1 (MQ=255)
                                                                                                                               gCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATa                                                                       <  2:1099806/90‑1 (MQ=255)
                                                                                                                               gCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATa                                                                       >  1:1099806/1‑90 (MQ=255)
                                                                                                                                cTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCGTCCa          <  2:1318067/150‑1 (MQ=255)
                                                                                                                                        aaCATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCGTCCAGTACCGAc  <  1:757872/150‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                     
AATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACACGTGTCCCGCTTAGCGGAACATCGCGCTGATCGATTCTTCGTTGCTGATGCGGCGTACCGCTTCAGCGACAACCGGTGCGATATCCAGCTGGCGGATACGGTCACAGGCTTGAGCAGCGGCAGACAGCGGAACGGTGTTGGTCACCACCAGCTCGTCCAGTACCGAC  >  NC_002947/838630‑838915

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