Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 97,073 Δ1 bp 7.1% coding (1072/1083 nt) mraY → phospho‑N‑acetylmuramoyl‑pentapeptide transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091397,0730A. 7.1% 185.4 / 7.0 56coding (1072/1083 nt)mraYphospho‑N‑acetylmuramoyl‑pentapeptide transferase
Reads supporting (aligned to +/- strand):  ref base A (23/29);  new base . (2/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTG  >  NC_000913/96926‑97218
                                                                                                                                                   |                                                                                                                                                 
ggCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGaa                                                                                                                                                  <  7:358035/149‑1 (MQ=255)
 gCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGatg                                                                                                                                                 <  2:65779/149‑3 (MQ=255)
   tCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAgg                                                                                                                                                >  3:114587/1‑148 (MQ=255)
     cTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGt                                                                                                                                               <  7:246610/147‑1 (MQ=255)
      tttAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACg                                                                                                                                            >  5:179465/1‑149 (MQ=255)
        tAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGGTTCTGGATTATTTCGCTGATGCTGGGTCTGATTGGTCTGGCAACGCTGAAGGTACg                                                                                                                                            <  6:72195/147‑1 (MQ=255)
         aaaCTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGtt                                                                                                                                          <  7:195722/148‑1 (MQ=255)
          aaCTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTaa                                                                                                                                        <  6:146485/149‑1 (MQ=255)
                gcgGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATgg                                                                                                                                  <  4:3143/149‑1 (MQ=255)
                gcgGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATgg                                                                                                                                  <  3:307375/149‑1 (MQ=255)
                 cgGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGc                                                                                                                                 >  3:314564/1‑149 (MQ=255)
                  ggACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCt                                                                                                                                <  6:259353/149‑1 (MQ=255)
                            ttttCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCTCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAgg                                                                                                                      <  4:114587/149‑1 (MQ=255)
                            ttttCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAgg                                                                                                                      <  5:7763/149‑1 (MQ=255)
                                ccGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGt                                                                                                                                           <  4:57131/124‑1 (MQ=255)
                                ccGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGt                                                                                                                                           >  3:57131/1‑124 (MQ=255)
                                        cACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATgtcg                                                                                                          >  5:121935/1‑149 (MQ=255)
                                         aCCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATgtcgt                                                                                                         <  6:179465/149‑1 (MQ=255)
                                         aCCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATgtcgt                                                                                                         <  7:387354/149‑1 (MQ=255)
                                                 atcaCCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATg                                                                                                                                   <  5:32990/115‑1 (MQ=255)
                                                 atcaCCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATg                                                                                                                                   >  6:32990/1‑115 (MQ=255)
                                                 atcaCCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCa                                                                                                                        <  8:40827/126‑1 (MQ=255)
                                                 atcaCCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCa                                                                                                                        >  7:40827/1‑126 (MQ=255)
                                                 atcaCCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAGGTACGTTAATCATGGCTGATTATCAGGGTATAAATGTCGTCATTATCgg                                                                                                >  2:325989/1‑149 (MQ=255)
                                                    accacTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGcc                                                                                              >  3:12671/1‑149 (MQ=255)
                                                     ccacTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTg                                                                                            >  6:74271/1‑149 (MQ=255)
                                                      cacTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTgg                                                                                           <  5:74271/149‑1 (MQ=255)
                                                           tGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGc                                                                                         <  3:218702/147‑1 (MQ=255)
                                                              aCTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTTTAAGGGTAAAAATGTCGTCATTATCGGACTGGGCCacac                                                                                    >  8:168414/1‑149 (MQ=255)
                                                              aCTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCAc                                                                                    <  6:8771/149‑1 (MQ=255)
                                                              aCTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGCTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCAc                                                                                    >  5:8771/1‑149 (MQ=255)
                                                                  aaaGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCggg                                                                                >  5:319743/1‑149 (MQ=255)
                                                                      gCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCttt                                                                            <  4:314564/149‑1 (MQ=255)
                                                                           ccGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGa                                                                                                                                                   >  8:260034/1‑73 (MQ=255)
                                                                           ccGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGa                                                                                                                                                   <  7:260034/73‑1 (MQ=255)
                                                                            cGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTAACGGCCTGGGCCTCACCGGGCTTTCCTGCg                                                                      >  1:9417/1‑149 (MQ=255)
                                                                                    cgcgTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAGGTACGTTAATCATGGCTGATTATCAGGGTATAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACtttt                                                             <  1:325989/149‑1 (MQ=255)
                                                                                      cgTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACttttt                                                            <  6:121935/149‑1 (MQ=255)
                                                                                       gTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTc                                                           >  3:238590/1‑149 (MQ=255)
                                                                                                        gATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACtttt                                                             <  6:249466/130‑1 (MQ=255)
                                                                                                        gATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACtttt                                                             >  5:249466/1‑130 (MQ=255)
                                                                                                        gATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGAc                                          >  7:89984/1‑149 (MQ=255)
                                                                                                        gATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGAc                                          >  7:307137/1‑149 (MQ=255)
                                                                                                          ttattTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGAcgc                                        >  8:195809/1‑149 (MQ=255)
                                                                                                            attTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGAcgcc                                       <  6:319743/148‑1 (MQ=255)
                                                                                                                 gCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGt                                 <  8:307137/149‑1 (MQ=255)
                                                                                                                 gCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGt                                 >  4:246797/1‑149 (MQ=255)
                                                                                                                  cTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGtt                                <  3:246797/149‑1 (MQ=255)
                                                                                                                      tGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGTTATgg                            <  7:195809/149‑1 (MQ=255)
                                                                                                                         tGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATTATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGTTATGGATa                         >  2:31998/1‑149 (MQ=255)
                                                                                                                           gTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTg                                                                    <  3:19722/104‑1 (MQ=255)
                                                                                                                           gTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTg                                                                    >  4:19722/1‑104 (MQ=255)
                                                                                                                           gTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACttttt                                                            >  3:95455/1‑112 (MQ=255)
                                                                                                                           gTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACttttt                                                            <  4:95455/112‑1 (MQ=255)
                                                                                                                                      tCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGTTATGGATACGCGTATGACAcc            >  8:294990/1‑149 (MQ=255)
                                                                                                                                              cGCTGAAGGTACGTTAATTATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCgcctggcc    <  1:31998/149‑1 (MQ=255)
                                                                                                                                               gCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCgcctggcct   <  4:12671/149‑1 (MQ=255)
                                                                                                                                               gCTGAAGGTACGTTAAGCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCgcctggcct   <  7:156317/149‑1 (MQ=255)
                                                                                                                                                cTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTAACGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCgcctggcctg  <  2:9417/149‑1 (MQ=255)
                                                                                                                                                 tGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCgcctggcctg  >  3:308365/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGCTCCTTTAAACTGCGCGGACAACGTATTTTCCGCATGGCACCGATTCATCACCACTATGAACTGAAAGGCTGGCCGGAACCGCGCGTCATTGTGCGTTTCTGGATTATTTCGCTGATGCTGGTTCTGATTGGTCTGGCAACGCTGAAGGTACGTTAATCATGGCTGATTATCAGGGTAAAAATGTCGTCATTATCGGCCTGGGCCTCACCGGGCTTTCCTGCGTGGACTTTTTCCTCGCTCGCGGTGTGACGCCGCGCGTTATGGATACGCGTATGACACCGCCTGGCCTG  >  NC_000913/96926‑97218

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