Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,524,569 T→C S567P (TCC→CCC)  mrcA → peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,524,5690TC100.0% 209.8 / NA 59S567P (TCC→CCC) mrcApeptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/28);  total (31/28)

TGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCC  >  NC_000913/3524432‑3524703
                                                                                                                                         |                                                                                                                                      
tGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGccc                                                                                                                                      >  1:455653/1‑140 (MQ=255)
 gggTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCt                                                                                                                                     >  2:1313820/1‑140 (MQ=255)
 gggTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCt                                                                                                                                     <  2:84246/140‑1 (MQ=255)
 gggTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCt                                                                                                                                     >  1:84246/1‑140 (MQ=255)
 gggTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCt                                                                                                                                     >  2:692813/1‑140 (MQ=255)
  ggTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCtt                                                                                                                                    <  1:14315/140‑1 (MQ=255)
  ggTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCtt                                                                                                                                    <  2:6514/140‑1 (MQ=255)
   gTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTc                                                                                                                                   <  1:1502022/140‑1 (MQ=255)
     cAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCAc                                                                                                                                 >  2:1257765/1‑140 (MQ=255)
     cAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCAc                                                                                                                                 >  2:1227036/1‑140 (MQ=255)
     cAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCAc                                                                                                                                 <  1:26834/140‑1 (MQ=255)
        tCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCAcccc                                                                                                                              >  1:352788/1‑140 (MQ=255)
        tCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCCTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACACTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCAcccc                                                                                                                              <  2:1117132/140‑1 (MQ=255)
            aaaaCGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATg                                                                                                                          <  1:1091958/140‑1 (MQ=255)
                         ggTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGcgcgcg                                                                                                             <  1:532745/140‑1 (MQ=255)
                           tACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGc                                                                                                           >  1:396176/1‑140 (MQ=255)
                             cgcgCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTa                                                                                                         >  2:139014/1‑140 (MQ=255)
                                gcAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTAcgc                                                                                                      <  2:352788/140‑1 (MQ=255)
                                    tGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTc                                                                                                  >  1:336642/1‑140 (MQ=255)
                                              gggCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAAcg                                                                                        >  1:205627/1‑140 (MQ=255)
                                              gggCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAAcg                                                                                        >  2:214998/1‑140 (MQ=255)
                                               ggCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAAcgg                                                                                       <  2:455653/140‑1 (MQ=255)
                                                gCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAAcggc                                                                                      >  1:1623310/1‑140 (MQ=255)
                                                   tCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCTTGGCGAACggcggc                                                                                   >  2:1409665/1‑140 (MQ=255)
                                                   tCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACggcggc                                                                                   <  2:220329/140‑1 (MQ=255)
                                                     gACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCtt                                                                                 >  1:1176690/1‑140 (MQ=255)
                                                          cGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCtgg                                                                            >  1:771302/1‑140 (MQ=255)
                                                          cGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCtgg                                                                            <  2:338031/140‑1 (MQ=255)
                                                              gCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGa                                                                        <  1:1227036/140‑1 (MQ=255)
                                                                     atCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTgg                                                                 >  2:1332751/1‑140 (MQ=255)
                                                                      tCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGt                                                                >  1:1273507/1‑140 (MQ=255)
                                                                        tGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGttt                                                              <  2:610614/140‑1 (MQ=255)
                                                                         gCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTa                                                             >  2:772789/1‑140 (MQ=255)
                                                                              gCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGc                                                        <  1:1505287/140‑1 (MQ=255)
                                                                                ttCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTACTGGTGGACCCGTGGTTTATCAGCaa                                                      >  2:539503/1‑140 (MQ=255)
                                                                                  cGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCaaaa                                                    >  2:502646/1‑140 (MQ=255)
                                                                                   ggCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAAt                                                   <  1:1300028/140‑1 (MQ=255)
                                                                                    gCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAAtt                                                  >  1:1104605/1‑140 (MQ=255)
                                                                                       tCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGaa                                               <  2:1498835/140‑1 (MQ=255)
                                                                                        cccGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGaaa                                              <  1:141938/140‑1 (MQ=255)
                                                                                         ccGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGaaaa                                             >  2:565116/1‑140 (MQ=255)
                                                                                            gCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGa                                          <  2:205627/140‑1 (MQ=255)
                                                                                            gCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGa                                          <  2:1176690/140‑1 (MQ=255)
                                                                                             cacaAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGAt                                         >  1:474588/1‑140 (MQ=255)
                                                                                                 aaaCATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAgg                                     <  2:474588/140‑1 (MQ=255)
                                                                                                  aaCATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAggg                                    <  2:1104605/140‑1 (MQ=255)
                                                                                                        gTCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGt                              >  1:547831/1‑140 (MQ=255)
                                                                                                         tCCACACCGAATCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGGGGTTTATCGGCAAAATTGAAAACGACAAGggcgggggg                             >  2:32927/1‑136 (MQ=255)
                                                                                                                  aaTCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGaa                    >  1:882638/1‑140 (MQ=255)
                                                                                                                  aaTCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGaa                    <  1:867451/140‑1 (MQ=255)
                                                                                                                   aTCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTATGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAg                   >  1:950205/1‑140 (MQ=255)
                                                                                                                   aTCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAg                   >  1:608180/1‑140 (MQ=255)
                                                                                                                    tCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTATGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGc                  <  2:950205/140‑1 (MQ=255)
                                                                                                                    tCGCTGGCGCTGGGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGc                  <  2:608180/140‑1 (MQ=255)
                                                                                                                     cgctggcgctggGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCg                 <  1:772789/140‑1 (MQ=255)
                                                                                                                        tggcgctggGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGaaa              >  2:614526/1‑140 (MQ=255)
                                                                                                                        tggcgctggGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGaaa              >  2:779128/1‑140 (MQ=255)
                                                                                                                         ggcgctggGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAAc             <  2:882638/140‑1 (MQ=255)
                                                                                                                              tggGTTCAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGaaa        <  1:625140/140‑1 (MQ=255)
                                                                                                                                    cAGCGCCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGcc  >  1:1533087/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                      
TGGGTCAGTCGAAAAACGTGGTGATGGTACGCGCAATGCGGGCGATGGGCGTCGACTACGCTGCAGAATATCTGCAACGCTTCGGCTTCCCGGCACAAAACATTGTCCACACCGAATCGCTGGCGCTGGGTTCAGCGTCCTTCACCCCAATGCAGGTGGCGCGCGGCTACGCGGTCATGGCGAACGGCGGCTTCCTGGTGGACCCGTGGTTTATCAGCAAAATTGAAAACGATCAGGGCGGCGTGATTTTCGAAGCGAAACCGAAAGTAGCC  >  NC_000913/3524432‑3524703

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