Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,367,598 C→A E244* (GAG→TAG)  aspA ← aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,367,5980CA100.0% 209.0 / NA 59E244* (GAG→TAG) aspAaspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (31/28);  total (31/28)

CCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGT  >  NC_000913/4367466‑4367733
                                                                                                                                    |                                                                                                                                       
ccAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGc                                                                                                                                  <  2:994085/140‑1 (MQ=255)
ccAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGc                                                                                                                                  <  2:783766/140‑1 (MQ=255)
ccAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGc                                                                                                                                  <  2:1161475/140‑1 (MQ=255)
ccAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGc                                                                                                                                  <  1:657868/140‑1 (MQ=255)
  aGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCgt                                                                                                                                >  2:1153043/1‑140 (MQ=255)
  aGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCgt                                                                                                                                >  1:637596/1‑140 (MQ=255)
   ggCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGc                                                                                                                                  >  2:890538/1‑137 (MQ=255)
   ggCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGc                                                                                                                                  <  1:890538/137‑1 (MQ=255)
     cGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGtt                                                                                                                             <  1:396791/140‑1 (MQ=255)
          cAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGAc                                                                                                                        <  1:983000/140‑1 (MQ=255)
            gcgcgcCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCa                                                                                                                      <  1:721979/140‑1 (MQ=255)
              gcgcCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGt                                                                                                                    >  1:983355/1‑140 (MQ=255)
              gcgcCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGt                                                                                                                    <  2:498917/140‑1 (MQ=255)
                      aaCCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTg                                                                                                            <  2:351055/140‑1 (MQ=255)
                           taacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGCACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGtt                                                                                                       <  2:699670/140‑1 (MQ=255)
                            aacataaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTTACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTg                                                                                                      >  1:674621/1‑140 (MQ=255)
                                taaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCAcc                                                                                                  >  2:608366/1‑140 (MQ=255)
                                taaGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCAcc                                                                                                  >  2:340054/1‑140 (MQ=255)
                                   gcgcCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCTGACCAGTACCGATTGCTGTTTCACCAAg                                                                                               >  2:1092639/1‑140 (MQ=255)
                                   gcgcCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAg                                                                                               >  2:709703/1‑140 (MQ=255)
                                     gccgcAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGt                                                                                             <  2:1191016/140‑1 (MQ=255)
                                        gcAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTaa                                                                                          >  1:570433/1‑140 (MQ=255)
                                        gcAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTaa                                                                                          <  2:983355/140‑1 (MQ=255)
                                         cAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAAc                                                                                         >  1:1252678/1‑140 (MQ=255)
                                             cAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTcc                                                                                     <  2:1252678/140‑1 (MQ=255)
                                                agGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCcagc                                                                                  <  1:608366/140‑1 (MQ=255)
                                                  gTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCcagcag                                                                                <  1:437613/140‑1 (MQ=255)
                                                     gCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGTCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCcagcagcag                                                                             >  1:1191597/1‑140 (MQ=255)
                                                     gCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCcagcagcag                                                                             >  2:146770/1‑140 (MQ=255)
                                                      cTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGt                                                                            >  1:679614/1‑140 (MQ=255)
                                                      cTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGt                                                                            >  2:131248/1‑140 (MQ=255)
                                                       ttCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGtt                                                                           <  2:679614/140‑1 (MQ=255)
                                                             cAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCgg                                                                     <  1:1153043/140‑1 (MQ=255)
                                                                     cAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTgg                                                             <  1:146770/140‑1 (MQ=255)
                                                                              cGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGTTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTc                                                    <  2:137741/140‑1 (MQ=255)
                                                                                  tGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGCACTATTTCGGCGTGTTCAGACCAGTACCGATTCCTGTTGCACCAAGGTTAACTTCCAGCAGCCGTTCAGCGGTACGTTGGATGTTTTTCActt                                                >  2:779602/1‑140 (MQ=255)
                                                                                      aaGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCActtcttc                                            >  1:420870/1‑140 (MQ=255)
                                                                                      aaGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCActtcttc                                            >  1:184851/1‑140 (MQ=255)
                                                                                      aaGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCActtcttc                                            >  1:8924/1‑140 (MQ=255)
                                                                                           aGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTca                                       >  1:1046009/1‑140 (MQ=255)
                                                                                             tAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTcagc                                     <  2:8924/140‑1 (MQ=255)
                                                                                             tAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACATTGGATGTTTTTCACTTCTTCTTTcagc                                     >  1:411682/1‑140 (MQ=255)
                                                                                               aCTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTcagcag                                   >  1:683234/1‑140 (MQ=255)
                                                                                                 ttCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGa                                 <  2:544382/140‑1 (MQ=255)
                                                                                                   cAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCttt                                         >  2:70543/1‑130 (MQ=255)
                                                                                                   cAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCttt                                         <  1:70543/130‑1 (MQ=255)
                                                                                                      ccAGTTTTTTCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTg                            >  1:227338/1‑140 (MQ=255)
                                                                                                          ttttttCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAg                        <  2:184851/140‑1 (MQ=255)
                                                                                                          ttttttCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAg                        <  2:1191597/140‑1 (MQ=255)
                                                                                                          ttttttCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAg                        >  1:260193/1‑140 (MQ=255)
                                                                                                           tttttCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgc                       >  1:422349/1‑140 (MQ=255)
                                                                                                           tttttCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgc                       >  1:527467/1‑140 (MQ=255)
                                                                                                            ttttCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAgcg                      <  2:637596/140‑1 (MQ=255)
                                                                                                              ttCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGcgcg                    <  1:450016/140‑1 (MQ=255)
                                                                                                              ttCACTGCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGcgcg                    <  2:420870/140‑1 (MQ=255)
                                                                                                                    gCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTc              >  2:81731/1‑140 (MQ=255)
                                                                                                                    gCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTc              >  1:757136/1‑140 (MQ=255)
                                                                                                                    gCCAGCGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTc              >  1:1125708/1‑140 (MQ=255)
                                                                                                                         cGGAGAGTACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGAc         >  1:975817/1‑140 (MQ=255)
                                                                                                                                tACTATTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGt  <  2:975817/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
CCAGGCGTTTCAGCGCGCCGTGAACCATAACATAAGCGCCGCAGTCAGAGGTCGCTTCGATCAGGTCTTCAGCCGGTACGCATGGGAAGCCAGTAACTTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGT  >  NC_000913/4367466‑4367733

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