Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,747,020 G→A A259V (GCG→GTG)  ffh ← Signal Recognition Particle (SRP) component with 4.5S RNA (ffs)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,747,0200GA96.5% 184.0 / NA 57A259V (GCG→GTG) ffhSignal Recognition Particle (SRP) component with 4.5S RNA (ffs)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/27);  minor base C (0/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCGCCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAACCGGGTTAAT  >  NC_000913/2746882‑2747162
                                                                                                                                          |                                                                                                                                              
tCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCAccgcc                                                                                                                                      <  2:426408/149‑1 (MQ=255)
  gATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCgg                                                                                                                                  >  2:250339/1‑151 (MQ=255)
  gATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCgg                                                                                                                                  >  1:426618/1‑151 (MQ=255)
    tCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGc                                                                                                                                >  2:10423/1‑151 (MQ=255)
    tCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGc                                                                                                                                >  2:45065/1‑151 (MQ=255)
      aGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCAt                                                                                                                              <  2:178817/151‑1 (MQ=255)
        cGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCg                                                                                                                            <  1:185898/151‑1 (MQ=255)
         gACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCg                                                                                                                            >  2:657713/1‑150 (MQ=255)
         gACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGc                                                                                                                           >  1:102918/1‑151 (MQ=255)
          aCAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGcc                                                                                                                          >  1:850656/1‑151 (MQ=255)
              tACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTc                                                                                                                       <  2:162461/150‑1 (MQ=255)
                     cccATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTgg                                                                                                               >  2:250523/1‑151 (MQ=255)
                          gCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCaa                                                                                                           <  1:221496/150‑1 (MQ=255)
                           ccGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAac                                                                                                          >  2:716769/1‑150 (MQ=255)
                                     gcgacgcgaTGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGt                                                                                               <  1:982533/151‑1 (MQ=255)
                                     gcgacgcgaTGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGt                                                                                               <  1:349284/151‑1 (MQ=255)
                                       gacgcgaTGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTaa                                                                                             <  1:406128/151‑1 (MQ=255)
                                        acgcgaTGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAg                                                                                            <  1:965501/151‑1 (MQ=255)
                                          gcgaTGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGc                                                                                           <  1:2874/150‑1 (MQ=255)
                                          gcgaTGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCg                                                                                          >  1:172598/1‑151 (MQ=255)
                                          gcgaTGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCg                                                                                          >  2:631010/1‑151 (MQ=255)
                                              tGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTa                                                                                       <  1:318087/150‑1 (MQ=255)
                                               gCGGTCCGGATGGATCGCCTCCAGCGTCTCAGTTGTCACGGCAGCGCCGACGAATTTGGCCGGTTTCCAGGTAATGTGACGAATAGAGAGCACGGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGc                                                                                           <  2:221153/145‑1 (MQ=255)
                                               gCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCCCCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTaa                                                                                      <  2:667903/150‑1 (MQ=255)
                                                 ggTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGc                                                                                   >  2:906770/1‑151 (MQ=255)
                                                      ggATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCAt                                                                              >  1:810029/1‑151 (MQ=255)
                                                        aTGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTg                                                                            <  1:157394/151‑1 (MQ=255)
                                                        aTGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTg                                                                            <  1:989401/151‑1 (MQ=255)
                                                          ggAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGa                                                                           >  2:301473/1‑150 (MQ=255)
                                                               ggCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCtt                                                                     >  2:475289/1‑151 (MQ=255)
                                                               ggCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCtt                                                                     >  1:892691/1‑151 (MQ=255)
                                                                gCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCttt                                                                    >  2:613276/1‑151 (MQ=255)
                                                                   ccAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTg                                                                  <  2:640208/150‑1 (MQ=255)
                                                                   ccAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGACGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCtttt                                                                   <  2:209405/149‑1 (MQ=255)
                                                                          cTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTg                                                          >  1:876984/1‑151 (MQ=255)
                                                                                 ttCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCa                                                    <  2:532114/150‑1 (MQ=255)
                                                                                   ctcGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGCGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATcc                                                 >  2:44277/1‑151 (MQ=255)
                                                                                      gCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGa                                              <  1:246861/151‑1 (MQ=255)
                                                                                        cAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGAcc                                            >  1:925028/1‑151 (MQ=255)
                                                                                           cacCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCgg                                          >  2:5371/1‑150 (MQ=255)
                                                                                             ccGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGt                                         <  1:499302/149‑1 (MQ=255)
                                                                                             ccGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTc                                        >  2:598125/1‑150 (MQ=255)
                                                                                             ccGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCa                                       <  2:988540/151‑1 (MQ=255)
                                                                                              cGAGGAACTTGATCGGTTTGCCAGTGAGCTGACGAATAGAGAGCCCCTCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACTCTCCGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCAt                                      <  1:199976/151‑1 (MQ=255)
                                                                                                    aCTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTc                                <  2:529167/151‑1 (MQ=255)
                                                                                                      ttGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCa                               >  2:730977/1‑150 (MQ=255)
                                                                                                       tGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCa                               <  2:325173/149‑1 (MQ=255)
                                                                                                        gATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAAcc                            <  1:956425/151‑1 (MQ=255)
                                                                                                        gATCGGTTTGCCAGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAAcc                            <  1:36751/151‑1 (MQ=255)
                                                                                                                    aGTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGtt                >  1:684996/1‑151 (MQ=255)
                                                                                                                     gTGATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGttt               <  2:695161/151‑1 (MQ=255)
                                                                                                                       gATGTGACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTc              <  2:277125/150‑1 (MQ=255)
                                                                                                                            gACGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAACCgg        >  1:916937/1‑151 (MQ=255)
                                                                                                                              cGAATAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAACCGGGt      <  1:10425/151‑1 (MQ=255)
                                                                                                                                aaTAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAACCGGGtt     >  1:787198/1‑150 (MQ=255)
                                                                                                                                aaTAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAACCGGGtt     >  1:42648/1‑150 (MQ=255)
                                                                                                                                  tAGAGAGCACCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAACCGGGTTAAt  >  2:142393/1‑151 (MQ=255)
                                                                                                                                          |                                                                                                                                              
TCGATCAGCGACAGTACGTCGCCCATGCCGAGAATACGCGACGCGATGCGGTCCGGATGGAACGGCTCCAGCGCCTCAGTTTTCTCGCCAACACCGAGGAACTTGATCGGTTTGCCAGTGATGTGACGAATAGAGAGCGCCGCACCGCCGCGGGCATCGCCGTCCACTTTGGTCAACACTACGCCGGTAAGCGGTAACGCTTCATTGAATGCTTTTGCCGTATTGGCCGCATCCTGACCGGTCATGGCGTCAACCACAAACAGGGTTTCAACCGGGTTAAT  >  NC_000913/2746882‑2747162

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