Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,391 +C coding (474/516 nt) b4555 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,3911.C100.0% 197.2 / NA 56K158N (AAG→AACb4555hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (18/38);  total (18/38)

CTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008300‑4008458
                                                                                              |                                                                            
cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAg                                                                               >  1:616201/1‑93 (MQ=255)
cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAg                                                                               <  1:5443504/93‑1 (MQ=255)
cTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAg                                                                               >  2:1915592/1‑93 (MQ=255)
 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTTGGAgg                                                                             >  1:5331944/1‑93 (MQ=255)
 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                              <  2:4835183/93‑2 (MQ=255)
 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                              <  2:2367467/93‑2 (MQ=255)
 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                              <  1:1823767/93‑2 (MQ=255)
 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                              <  1:61034/93‑2 (MQ=255)
 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                              <  1:4510488/93‑2 (MQ=255)
 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgg                                                                             >  1:1238839/1‑93 (MQ=255)
 tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGGAgg                                                                             >  2:246032/1‑93 (MQ=255)
  cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcg                                                                             <  1:1055401/93‑3 (MQ=255)
  cGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcg                                                                             <  1:4583644/93‑3 (MQ=255)
   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                            >  1:378179/1‑90 (MQ=255)
   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                            <  2:4594730/93‑4 (MQ=255)
   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcgg                                                                            <  1:4328415/93‑4 (MQ=255)
   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGGGCTGGGAgcgg                                                                            >  2:3886453/1‑90 (MQ=255)
    ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGGCATGCGCTGGAAgcggg                                                                           >  2:4431194/1‑89 (MQ=255)
    ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                           <  1:5689436/93‑5 (MQ=255)
    ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                           <  1:4658164/93‑5 (MQ=255)
    ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                           <  2:2513188/93‑5 (MQ=255)
    ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                           <  2:3513598/93‑5 (MQ=255)
         gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                  <  2:1055363/93‑1 (MQ=255)
           gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                >  2:5666223/1‑92 (MQ=255)
           gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                >  1:2245207/1‑92 (MQ=255)
           gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                >  1:383212/1‑92 (MQ=255)
            aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCCAAAGGCCATGCGCTGGGAGCGGTGGCG‑‑‑‑aaga                                                               >  2:5218446/1‑91 (MQ=255)
            aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               <  1:4210164/93‑3 (MQ=255)
            aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               >  2:3437448/1‑91 (MQ=255)
            aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGAGATTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               >  2:5286486/1‑91 (MQ=255)
             tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                              <  2:1335952/93‑4 (MQ=255)
             tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                              >  1:6051407/1‑90 (MQ=255)
             tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                              <  1:2952901/93‑4 (MQ=255)
                                              tCTGGAAGTGATGGCGGGCGGCGTTTCAAAATGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  1:1901182/93‑1 (MQ=255)
                                              tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 >  2:3190494/1‑93 (MQ=255)
                                              tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  2:2842303/93‑1 (MQ=255)
                                               cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                <  1:4465851/93‑2 (MQ=255)
                                               cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                <  1:2334120/93‑2 (MQ=255)
                                                tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:4989388/93‑3 (MQ=255)
                                                tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:1737083/93‑3 (MQ=255)
                                                  gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  2:226086/93‑5 (MQ=255)
                                                  gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  2:6069960/93‑5 (MQ=255)
                                                  gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  1:4817457/93‑5 (MQ=255)
                                                         tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  1:2302993/93‑1 (MQ=255)
                                                         tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  2:2950850/93‑1 (MQ=255)
                                                         tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  2:4924566/93‑1 (MQ=255)
                                                          ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  1:1022573/93‑1 (MQ=255)
                                                          ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  2:2143856/93‑1 (MQ=255)
                                                          ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  2:743236/93‑1 (MQ=255)
                                                            ggggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                   <  2:3840682/92‑1 (MQ=255)
                                                            cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                   <  1:1615914/93‑1 (MQ=255)
                                                            cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                   <  2:3318504/93‑1 (MQ=255)
                                                             gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  >  1:3361384/1‑93 (MQ=255)
                                                             gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  2:1588234/93‑1 (MQ=255)
                                                              ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  2:4580266/93‑1 (MQ=255)
                                                              ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  2:4815600/93‑1 (MQ=255)
                                                              ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  1:2267000/93‑1 (MQ=255)
                                                               gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGCTTGTATTGCGTTTGGTGACGGGATGAacg                >  2:4396292/1‑93 (MQ=255)
                                                               gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                <  2:2466906/93‑1 (MQ=255)
                                                               gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                <  2:4018872/93‑1 (MQ=255)
                                                                cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  1:4477690/93‑1 (MQ=255)
                                                                cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  2:670750/93‑1 (MQ=255)
                                                                cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  2:4254565/93‑1 (MQ=255)
                                                                cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAGGAacga               <  1:5107890/93‑1 (MQ=255)
                                                                 ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGTGTTTGGTGACGGGAGGAacgac              >  2:2158791/1‑93 (MQ=255)
                                                                  gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  1:3018123/1‑93 (MQ=255)
                                                                  gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  2:2336717/1‑93 (MQ=255)
                                                                   cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            >  2:3055681/1‑93 (MQ=255)
                                                                   cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            <  1:4718515/93‑1 (MQ=255)
                                                                    gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           <  2:6070951/93‑1 (MQ=255)
                                                                        tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        >  2:1177485/4‑93 (MQ=255)
                                                                         caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:3624707/91‑1 (MQ=255)
                                                                         caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       >  2:4662578/3‑93 (MQ=255)
                                                                         caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:824842/91‑1 (MQ=255)
                                                                          aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:4072808/2‑93 (MQ=255)
                                                                          aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:4803706/2‑93 (MQ=255)
                                                                          aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:2985712/2‑93 (MQ=255)
                                                                          aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:1681157/2‑93 (MQ=255)
                                                                          aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      <  1:2804031/92‑1 (MQ=255)
                                                                          aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:5024781/2‑93 (MQ=255)
                                                                           aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAtc     >  2:4531756/1‑92 (MQ=255)
                                                                           aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  2:662695/1‑93 (MQ=255)
                                                                             aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:4274699/93‑1 (MQ=255)
                                                                             aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:910103/93‑1 (MQ=255)
                                                                              ggCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  >  2:5616454/1‑93 (MQ=255)
                                                                                              |                                                                            
CTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAGGCCATGCGCTGGAAGGCTGGCG‑‑‑‑AACTGGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008300‑4008458

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