Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,968,712 G→C 9.4% A50A (GCC→GCG gph ← phosphoglycolate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,968,7120GC 9.4% 157.6 / 8.5 55A50A (GCC→GCGgphphosphoglycolate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (26/23);  new base C (2/3);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

CGATGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTC  >  NC_002947/1968564‑1968851
                                                                                                                                                    |                                                                                                                                           
cGATGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCggg                                                                                                                                            <  2:571241/150‑1 (MQ=255)
  aTGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGgcg                                                                                                                                          >  2:779485/1‑150 (MQ=255)
     tCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCa                                                                                                                                       <  1:668042/150‑1 (MQ=255)
       ggCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCagc                                                                                                                                     >  1:1214646/1‑148 (MQ=255)
       ggCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCCCGCGCCatg                                                                                                                                     >  1:106055/1‑148 (MQ=255)
         ccAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCCCGCGCCatggt                                                                                                                                   <  2:880668/150‑5 (MQ=255)
          cagcagCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCCCGCGCCatggtt                                                                                                                                  <  2:106055/150‑6 (MQ=255)
            gcagcTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCgaga                                                                                                                                >  1:216149/1‑150 (MQ=255)
              agcTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCCCGGCGCGAGCACCCCCCGCGATc                                                                                                                              >  2:332498/1‑150 (MQ=255)
               gcTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCggg                                                                                                                                            >  2:115136/1‑135 (MQ=255)
               gcTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCggg                                                                                                                                            <  1:115136/135‑1 (MQ=255)
                 tCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCCCGCGCCatgg                                                                                                                                    >  2:500307/1‑138 (MQ=255)
                 tCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCCCGCGCCatgg                                                                                                                                    <  1:500307/141‑4 (MQ=255)
                           ccTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCa                                                                                                                 <  1:1070996/150‑1 (MQ=255)
                             tCGAACAGCTTGCTGGGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCGAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAgg                                                                                                               <  1:452209/150‑1 (MQ=255)
                               gAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGtt                                                                                                             <  2:200849/150‑1 (MQ=255)
                                aaCAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTg                                                                                                            >  1:593948/1‑150 (MQ=255)
                                aaCAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTAACCCTCCGCCGCTCGCGCGGAGGCGGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTg                                                                                                            >  1:610787/1‑150 (MQ=255)
                                         cTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAAcg                                                                                                   <  2:216149/150‑1 (MQ=255)
                                          tgtgCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGtt                                                                                                             <  1:465067/139‑1 (MQ=255)
                                          tgtgCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGtt                                                                                                             >  2:465067/1‑139 (MQ=255)
                                            tgCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAAcggcc                                                                                                <  1:580818/150‑1 (MQ=255)
                                                 cgcgCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAg                                                                                           >  1:917432/1‑150 (MQ=255)
                                                 cgcgAAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAg                                                                                           >  2:720383/1‑150 (MQ=255)
                                                    gcAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGTAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAAcggccgg                                                                                              >  1:1057205/1‑144 (MQ=255)
                                                    gcAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGTAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAAcggccgg                                                                                              <  2:1057205/144‑1 (MQ=255)
                                                                gTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAg                                                                            >  2:340777/1‑150 (MQ=255)
                                                                       tCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGcc                                                                     >  1:203889/1‑150 (MQ=255)
                                                                         gAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAAcg                                                                                                   <  1:171137/118‑1 (MQ=255)
                                                                         gAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAAcg                                                                                                   >  2:171137/1‑118 (MQ=255)
                                                                          aGGAACTCCAGGCGCAGGGGCTCGAAGCCGTAAGAATCGGGATCCGTGGCGGACGATAAGGGAACCATGGCGCGGGCGGCACCCGAGATCACGCCACGCATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTgg                                                                  >  2:468934/1‑150 (MQ=255)
                                                                              aCTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCaaa                                                              >  2:675658/1‑150 (MQ=255)
                                                                                 ccAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATgg                                                           <  2:1224625/150‑1 (MQ=255)
                                                                                  cAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGc                                                          >  2:701296/1‑150 (MQ=255)
                                                                                   aGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCg                                                         >  1:879619/1‑150 (MQ=255)
                                                                                     gcgcAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCTCGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGAt                                                       <  1:477587/150‑1 (MQ=255)
                                                                                                    gCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCgt                                        <  2:203889/150‑1 (MQ=255)
                                                                                                      cGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGa                                                                             <  2:805281/111‑1 (MQ=255)
                                                                                                      cGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGa                                                                             >  1:805281/1‑111 (MQ=255)
                                                                                                       gTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTc                                     >  2:445624/1‑150 (MQ=255)
                                                                                                         cAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACCGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGGCGTGTCCa                                   <  2:1214646/150‑1 (MQ=255)
                                                                                                           gCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCgtatttgct                                 <  1:701296/150‑8 (MQ=255)
                                                                                                             cTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCaa                               >  1:711712/1‑150 (MQ=255)
                                                                                                              tCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAg                              >  1:914334/1‑150 (MQ=255)
                                                                                                                cGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGt                            <  2:711712/150‑1 (MQ=255)
                                                                                                                    tCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAGGTCCGGCGCCGTGTCCAGCAAGGTGCCg                        <  1:675658/150‑1 (MQ=255)
                                                                                                                               cGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGaa             <  1:445624/150‑1 (MQ=255)
                                                                                                                                    gggCAACCATGCCGCGGCCGCCACCCGAGATCACGCCACGGCTCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTa        <  1:332498/149‑1 (MQ=255)
                                                                                                                                       cAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTg     >  1:554509/1‑150 (MQ=255)
                                                                                                                                        aaCCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTc                                                                                  >  2:102231/1‑72 (MQ=255)
                                                                                                                                        aaCCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTc                                                                                  <  1:102231/72‑1 (MQ=255)
                                                                                                                                        aaCCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCAt                   >  2:113157/1‑135 (MQ=255)
                                                                                                                                        aaCCATGGCGCGGGCGCCACCCGAGATCACGCCACGGAGCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCAt                   <  1:113157/135‑1 (MQ=255)
                                                                                                                                          ccATGGCGCGGGCGCCACCTGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTc  >  2:799069/1‑150 (MQ=255)
                                                                                                                                          ccATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGaa             >  1:217383/1‑139 (MQ=255)
                                                                                                                                          ccATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGaa             <  2:217383/139‑1 (MQ=255)
                                                                                                                                            aTGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGAt                                                       >  1:1213667/1‑95 (MQ=255)
                                                                                                                                            aTGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGAt                                                       <  2:1213667/95‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
CGATGTCGGCCAGCAGCTCGGCCATGCCCTCGAACAGCTTGCTGTGCACCGCGCAGTCGCGCTGGTAACGCTCGAGGAACTCCAGGCGCAGGGCCTCGAAGCCGTCAGCCTCCGGGTCCATGGCGAACGCAGTGGCAACCATGGCGCGGGCGCCACCCGAGATCACGCCACGGATCAGGTTGTCGTCAACGGCCGGCAGCCCGCGTTCGGCGAGCATGGCCTGGCAAATGGCGATGAAGTCCGGCGCCGTGTCCAGCAAGGTGCCGTCCATGTCGAAGAGTACTGCTC  >  NC_002947/1968564‑1968851

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