Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,723,065 G→A A113A (GCG→GCA pssA → phosphatidylserine synthase, CDP‑diacylglycerol‑serine O‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,723,0650GA100.0% 227.4 / NA 69A113A (GCG→GCApssAphosphatidylserine synthase, CDP‑diacylglycerol‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/34);  total (35/34)

ACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGT  >  NC_000913/2722932‑2723199
                                                                                                                                     |                                                                                                                                      
aCGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGa                                                                                                                                   <  2:475372/139‑1 (MQ=255)
    gTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAAt                                                                                                                               >  1:1433470/1‑139 (MQ=255)
     ttGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTTGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATc                                                                                                                              <  2:1336706/139‑1 (MQ=255)
       gTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCg                                                                                                                            >  2:1335937/1‑139 (MQ=255)
         aTGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCggg                                                                                                                          <  1:413529/139‑1 (MQ=255)
         aTGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCggg                                                                                                                          >  2:1471250/1‑139 (MQ=255)
              gCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCggg                                                                                                                          <  1:538197/134‑1 (MQ=255)
              gCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCggg                                                                                                                          >  2:538197/1‑134 (MQ=255)
                tAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCggg                                                                                                                          >  2:321081/1‑132 (MQ=255)
                tAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCggg                                                                                                                          <  1:321081/132‑1 (MQ=255)
                  aaaGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGt                                                                                                                 >  2:161293/1‑139 (MQ=255)
                  aaaGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGt                                                                                                                 >  1:1125167/1‑139 (MQ=255)
                        aGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGt                                                                                                           >  1:835805/1‑139 (MQ=255)
                            tCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTAt                                                                                                       >  2:934050/1‑139 (MQ=255)
                             ccGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATg                                                                                                      <  2:1493253/139‑1 (MQ=255)
                                gAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCg                                                                                                   <  2:1211656/139‑1 (MQ=255)
                                  aCTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGtt                                                                                                 <  1:1257246/139‑1 (MQ=255)
                                      gATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaat                                                                                            <  2:544662/139‑1 (MQ=255)
                                       aTGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaat                                                                                            <  1:103729/139‑1 (MQ=255)
                                        tgtgCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCgg                                                                                                                           >  2:630135/1‑107 (MQ=255)
                                        tgtgCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCgg                                                                                                                           <  1:630135/107‑1 (MQ=255)
                                           gCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaatcaat                                                                                        >  1:969316/1‑139 (MQ=255)
                                                      tCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCTGGGCGTAGATGTTCCGGTTTa                                                                                                        <  2:902569/112‑1 (MQ=255)
                                                      tCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCTGGGCGTAGATGTTCCGGTTTa                                                                                                        >  1:902569/1‑112 (MQ=255)
                                                          cTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATcc                                                                                                                             <  1:1128625/87‑1 (MQ=255)
                                                          cTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATcc                                                                                                                             >  2:1128625/1‑87 (MQ=255)
                                                             gCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCtt                                                                         <  2:383731/136‑1 (MQ=255)
                                                             gCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCtt                                                                         >  1:383731/1‑136 (MQ=255)
                                                                         aCGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACttt                                                          >  2:630000/1‑139 (MQ=255)
                                                                          cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGtt                                                                                                 <  1:1342885/99‑1 (MQ=255)
                                                                          cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGtt                                                                                                 >  2:1342885/1‑99 (MQ=255)
                                                                          cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACttt                                                          <  2:1589582/138‑1 (MQ=255)
                                                                          cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACttt                                                          >  1:1589582/1‑138 (MQ=255)
                                                                          cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         >  1:1602995/1‑139 (MQ=255)
                                                                           gTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTaa                                                        <  2:1433470/139‑1 (MQ=255)
                                                                             ggACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTg                                                                        <  1:626235/121‑1 (MQ=255)
                                                                             ggACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTg                                                                        >  2:626235/1‑121 (MQ=255)
                                                                                    ttGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTat                                               >  2:679460/1‑139 (MQ=255)
                                                                                         gcTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATcg                                          >  1:66255/1‑139 (MQ=255)
                                                                                             cGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGt                                                                                                           >  1:894641/1‑70 (MQ=255)
                                                                                             cGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGt                                                                                                           <  2:894641/70‑1 (MQ=255)
                                                                                             cGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTa                                                                                                        >  1:1204456/1‑73 (MQ=255)
                                                                                             cGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTa                                                                                                        <  2:1204456/73‑1 (MQ=255)
                                                                                             cGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGt                                                                                  >  1:89189/1‑95 (MQ=255)
                                                                                             cGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGt                                                                                  <  2:89189/95‑1 (MQ=255)
                                                                                             cGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATcgacga                                      >  1:649124/1‑139 (MQ=255)
                                                                                               gCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATa                                    <  1:934050/139‑1 (MQ=255)
                                                                                                  tCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCg                                 <  1:1471250/139‑1 (MQ=255)
                                                                                                  tCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCg                                 <  1:1335937/139‑1 (MQ=255)
                                                                                                  tCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCg                                 <  2:835805/139‑1 (MQ=255)
                                                                                                         cTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGTGCGTAGATGCTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCg                                 >  2:1021079/1‑132 (MQ=255)
                                                                                                         cTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCg                                 <  1:1021079/132‑1 (MQ=255)
                                                                                                          tAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTtat                         >  2:880966/1‑139 (MQ=255)
                                                                                                            aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAg                       <  1:801021/139‑1 (MQ=255)
                                                                                                               cTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCgg                    >  1:970131/1‑139 (MQ=255)
                                                                                                                tGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         <  2:472968/101‑1 (MQ=255)
                                                                                                                tGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         >  1:472968/1‑101 (MQ=255)
                                                                                                                   cTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGcc                <  1:630000/139‑1 (MQ=255)
                                                                                                                      gTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGc                                  >  2:597392/1‑118 (MQ=255)
                                                                                                                      gTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGc                                  <  1:597392/118‑1 (MQ=255)
                                                                                                                            ccGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTg                                                                        >  1:570838/1‑74 (MQ=255)
                                                                                                                            ccGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTg                                                                        <  2:570838/74‑1 (MQ=255)
                                                                                                                            ccGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCg                                 <  2:278368/113‑1 (MQ=255)
                                                                                                                            ccGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCg                                 >  1:278368/1‑113 (MQ=255)
                                                                                                                               cATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTg                                                                                 >  2:1561872/1‑62 (MQ=255)
                                                                                                                               cATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTg                                                                                 <  1:1561872/62‑1 (MQ=255)
                                                                                                                                 tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTg                                                                >  2:717701/1‑77 (MQ=255)
                                                                                                                                 tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTg                                                                <  1:717701/77‑1 (MQ=255)
                                                                                                                                 tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGt  >  1:1140412/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                      
ACGCGTTGTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGT  >  NC_000913/2722932‑2723199

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