Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,921,550 T→A 13.8% intergenic (+25/+36) PP_1720 → / ← PP_1721 alcohol dehydrogenase/phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,921,5500TA13.8% 128.4 / 12.7 58intergenic (+25/+36)PP_1720/PP_1721alcohol dehydrogenase/phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base A (6/2);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

CGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAG  >  NC_002947/1921403‑1921691
                                                                                                                                                   |                                                                                                                                               
cGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTa                                                                                                                                                <  5:32551/149‑1 (MQ=255)
        gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGGCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGaa                                                                                                                                        <  1:241680/149‑1 (MQ=255)
            cAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAgggg                                                                                                                                   <  7:193649/149‑1 (MQ=255)
                tGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCgg                                                                                                                               <  4:104945/149‑1 (MQ=255)
                 ggTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTAGCGGGTTTACCCGCGAAGGGGCCGGc                                                                                                                              >  5:240471/1‑149 (MQ=255)
                      acacCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGc                                                                                                                        >  3:208622/1‑149 (MQ=255)
                      acacCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGc                                                                                                                        >  3:285147/1‑149 (MQ=255)
                      acacCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGc                                                                                                                        >  7:236737/1‑149 (MQ=255)
                       cacCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAggg                                                                                                                                    <  6:208708/137‑1 (MQ=255)
                                  tGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCt                                                                                                            <  8:236737/149‑1 (MQ=255)
                                   ggAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                           <  3:86932/149‑1 (MQ=255)
                                    tcgTTTGCCCCGGCGGCGTATGCCGAGCTGGAGACCAACCCGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTACGGGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTc                                                                                                          <  4:285147/147‑1 (MQ=255)
                                    taGTTTGCCCAGGCGGCTTATGCCGAGATGGAGACCAACCAGGTGTCGGGAAAGCTGGTGATGTTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGTTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTc                                                                                                          <  4:200636/148‑1 (MQ=255)
                                      gTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcag                                                                                                        >  3:45603/1‑149 (MQ=255)
                                      gTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTAGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTc                                                                                                          >  3:200636/1‑147 (MQ=255)
                                       tttGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagc                                                                                                       >  7:255265/1‑149 (MQ=255)
                                           cccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctcc                                                                                                    <  6:309959/148‑1 (MQ=255)
                                           cccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctc                                                                                                     >  7:192646/1‑147 (MQ=255)
                                            ccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGTGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctccag                                                                                                  >  1:4715/1‑149 (MQ=255)
                                            ccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGaa                                                                                                                                        <  5:136785/113‑1 (MQ=255)
                                            ccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGaa                                                                                                                                        >  6:136785/1‑113 (MQ=255)
                                                      tATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCg                                                                                        <  8:192646/149‑1 (MQ=255)
                                                        tGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACACGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGc                                                                                      >  3:68050/1‑149 (MQ=255)
                                                          ccGAGCTGGATACAACCCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGcc                                                                                     <  8:255265/148‑1 (MQ=255)
                                                           tgAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCAc                                                                                   >  2:96484/2‑149 (MQ=255)
                                                           cGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCAc                                                                                   <  4:208622/149‑1 (MQ=255)
                                                                  gagaCCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCgg                                                                                       >  1:303059/1‑138 (MQ=255)
                                                                  gagaCCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGa                                                                            >  4:334041/1‑149 (MQ=255)
                                                                  gagaCCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGa                                                                            >  4:344184/1‑149 (MQ=255)
                                                                  gagaCCAACATGGTGTGGGGCAAGGTGGTGATGTTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTGGCGGGTTTACCGGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCgg                                                                                       <  2:303059/138‑1 (MQ=255)
                                                                    gaCCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAgggg                                                                                                                                   >  6:2232/1‑93 (MQ=255)
                                                                           caGGTGGCGGGCATGCTGGTGATGGTGATCGCCCCTAGCCGGGTGAAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCAt                                                                   <  8:138264/149‑1 (MQ=255)
                                                                                gtCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTcg                                                              >  7:27444/1‑149 (MQ=255)
                                                                                 tCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGTGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTcgc                                                             <  2:4715/149‑1 (MQ=255)
                                                                                  cGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTcgcg                                                            <  3:334041/149‑1 (MQ=255)
                                                                                   gggCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctccagctc                                                                                               >  4:151647/1‑113 (MQ=255)
                                                                                   gggCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctccagctc                                                                                               <  3:151647/113‑1 (MQ=255)
                                                                                    ggCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCATCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCATCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCCGCAAGCTcgccc                                                           >  4:131610/1‑146 (MQ=255)
                                                                                     gCAAGGTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCaa                                                         <  3:131610/149‑1 (MQ=255)
                                                                                                        cGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGc                                                                 <  8:188617/122‑1 (MQ=255)
                                                                                                        cGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGc                                                                 >  7:188617/1‑122 (MQ=255)
                                                                                                          aCCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGAGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcacca                                     >  1:91706/1‑149 (MQ=255)
                                                                                                             cTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCgg                                                                                       >  5:103051/1‑95 (MQ=255)
                                                                                                                 ccTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCCCCCATTGATGTCGGCATGCTCGCGCAACACTGCGTCCGGGTTGACCACATGCcgg                             >  8:224239/1‑149 (MQ=255)
                                                                                                                 ccTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCggg                             >  2:63300/1‑149 (MQ=255)
                                                                                                                   tGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGt                                                 <  4:163602/127‑1 (MQ=255)
                                                                                                                   tGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGt                                                 >  3:163602/1‑127 (MQ=255)
                                                                                                                      tgtAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATc                        <  2:91706/149‑1 (MQ=255)
                                                                                                                          aGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCATCACTGCGTCCGGGTTCACCACATGCGGGTGATCCAct                    <  1:63300/149‑1 (MQ=255)
                                                                                                                                    gCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTcaacaa            <  8:27444/147‑1 (MQ=255)
                                                                                                                                      cccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCaacaacaac        >  8:238655/1‑149 (MQ=255)
                                                                                                                                       ccTTCGCGGGTAAACCCGCGAAGGGGCCAGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACg       >  5:3105/1‑149 (MQ=255)
                                                                                                                                       ccTTCGCGGGTAAACCCGCGAAGGGGCCAGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACg       >  1:122524/1‑149 (MQ=255)
                                                                                                                                        cTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACg       >  1:47626/1‑149 (MQ=255)
                                                                                                                                        cTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCggg                                            >  6:298248/1‑111 (MQ=255)
                                                                                                                                        cTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCggg                                            <  5:298248/111‑1 (MQ=255)
                                                                                                                                         ttCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACgg      >  6:269504/1‑148 (MQ=255)
                                                                                                                                             cgGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAg  >  4:306812/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAG  >  NC_002947/1921403‑1921691

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