Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,221,097 2 bp→GC coding (824‑825/3093 nt) b3076 → evolved beta‑D‑galactosidase alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,221,0970CG98.7% 371.2 / NA 151T275S (ACG→AGC) b3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/2);  major base G (77/70);  minor base T (0/2);  total (77/74)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,221,0980GC98.7% 385.5 / NA 150T275S (ACG→AGCb3076evolved beta‑D‑galactosidase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base C (77/70);  minor base T (0/2);  total (77/73)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

TGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCACGTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATG  >  NC_000913/3221029‑3221167
                                                                    ||                                                                      
tGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCt                                                                 >  2:4277411/1‑76 (MQ=255)
tGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCt                                                                 <  2:1406648/76‑1 (MQ=255)
 gATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtt                                                                <  2:4788546/76‑1 (MQ=255)
   tGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtttt                                                              <  2:4502840/76‑1 (MQ=255)
   tGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtttt                                                              <  2:4764863/76‑1 (MQ=255)
   tGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtttt                                                              <  1:5816548/76‑1 (MQ=255)
   tGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATCTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtttt                                                              <  1:5657014/76‑1 (MQ=255)
   gggCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCtttt                                                              <  2:5595854/75‑1 (MQ=255)
     gCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTAc                                                            <  2:2310896/76‑1 (MQ=255)
      cGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACt                                                           >  2:2486996/1‑76 (MQ=255)
      cGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACt                                                           >  2:361848/1‑76 (MQ=255)
        aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTgc                                                         >  2:2325108/1‑75 (MQ=255)
        aaCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGt                                                         <  2:4573054/76‑1 (MQ=255)
        aaCGCGTGATGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGt                                                         <  2:2187627/76‑1 (MQ=255)
         aCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTc                                                        >  1:4134725/1‑76 (MQ=255)
         aCGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTc                                                        <  2:4443753/76‑1 (MQ=255)
         aCGCGTGGGGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTc                                                        <  2:843926/76‑1 (MQ=255)
          cgcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCg                                                       <  2:1019947/76‑1 (MQ=255)
          cgcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCg                                                       >  1:2543239/1‑76 (MQ=255)
           gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGa                                                      >  1:1024983/1‑76 (MQ=255)
           gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGa                                                      >  1:3197238/1‑76 (MQ=255)
           gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGa                                                      <  1:515413/76‑1 (MQ=255)
           gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGa                                                      >  1:5178587/1‑76 (MQ=255)
           gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGa                                                      >  1:20194/1‑76 (MQ=255)
           gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGa                                                      >  1:1607807/1‑76 (MQ=255)
           gcgTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCCTTGCTTTTACTGTCGa                                                      >  2:5365217/1‑76 (MQ=255)
             gtggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAAc                                                    <  2:4005579/76‑1 (MQ=255)
              gggtgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACa                                                   <  2:5210138/75‑1 (MQ=255)
                 tgCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGcc                                                <  1:1131483/76‑1 (MQ=255)
                  gCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTCTTACTGTCGAACAGCCg                                               >  1:2861489/1‑76 (MQ=255)
                  gCACAGAAGCGCAATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGccn                                               <  1:4464482/76‑2 (MQ=255)
                   cacaTCACCGCCATTGATCATTTGGCAATTGAAAACCGGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgc                                              <  2:1585440/76‑1 (MQ=255)
                   cacaGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgc                                              >  2:1532230/1‑76 (MQ=255)
                   cacaGCAGCGCCATTGATCATTGGGCAATTGAAAAACGGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgc                                              <  2:758999/76‑1 (MQ=255)
                   cacaGCAGCGCAATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCgc                                              <  2:3248307/76‑1 (MQ=255)
                      agcagcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCcgccgc                                           >  2:5813365/1‑76 (MQ=255)
                        cagcGCCATTGATCATTTGGCAATGGAAAAACTGACCAGCGCCGGCTTTGCTTTTACTGTCGAACAGCCGCAGCaa                                         <  2:2799837/76‑1 (MQ=255)
                         agcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGCACAGCCGCCGCAAt                                        >  2:5798793/1‑76 (MQ=255)
                         agcGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAAt                                        >  1:5634792/1‑76 (MQ=255)
                          gcgcCATGGATCATTTGGCAATTGAAAACCTGACCAGCGCCAGCTTTGTTTTTACTGTCGAACAGCCGCAGCAATg                                       <  2:2875315/76‑1 (MQ=255)
                           cgcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATgg                                      >  1:3935147/1‑76 (MQ=255)
                            gcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATggc                                     >  2:4472441/1‑75 (MQ=255)
                            gcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGt                                     >  2:4100853/1‑76 (MQ=255)
                            gcCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGt                                     <  1:245048/76‑1 (MQ=255)
                             ccATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTc                                    >  1:5154287/1‑76 (MQ=255)
                             ccATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTc                                    <  2:439165/76‑1 (MQ=255)
                              cATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTca                                   <  1:986136/76‑1 (MQ=255)
                              cATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTca                                   >  2:3933060/1‑76 (MQ=255)
                               aTTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTccg                                  >  2:703462/1‑74 (MQ=255)
                               aTTGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcag                                  >  1:2812795/1‑76 (MQ=255)
                                ttGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACCGCCGCAGCCATGGCcaga                                 >  2:141283/1‑75 (MQ=255)
                                ttGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagc                                 >  1:1819908/1‑76 (MQ=255)
                                 tGATCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagca                                <  1:3983218/76‑1 (MQ=255)
                                 tGATCATTTGACAATTGAAAAACTGCCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagca                                <  1:2704475/76‑1 (MQ=255)
                                  gATCATTTGGAAATTGAAAAACTGCCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTcagcag                               <  1:1531837/76‑1 (MQ=255)
                                   aTCATTTGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGa                              >  1:1745263/1‑76 (MQ=255)
                                   aTCATTTGGAAAGTGAAAAACGGACCAGCGCCGGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGa                              <  2:2037227/76‑1 (MQ=255)
                                     cATTTGGCAATTGAAAAACTGCCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAAt                            <  1:1032591/76‑1 (MQ=255)
                                     cATTGGGCAATTAAAAAACTGCCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAAt                            <  1:4897402/76‑1 (MQ=255)
                                       tttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATcc                          >  2:3823441/1‑76 (MQ=255)
                                       ttgggCAATGGAAAAACTGCCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATcc                          <  2:5259906/73‑1 (MQ=255)
                                        ttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATccc                         >  2:4390641/1‑76 (MQ=255)
                                        ttGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGAAATGGTCAGCAGAATcac                         >  2:1589393/1‑74 (MQ=255)
                                         tGGCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGCCAGCAGAATcccc                        >  2:1342796/1‑76 (MQ=255)
                                          ggCATTTGAAAACCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCt                       <  1:1521777/76‑1 (MQ=255)
                                          ggCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCt                       >  1:1194714/1‑76 (MQ=255)
                                          ggCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCt                       >  1:2029531/1‑76 (MQ=255)
                                          ggCAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCt                       <  1:3967508/76‑1 (MQ=255)
                                          ggCAATTGAAAAACTGACCAGCGCCAGCTGTGCTTTTACTCTCGAACAGCCGCAGCAATGGGCAGCAGAATcaccc                       >  1:2318359/1‑75 (MQ=255)
                                          ggCAATTAAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGCCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCt                       <  2:1510038/76‑1 (MQ=255)
                                          cgCACTTGAAAAGCTGTCCATCGCCATCTTTGTTTTTTCTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCt                       <  1:3803934/75‑1 (MQ=255)
                                            cAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATcccccta                     <  1:3001055/76‑4 (MQ=255)
                                            cAATTGAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTa                     <  1:1219725/76‑1 (MQ=255)
                                            cAAGTGAAAACCGGCCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTa                     <  2:3207264/76‑1 (MQ=255)
                                             aaTTGAAACCCCGCCCCGCGCCGGCTTTGCTTTTCCGGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTAt                    <  2:306219/76‑1 (MQ=255)
                                               tgaaaaacCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCt                  <  1:5022482/73‑1 (MQ=255)
                                                tGAAAAACGGCCCAGCGCCAGCTTTGTTTTTACTGTCGACCAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCtt                 <  2:5059465/76‑1 (MQ=255)
                                                 gAAAAACTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCttt                <  2:121364/76‑1 (MQ=255)
                                                 aaaacaCTGACCCGCGCACGTTTTGCTTTTACTGTCGACCAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCttt                <  1:4489348/75‑1 (MQ=255)
                                                  aaaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGGCCGCAGAATCCCCTTATCtttc               >  2:5016638/1‑75 (MQ=255)
                                                   aaacCTGTCCGGCGCCACCTTTGCTTTTGCTGCCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAc              <  1:1367469/76‑1 (MQ=255)
                                                   aaaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAc              >  2:3436059/1‑76 (MQ=255)
                                                    aaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAcc             >  1:1256889/1‑76 (MQ=255)
                                                    aaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAcc             >  1:4717790/1‑76 (MQ=255)
                                                    aaaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTAcc             >  1:781521/1‑76 (MQ=255)
                                                     aaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCa            >  2:3564954/1‑76 (MQ=255)
                                                     aaCTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCa            >  1:5406448/1‑76 (MQ=255)
                                                      aCTGCCCAGCGCCAGCTTTGCTTTTACTGTCGACCAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCAt           <  1:3576063/76‑1 (MQ=255)
                                                       cTGACCAGCGCCAGCTTTGCTTTTCCTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATc          >  1:2724818/1‑76 (MQ=255)
                                                       cTGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTCCCATc          >  2:403284/1‑76 (MQ=255)
                                                        tGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCCGCAGAATCCCCTTATCTTTACCATCt         >  2:3942523/1‑76 (MQ=255)
                                                        tGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt         >  1:4277015/1‑76 (MQ=255)
                                                        tGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt         >  1:2646425/1‑76 (MQ=255)
                                                        tGACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCt         >  2:2424770/1‑76 (MQ=255)
                                                         gACCAGCGCCAGCTTTGCTTTTACTGTCGAACATCCGCAGCAATGGGCAGCAGAATCCCCTTATCTATCCCATCTg        >  1:3193110/1‑76 (MQ=255)
                                                         gACCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTg        >  1:4067525/1‑76 (MQ=255)
                                                          cccAGCGCCAGCTTTGCTTTTACGGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTgg       <  2:2853936/75‑1 (MQ=255)
                                                          aCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTgg       >  1:2216006/1‑76 (MQ=255)
                                                          aCCAGCGCCAGCTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTgg       <  1:5036963/76‑1 (MQ=255)
                                                               cgcCAGCTTTGCTTTTACTGCCAAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATg  <  2:1281971/76‑1 (MQ=255)
                                                                    ||                                                                      
TGATGGCGAACGCGTGGTGCACAGCAGCGCCATTGATCATTTGGCAATTGAAAAACTGACCAGCGCCACGTTTGCTTTTACTGTCGAACAGCCGCAGCAATGGTCAGCAGAATCCCCTTATCTTTACCATCTGGTCATG  >  NC_000913/3221029‑3221167

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