Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,980,918 2 bp→CG coding (348‑349/1236 nt) aslB → putative arylsulfatase regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,9180GC99.4% 430.8 / ‑4.9 158W116C (TGG→TGCaslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (81/76);  minor base T (1/0);  total (82/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01
*NC_0009133,980,9190CG99.4% 428.2 / NA 158P117A (CCG→GCG) aslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base G (81/75);  minor base T (0/1);  total (82/76)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

TATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCA  >  NC_000913/3980850‑3980986
                                                                    ||                                                                   
tatGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtc                                                               <  1:2742936/76‑1 (MQ=255)
 atGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtct                                                              <  2:4329527/76‑1 (MQ=255)
  tGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGATGAATGGTGCGCGTTtctc                                                             <  2:5041113/76‑1 (MQ=255)
  tGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                             <  1:1914286/76‑1 (MQ=255)
  tGGTGCGGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                             <  2:1445094/76‑1 (MQ=255)
    gTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTcgc                                                           <  2:1690522/76‑1 (MQ=255)
    gTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTcgc                                                           <  2:4285969/76‑1 (MQ=255)
      gCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCgg                                                         <  1:5156375/76‑1 (MQ=255)
      gCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCgg                                                         >  1:5753024/1‑76 (MQ=255)
       cTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                        >  1:4455123/1‑76 (MQ=255)
       cTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                        <  1:4047138/76‑1 (MQ=255)
       cTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                        <  2:2890416/76‑1 (MQ=255)
       cTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                        <  2:4024482/76‑1 (MQ=255)
       cTGGCAGGAAGATAAGCAACAGCTTCCAGACTAACGGCGTGCTGATGGGTGACGAATGGTGCGCGTTTCTCGCGGa                                                        >  1:2547902/1‑76 (MQ=255)
       atGGCAGGAAGATAAGTAACAGCTTCCAGACTAGCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                        <  1:4277401/75‑1 (MQ=255)
        tGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTTCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAg                                                       <  1:116044/76‑1 (MQ=255)
        tGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAg                                                       <  2:1415717/76‑1 (MQ=255)
        tGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAg                                                       >  2:724280/1‑76 (MQ=255)
         ggCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCgg                                                         >  1:2562041/1‑73 (MQ=255)
         ggCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCgg                                                         <  2:2562041/73‑1 (MQ=255)
         ggCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGc                                                      <  1:2100269/76‑1 (MQ=255)
          gCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGca                                                     <  2:41501/76‑1 (MQ=255)
          gCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGca                                                     >  1:181693/1‑76 (MQ=255)
           cAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcat                                                    >  1:1711568/1‑76 (MQ=255)
           cAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcat                                                    <  2:2076052/76‑1 (MQ=255)
            aGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatc                                                   >  1:4318298/1‑76 (MQ=255)
            aGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatc                                                   >  1:1042562/1‑76 (MQ=255)
              gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatcat                                                 >  1:4993684/1‑76 (MQ=255)
               aaGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAtt                                                >  2:2591992/1‑76 (MQ=255)
               aaGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAtt                                                >  1:4310456/1‑76 (MQ=255)
               aaGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAtt                                                <  1:2345060/76‑1 (MQ=255)
                aGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCCTCAttt                                               >  2:4286720/1‑76 (MQ=255)
                aGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAttt                                               <  1:2681388/76‑1 (MQ=255)
                aGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAttt                                               >  2:3070250/1‑76 (MQ=255)
                aGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAttt                                               >  1:2228270/1‑76 (MQ=255)
                 gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGCGCATCAtttt                                              >  2:2881285/1‑76 (MQ=255)
                 gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAtttt                                              >  2:3954892/1‑76 (MQ=255)
                  aTAAGTACCAGCTTCCAGACTCACGGCGTGCTGCTGGATGCCGCATGGTGCGCGTTTCCCGCGGAGCATCATTTTc                                             >  1:944932/1‑76 (MQ=255)
                  aTAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCCTCAttttt                                             >  2:159637/1‑75 (MQ=255)
                  aTAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTc                                             >  1:915130/1‑76 (MQ=255)
                   tcaGTGACAGCCGCCAGACTAAGTGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCt                                            <  2:676290/74‑1 (MQ=255)
                   tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCt                                            >  2:3359170/1‑76 (MQ=255)
                   tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGATTCTCGCGGAGCATCATATTCt                                            >  1:4681488/1‑76 (MQ=255)
                    aaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCtt                                           <  1:4186703/76‑1 (MQ=255)
                    aaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCtt                                           >  1:5530954/1‑76 (MQ=255)
                     aGTAACAGCTTCCAGACTACCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttg                                          <  1:4926542/76‑1 (MQ=255)
                     aGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttg                                          <  2:2537518/76‑1 (MQ=255)
                      gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgt                                         <  1:843189/76‑1 (MQ=255)
                       tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgtt                                        >  1:3000108/1‑76 (MQ=255)
                        aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttggtg                                       >  1:4083399/1‑76 (MQ=255)
                          cAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGt                                     <  2:1179434/76‑1 (MQ=255)
                          cAGCTTCCAGACTAACGGCGTGCTGCGGGATGCCGAATGGTGCGCGTTTCTAGCGGAGCATCATTTTCTTGTTGGt                                     >  1:4403298/1‑76 (MQ=255)
                           aGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGCTTCTCGCGGAGCATCATTTTCTTGTTGGtt                                    >  2:3113637/1‑76 (MQ=255)
                            gCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGttt                                   <  2:4408160/76‑1 (MQ=255)
                               tCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATc                                <  2:1046644/76‑1 (MQ=255)
                               tCCAGACTAACGGCGTGCGGCTGGATGACGAATGGTGCGCGTTTCTCGCGGTGCATCATTTTCTCGTTGGCTTATc                                >  2:4333699/1‑76 (MQ=255)
                                ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCg                               <  2:1059645/76‑1 (MQ=255)
                                  aGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCt                             >  1:1924667/1‑76 (MQ=255)
                                  aGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCt                             <  1:5633435/76‑1 (MQ=255)
                                  aGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTATATCGCt                             >  1:4380140/1‑76 (MQ=255)
                                    aCTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTgg                           >  1:4208951/1‑76 (MQ=255)
                                    aCTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTgg                           >  2:5848565/1‑76 (MQ=255)
                                      tAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGAt                         <  2:1770648/76‑1 (MQ=255)
                                      tAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGAt                         <  1:2398948/76‑1 (MQ=255)
                                      tAACGGCGTGCTGCTGGATGACGAATGGAGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGag                         >  2:3770682/1‑75 (MQ=255)
                                       aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATg                        >  2:5696120/1‑76 (MQ=255)
                                       aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATg                        <  2:2708406/76‑1 (MQ=255)
                                        aCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCtt                                           <  1:4730353/56‑1 (MQ=255)
                                        aCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCtt                                           >  2:4730352/1‑56 (MQ=255)
                                        aCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATgg                       <  1:5726249/76‑1 (MQ=255)
                                          ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc                     >  1:3198472/1‑76 (MQ=255)
                                          ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc                     >  1:1506418/1‑76 (MQ=255)
                                          ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc                     >  1:2383788/1‑76 (MQ=255)
                                           gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGccc                    >  2:2511490/1‑76 (MQ=255)
                                             gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGc                  >  2:5172481/1‑76 (MQ=255)
                                             gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGATGTATGTCACGc                  >  1:1781559/1‑76 (MQ=255)
                                             atgctgcTGGATGACGAATGGGGCGCGGTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGc                  >  2:1928890/2‑76 (MQ=255)
                                              tgctgcTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGcc                 >  1:3053599/1‑76 (MQ=255)
                                               gctgctGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCt                <  1:1023564/76‑1 (MQ=255)
                                               gctgctGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCt                >  2:3577654/1‑76 (MQ=255)
                                               gctgctGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCt                <  2:3407953/76‑1 (MQ=255)
                                                 tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTga              >  2:45548/1‑76 (MQ=255)
                                                 tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTga              >  1:1409573/1‑76 (MQ=255)
                                                 tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTga              <  2:4084279/76‑1 (MQ=255)
                                                 gtctgTATTACGAATGGTTCTCGTGTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTga              <  2:4060311/74‑1 (MQ=255)
                                                   ctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTgaga            >  2:1413341/1‑76 (MQ=255)
                                                    tgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGAt           >  1:3344733/1‑76 (MQ=255)
                                                      gATGACGAATGGTGCGCGTTTCTCGCGGCGCATCCTTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGGGCTcc         >  2:5334776/1‑76 (MQ=255)
                                                      gATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc         >  2:5282026/1‑76 (MQ=255)
                                                       aTGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCca        <  2:3344733/76‑1 (MQ=255)
                                                        tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTCTCGCTGGATGGCGCGCCTGAGATCcac       <  1:5282225/76‑1 (MQ=255)
                                                         gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCcaca      <  2:1160819/76‑1 (MQ=255)
                                                         gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCcaca      <  1:1223470/76‑1 (MQ=255)
                                                         gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCcaca      <  1:37212/76‑1 (MQ=255)
                                                          aCGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaa     >  1:4834866/1‑76 (MQ=255)
                                                            gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatc   >  1:246309/1‑76 (MQ=255)
                                                            gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatc   <  1:2028187/76‑1 (MQ=255)
                                                             gaTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatca  <  2:544642/75‑1 (MQ=255)
                                                             aaTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatca  >  1:2348674/1‑76 (MQ=255)
                                                             aaTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatca  <  2:911380/76‑1 (MQ=255)
                                                                    ||                                                                   
TATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCA  >  NC_000913/3980850‑3980986

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